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The grant period was for two years: from January 1994 to 
December 1995. The grant, basically, supported the salary o a 
senior research associate and provided some travel 
attend conference and 25% salary of a graduate student for one 
year. No salary money was allocated to the principa 
investigator (PI). Dr. Ram Tripathi was selected to be the 
senior research associate. Although the grant was initially 
approved for three years subject to availability ° f : ^ d 

discontinued after two years because of (a) a lack of funds and 
(b) reorganization within NASA, which eliminated the branch and 

the research program. 


Despite the short duration, 
productive. There have been five 
as shown in attachment 1 . 


the grant has been very 
following short communications: 



1 Bose condensation of nuclei in heavy-ionco,,,smn fimte b * L W 

Townsend and R.K. Tripathi, Bull. Am. Phys. Soc. 39, 1395 ( ) 

2. Observation ofBose condensation of nuclei in hea^-ion collision by R K Tripathi and 
L.W. Townsend, Bull. Am. Phys. Soc. 39, 1395 (19 ). 

3 Evaporation of strange panicles in anti-proton nucteus ^.ions by F B Malik, L.W. 
Townsend and R K Tripathi, Bull. Am Phys.Soc. 39, 1426 (1994). 

4 . Nuclear fragmentation of heavy-ion by protons, by L.W. Townsend, and R.K. Tnpatht, 
Bull. Am.Phys. Soc. 39, 1426 1(1994). 

1632(1995). 

In addition, eleven full length ariicles have “Pf" ed “ ° r ?*** ° r bee " 

submitted for publication These are listed below and attached in Attachment 


1 . 


3. 


5. 


6 . 


7 . 


Accurate universal parameterization of absorption cross s^ion by RK Tripathi, F A 
Cucinotta and J.W. Wilson, accepted for publication in N.I.M-B (1 )• 

Universal parameterization of absorption cross sections by R.K. , Tripathi, F A. Cucinot 
and J.W. Wilson, NASA Report: N AS A-TP- 17580, (1996). 

Emission of light fragments as an indicator ofequilibratedpopulaUonsin nuclear collisions 
by R.K, Tripathi and L.W. Townsend (submitted for publication 1996). 

Transport of light ions in matter by J.W. Wilson, F. A. Cucinotta, H. Tai, J.L. Shinn, S Y 
Chun, R.K. Tripathi and L. Sihver (submitted for publication). 

Bose-Einstein condensation of nuclei by L W. Townsend and R K Tripathi (Accepted for 
publication in Condensed Matter Theories Vol. 1 1, 199 ). 

Theoretical model of HZE particle fragmentation by hydrogen targets _by LW. Townsend, 
F A. Cucinotta, R. Bagga and R. K. Tripath,, Adv. Space Research 12, 109 (1996). 

Liquid drop model considerations in HZE particle fragmentation by hydrogen by L* 
Townsend R.K Tripathi, F A. Cucinotta, and R. Bagga (submitted for pubhcat.on, 1996). 

NUCFRGZ An evaluation of the semi-empirical nuclear fragmentation database by J.W. 
^fon R K Tripathi, F A Cucinotta, J.L. Shmn F F Badavi, S.YXhun, J.W. Norbuot, 
C. J. Zeitlin, H. Heilbronn and J. Miller. NASA Technical Report 3533. 



9 Atomic electron correlation and particle and anti-particle induced single and double- 
ionization by T. Das and F. B. Malik (Accepted for publication in Condensed Matter 
Theories, vol 11 (1996)). 

10. On the ionization of He by proton and anti-proton: I. Single ionization by T Das and F. B 
Malik (submitted for publication, 1996). 

11 On the ionization of He by proton and anti-proton: II. Double ionization by T. Das and F. 
B. Malik (submitted for publication, 1996). 

The research done under the grant is recognized to be important and Dr. Townsend, one 
of R K Tripathi’s collaborators and F. B. Malik were invited to present their research done under 
the grant in plenary sessions at XIX International Workshop on Condensed Matter Theories held 
in Caracas, Venezuela in June 1995. 
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THURSDAY AFTERNOON 


of grid points in proportion to dependent variable. 
The three-dimensional results presented are for ult- 
ra relativistic condition, with a steep density grad- 
ient in initial density distribution. An elliptic gen- 
eration system was solved by SOR iteration techn- 
ique, to control grid skewness, by smoothing grid. 
Results of this simulation were compared with the 
available simulations. 


)ense Matter G. Q. LI, B. A. LI and C. M. KO, 
i^-Kaon flow in heavy-ion collisions is studied in the 


relativistic transport model and is found to be sensitive to the 
kaon mean-field potential in dense matter. With an attractive kaon 
mean-field potential, kaon flow is in the same direction as the nu- 
cleon flow, while with a repulsive mean- field potential, it is opposite 
to that of nucleons. Kaon flow in heavy-ion collisions is thus a useful 
observable for studying kaon properties^ dense matter. 

•Supported by NSF Grant No. 9212209 and the Robert A. Weld 
Foundation Grant No. A- 11 10 1 


AE7 Bose Condensation of N uclei in Heavy Ion 
Collisions: Finite System Effects , h. W. TOWNSEND, 
NASA Langlev Research Center, and R. K. TRIPATHI, 
foil them Illinois University *. —In previous work the 
possibility of Bose condensation of nuclei in heavy ion 
collisions was demonstrated 1 by using quantum statistic* 
to model an equilibrated system consisting of a finite 
number of nucleons. The existence of a concomitant toi al 
phase transition, however, was affected by the relatively 
small numbers of nucleons involved in the collision 
(~ 80). In this work we investigate the effects of finite 
nucleon numbers on the observation of a possible total 
phase transition. 

•Supported by NASA Grant. I 

’R. K. Tripathi and L. W. Townsend, Phys. Rev. C (1994, ih 
press) 1 


14:54 

s S-Wave x-x Correlations in Cold Nuclear Matter. 

R. RAPP, Z. AOUISSAT, J. WAMB ACH, G. CHANFRAY 1 , P. 
SCHUCK 2 , KFA Julich. *IPN Lyon, 2 ISN Grenoble. — Start- 
ing from a standard treatment of the single pion selfenergy in 
terms of particle-hole and A-hole Lindhard functions we exam- 
ine x-x s-wave correlations in cold nuclear matter. When calcu- 
lating the in-medium x-x T-matrix we take into account the full 
off-shell properties of the pion selfenergy. The latter enters the 
scattering equation of Li ppmann- Schwinger type through the 
uncorrelated 2x propagator of the intermediate state. Using a 
well tested phenomenological x-x interaction (the Julich meson 
exchange model) we find a x-x pairing instability in slightly 
compressed nuclear matter (/> ® l-3po). This model, however, 
does not respect constraints from chiral symmetry which are 
known to govern the near threshold behaviour of x-x scattering. 
Additionally imposing chiral constraints on the Julich model in 
fact inhibits the instability leading to more moderate effects in 
the subthreshold region of the in-medium scattering amplitude. 


1 5*06 

AE 9 iMoiitron- N ^itron Correlation .Funct i ons. C. M. 
Mader R. S. Bennink and P. A. DeYoung, Hope College - The 
structure of two-proton correlation functions measured at small rel- 
ative angles is dominated by the 2 He resonance peak. The heignt of 
this peak is then related to the size of the emitting system. In the 
two-neutron system, the Coulomb interaction is jio longer a factor, 
and thus only the nuclear interaction and Pauli exclusion princi- 
ple wi.l play a role in the correlation function. It is expected that 
without a 2-neutron resonance, the correlation peak will shift to 
zero relative momentum and the source size will be determined by 
the width of the peak. We use the Koonin formalism to calculate 
the two-neutron correlation function. We explore the effects of the 
short-ranged nuclear potential as well as the parameterization o 
the neutron source and compare with recent experimental results 
for 210 MeV 16 0 + 77 A1 reactions. 


13,10 , 

AH io Observation of Bose Conde nsation of Nuclei yn 

Hpaw Inn Collisions: R. K. TRIPATHI, Southern Illinois 
University *, and L. W. TOWNSEND, NASA Langley 
Research Center -We have demonstrated 1 the presence of 
Bose condensed nuclei in heavy ion collisions. Here, we 
present results identifying the optimum conditions for 
their observations and focus on the ways of experimental 
verification of this new, exciting phenomenon. 
•Supported by NASA Grant. 

Ir. ic Tripathi and L. W. Townsend, Phys. Rev. C 50, R1 
(1994) 


15:30 v . 

AE 11 Self-consistent momentum-dependence m hot and dgps e 
nuclear matter, * VINOD K. MISHRA and G. FAI, 

— The momentum dependence of the nucleon-nucieon interaction 
and of the mean field has important consequences for the 
properties of nuclear matterfl] and for collision simulations^]. In 
particular, momentum-dependent interactions (such as the 
Momentum-Dependent Yukawa Interaction or MDYI[3]) lead to 
a stiffer equation of state (EOS) for nuclear matter than the 
corresponding momentum-independent ones. We have earned out 
in the past couple of years a study of local and global observables 
of hot and dense nuclear matter with the MDYl in successively 
better approximations. To complete this program we now examine 
the same observables in a self-consistent Hartree approximation 
which moves the project from a largely analytical domain to a 
purely numerical study. We use the earlier analytical 
approximations as input to the self-consistent calculation. We 
present results on the optical potential, the effective-mass, pressure 
(EOS), and the incompressibility for the Maxwell-Boitzmann 
distribution and compare them with the momentum-independent 
results and results of the earlier analytical approximations 
• Work supported by DOE Grant bE-FG02-86ER40251. 

1. V.K. Mishra, et al, PRC47, 1519 (1993). 

2. C, Gale, G.F. Bertsch and S. Das Gupta, PRQ5, 1666 (1987). 

3. G.F. Bertsch and S. Das Gupta, Phys. Rep. I£Q, 189 (1988). 


15:42 s 

AE 12 __ .. . 

Toroidal Structures in BUU-Predictio ns for Ar + — 5c_ 

at E/A=80 MeV viewed by Proton-Proton Correlations ' 
G J KUNDE, ' S.J. GAFF, C.K. GELBKE, T. GLASMACHEK, 
D.O. HANDZY, L. MARTIN, W. BAUER, F.C.DAFF1N 
NSCL, Michigan State University 

For central (b=0 £m) collisions BUU- predictions suggest that a 


Sm- 
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In first order the use of a free NN t-matrix accounts for the short-range 
interaction between the projectile and a target nucleon. The effect of the 
nuclear medium can be interpreted as higher order correction and arises 
from the difference between the free NN propagator and the propagator 
in the nuclear medium. A formal framework consistent with a multiple 
scattering expansion has been constructed to include these contributions 
by using nuclear mean field potentials, which when applied to p-nucleus 
scattering gives a significant improvement of the elastic scattering ob- 
servables at energies less than 200 MeVX The influence of these medium 
modifications on the free NN t-matrix will be discussed in detail for dif- 
ferent nuclei, together with its energy dependence between 100 and 800 
MeV and its influence on the off-shell structure of the NN t-matrix. 

* Supported by the U.S. Department of Energy, 
t C.R. Chinn, Ch. Elster, R.M. Thaler, Phys. Rev. C48 (1993), 2956. 


&hr6dinger equation has been extended to scattering involving Klein- 
Gordon equation with a view of determining complex potentials for 
scattering of spinless particle by a spinless target at relativistic 
energies. The formalism replaces the differential equation by 
difference equation and the inversion results into solving a continued 
fraction equation and provides number of points of the potential equal 
to the number of partial waves. In case the number of partial waves 
is large i.e., about 10 or more and the knowledge of logarithmic 
derivative is accurate, the potential is determined very well. A few 
model for determining potentials from phase shifts as well as 
logarithmic derivatives will be presented with a view of applying the 
method for the pion -nucleus case. 

[1] M.M. Alam and F.B. Malik, Phys. Lett. B 222, 14 (1990) and 
Nucl. Phys. A 524, 88 (1991) 


14:30 

EBI Carnality with Noncauaal Potentials ip the Many-Body 
Approach to Elastic Nuclear Scattering. * V.A. MADSEN, Ore- 
gon State University*- It has been shown by Bell and Squires 
[1] that the self-energy operator of the one-particle time or- 
dered Green’s function can be interpreted as the optical po- 
tential for nucleon-nucleus scat; -ring. However, Mahaux and 
Sartor [2] have demonstrated that this time nonlocal poten- I 
tial is noncausal in the sense that ^ 0 for f > t, / 

where t is the current time and if, the nonlocal time. We / 
shall demonstrate, mak: " use of the retarded Green’s func- / 
tion, that the time'depr Lippmann- Schwinger equation ( 

for the elastic scatter^ ..evertheiess carnal V 

* Research supported in part by the U.S. Department of En- 
ergy under Contract FG96-8SER4C283 at Oregon S*ate Uni- 
versity. 

1. J.S. Beil and E.J. Squires. Phys. Rev. Lett. S r 96 (1959). 

2. C. Mahaux and R.. Sartor, Nuclear Physica, A530 , 303 i 

(1991). / 


14:42 

EE 7 Validity of Sinha'* local apyroximati on to the j 
exchange nonlocnl j tr. D. LUX ACZZ? ai*d G.H. KAVCTSCKCL » f 

Urriv, g£ Connecticu t. Several method* for replacing J 

exchange nonlocalities with local equivalent potential* 
exiat m the literature. In testing the validity of / 
S inha's method [1], Geonriev et al [2] suggested an 
improved version which, includes the imaginary part of J 
the 1oc5a1 manenuja. In the present study we apply theiir 

[2] version to the case of n- O scattering at several! 
energies frora 20 to 100 h>eV. In a prior investigation ^ 

[3) , we obtained the n- lr O nonlocal potential by the 
folding of a c^piex g-raatrix and now use the 
corresponding pnase equivalent inversion potential [3] 
for car. -ri son with Sinha’s approx inat ion. Results will 
be presented and discussed. 

[1] B. Sinha, Phys. Reports 20 C, 1 (1975). 

[2] B. 0. Georg iev and R. S. Mackintosh, Fhys. Lett. B 

73, 250 (1978). 

[3] G. H. Rawitsche^, D. Lukaszek, R. S. Mackintosh, 
and S. G. Cooper, Fnys. Rev. C 49, 1621 (1994). 


An Inverse Scaring Theory at a meo fcnergv to: kicip- 
cuaucr. . 2. F. SHHHADBH and F. Bary MALIK, Smhzni 


lllmoi? Univ. a: Cari>ondr- )e.— 1 The inverse scattering formalism of 
Alam ana Iwik [Ij at a fixed energy for potential scattenng using 


15:06 

EB 9 Evaporation of Strange Particles in Antiproton I 
Nucleus Reactions : F. B. MALIK, Southern Illinois I 
University* . L. W. TOWNSEND, NASA Langlev Research 
Center, and R. K. TRIPATHI, Southern Illinois / 

University * -We have investigated the production of / 
neutral strange particles in antiproton - nucleus collisions, 
and found that an evaporation - like model provides I 
better explanation for their emission. New systematic* 
about their production will be discussed. Detailed / 
comparisons will be made with the available experinjental 
data. / 

S *$tr*pQrted by NASA Grant. 


15:18 

EB 10 Nuclear Fragmentation of Heavy Ions bv Protons 
R. BAGGA, Old Dominion University. L.W. TOWNSEND, 
NASA Langlev Research Center. R. K. TRIPATHI, Southern 
THinnis University . — Optical Model methods for calculating 
cross sections for the breakup of heavy nuclei by protons are 
presented. The fragmentation cross sections are calculated 
with a modified abrasion-ablation collision formalism where 
the abrasion stage is treated as a quantum-mechanical 
knockout process and the ablation stage incorporates 
excitation energy from particle-hole and frictional-spectator- j 
interaction (FSI) processes. Comparisons with recently / 
published cross section results will be presented. j 


SESSION EC: HEAVY-ION REACTIONS III: 
CORRELATIONS AND COLLECTIVE PHENOMENA 
Saturday afternoon, 29 October 1994 
Gen Hooker's Redoubt at 13:30 
T. C. Sangster, presiding 


13:30 

EC 1 Systematic Study of Hid 
with an Exclusive lmr>act-pa r ame 

D. R. BOWMAN, D. FOX, 

E. HAGBERG, AECL Rescarc 
BEAULIEU, R. LAFOREST, 
Univcrsite Lava) and T. LI, A. V.A 


tclustv? Impact-pa r ameter Trigger G.C. BALL, D. HORN, 
>WMAN, D. FOX, A. GALINDO-URIB ARRI and 
£RG, AECL Research. Chalk River Laboratories: L. 
U, R. LAFOREST, V. LAROCHELLE and R. ROY, 
Lava) and T. LI, A VANDER MOLEN, G D WESTFALL, 


J.S. WINFIELD, J. YEE AND S.J. YENELLO, Michigan State 
University — Svmmetric heavy -ion collisions were studied at TASCC . 
and NSCL with 22-80A MeV beams of “Ne 3U7 C1 and ^Ar, The 4ic 
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14:42 
EA 3 

Structure of Neutron Rich Nuclei.* 

I. Y. LEE, Lawrence Berkeley National Laboratory. 



‘ V-- 'yf. 


Neutron rich nuclei are of particular interest since they might reveal new aspects of nuclear structure 
associated with an excess of neutrons, such as new region of deformation, shell effects and modes of 
excitation. Deep-Inelastic reactions have been shown to produce neutron-rich nuclei with a high multi- 
plicity of gamma-rays. However, since these reactions produce many final nuclei. It has been very difficult 
to study them until now, due to the lack of sensitivity of the available gamma-ray detector arrays. We have 
carried out gamma-spectroscopy studies of neutron rich nuclei using the reaction 4S Ca+ lv6 Yb. A Si-strip 
detector was used to detect the projecrile-iike fragments and the coincident gamma rays were detected in 
the Gammasphere. States with spin as high as 20 in neutron-rich nuclei such ' * ; Yb and '*Tm were 
observed from 2- and 3-fold gamma-ray coincidence data. We will discuss results cn the variation of 


nuclear structure over a wide range of neutron numbers. For example, the study of pairing strength as a 
function of spin, and the variation of the interaction strength of the fust backbendingf , 

r v . V; i ■ : , v ‘ .i- ; 

"’♦'Work supported by the'TJlS. DOE under contract number DE-ACD3-76SF0098 \ 
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15:18 . * . 

EA 4 New Insights from Studies of Spontaneous Mission: Yields and Neutron Multiplicities. Cold Fiss isn. ind Structure 

of Ncutron-Rica Nuclei/ ‘ ■ 

J. H. HAMILTON, Vanderbilt University. 


An overview of many of the new insights extracted from v,y and y-y-y coincidence studies of spontaneous fission of 
^Cf and :42 Pu with large detector arrays at ORNL and Gammasphere will be presented. These include direct 
measurements of yields and neutron multiplicities from 0 to 10 v for five correlated pairs, identification of a number of 
cold fission (cluster radioacuvity) channels, identification of new structures and behaviors in many neutron- rich nuclei 
from w Sr to 160 Sm. From the yields of the correlated pairs, a new, second mode of SF is observed in the Mo-Ba pairs 
with much lower total kinetic energy. This mode goes via a hyperdeformed shape (3:1 axis ratio) for 144 3a, l45 3a, or 
146 Ba. Calculations indicate a third minima in 232 Cf at - 0.7 and p 2 ~ 0.9. The zero neutron cluster radioactivity 
yields for odd-odd Mo-Ba pairs are about four times larger than for the even-even ones in agreement with recent cluster 
radioactivity predictions. Identical bands are observed in 98 ' 100 Sr, l08,n0 Rn, 144 * 146 Ba, and A32 * :54 Nd nuclei and new types 
of identical bands in i36 * i58 * i60 s m nuclei, and the first identical octupole bands. Octupoie deformation is seen in N = 86, 
l42 Ba, 344 Ce and the new high spin octupole states to 19* exhibit the first backbending and quenching of the static octupoie 
strength as predicted by theory. Selected examples of the new physics being seen will be presented. 

*Worked carried out by Vanderbilt, Dubna, ORNL, LBL, INEL, Warsaw U., U. Tennessee, U. Pittsburgh, U. 
Califomia/Berkeiey, Inst. Phys. Bratislava collaboration. Work supported by U.S. Dept, of Energy, Russian Federal 
Foundation, NSF, SASc, Polish Comm. Sc. Res. 


SESSION EB: HEAVY IONS IV: NUCLEON-NUCLEON 
CORRELATIONS AND INTERMEDIATE ENERGY 
Saturday afternoon, 28 October 1995 
Oak Room at 13:30 
T. C. Sangster, presiding 

13:30 

EB 1 Proton-Proton Correlation Functions for 18 0 + 1#T Au colli- 
sions at E/A = 200 MeV S.J. GAFF, D.O. HANDZY, W. BAUER. F.C. 
DAFFIN, J.D. DINIUS, C.K. GELBKE, T. GLASMACHER, E. GUALTIERI, 
S. HANNUSCHKE, M.J. HUANG, W.C. HSI. G.J. KUNDE, M.A. LISA, W.J. 
LLOPE, W.G. LYNCH, L. MARTIN, C.P. MONTOYA, R. PAK, L. PHAIR, S. 
PRATT, C. SCHWARZ, N. STONE, M.B. TSANG, A.M. VANDER MOLEN, 
G.D. WESTFALL, J. YEE, AND S.J. YENNELLO, National Superconducting 
Cyclotron Laboratory, \fichigan State University — Using a high-rwoiution 
hodoscope, coincident protons were measured from the reaction 18 0 + 1,7 Au 
at E/A = 200 MeV. The experiment was preformed at the National Super- 
conducting Cyclotron Laboratory. The results to be presented include the 
correlation functions, boih integrated and gated on total momentum, and the 
proton energy spectra. These experimental results will be compared with BUU 
calculations. Consistent with other measurements, discrepenciea with BUU 
predictions will be shown. 1 

Supported in part by the National Science Foundation. 


13:42 

EB 2 Proton-proton correlations in peripheral colliiions. 

H.M.Xu, G.K. Ajupova, C.A. Gagharoi, Y.W. Lui, A.A. Reid and 
R.P. Schmitt Texas AfcM U." — We have measured proton-proton 
correlation functions in 14 N induced reactions on 37 AI and n *Sn 
targets at average angles = 10° and 9 a „ — 25° with the Texas 
A&M University Proton Spectrometer. The spectrometer data will 
be compared to earlier multi-dectector data and the dependence of 
proton-proton correlations on the emision angles and the momenta 
of the proton pairs will be discussed. 

"Work supported in part by the Department of Energy under Grant 
No. DE-FG03-93ER40773 and by the Robert A. Weich Foundation. 

13:54 
EB 3 

One Nucleon Removal Cross Sectio ns in Light and 
Medium Nuclei: R. K. TRIPATHI, Southern Illinois 

University *. J. W. WILSON, F. A. CUCINOT7A, J. L. SHINN, 
NASA Langiev Research Center. F. F. BADAVI, Christopher 
Newport University , and S. Y. CHUN, Old Dominion 


^ Vol. 40, No. 10 (1995) ^ ^ /}, 1632 
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\ /piversitv —Using a semiempirical nuclear fragmentation 
model developed 1 at NASA Langley Research Center 
(NUCFRG2), we analyze one nucleon removal cross sections 
in terms of nuclear and coulomb cross sections. We find 
there is a need for structure effects in light and medium 
nuclei. The interesting cases of one neutron and one proton 
removal cross sections in 40 Ar and 5b Fe will be discussed. 
Within a simple single particle shell model picture, both 
these nuclei can be described by 2p-2h states with respect to sd 
and fp shells, respectively. Experimentally, one neutron 
removal cross sections for these nuclei are approximately 
three times larger than one proton removal cross sections 
reflecting to the neutron shells to be "farther ' than proton 
shells. 

•Supported by NASA Grant. 

ij. W. Wilson, et al., NASA TP-3533, 1995 

14:06 

EB4 Nuclear Dissipation from Measurements of Light 
Charged Particle Emission in Fusion — Fission Reactions 
R.C. LEMMON, NSCL, Michigan State U. , D.J. HINDE, J.R. 
LEIGH, C.R. MORTON, J.O. NEWTON, M. DASGUPTA, 
Australian National U. . J.P. LESTONE, U. of Washington . 

At moderate excitation energies, the dynamics of nuclear fission are 
determined by the dissipative coupling of collective and internal sin- 
gle -particle degrees of freedom. One of the foremost problems in nu- 
clear physics is understanding the mechanism and magnitude of this 
dissipation. The alpha particle and proton multiplicities from the 
reactions 19 F+ lM Tm, M Si+ l60 Gd and ^S+^Sm have been mea- 
sured. Pre- and post-saddle timescales for the fission of 188 Pt were 
then extracted from these multiplicities using a statistical model 
modified to include fission delay times and a Kramers factor [1). 
The timescales have now been interpreted within the framework of 
a Langevin collective transport model to give values of the dissipa- 
tion as a function of deformation. These results will be presented 
and comparisons made with different microscopic models of nuclear 
dissipation [2,3]. 

[1] J.P. Lestone, Phya. Rev. Lett. TO, (1991) 2245 

[2] J.R. Nix and A.J. Sierk, Proc. International School- Seminar 
on Heavy Ion Physics, Dubna, 1986, Joint Institute for Nuclear 
Research Report No. JINR-D7-87-68 (1987) 

[3] P2 FVobrich, LI. Gontchar and N1D. Maviitov, Nuclear Physics 

A556 (1993) 281 ** ^ 

14:18 

EB 5 Characterizing the Two S t a ye Reaction Dynamics in I 
AGeV Projectile Fragmentatio n of Gold Nuclei XA. HAUGERL S. 
ALB ERGO 2 , F. BIESER 3 , FJ\ BRADY 4 , Z. CACOAT DA. 
CEBRA 4 , AJ3. CHACON 3 , J.L. CHANCE 4 , Y. CHOI 1 . S. COSTA‘. 
LB. ELLIOTT 1 , MX. GILKES 8 , A.S. HIRSCH 1 , E.L. HJORT 1 , A. 
INSOLIA 2 , M. JUSTICE 6 , D. KEANE 6 , J.C. KITNER 4 , V. 
LINDENSTRUTH 7 , M. A. USA 3 , U. LYNEN 7 , H.S. MATIS 3 , M. 
MCMAHAN 3 , C MCFARLAND 3 , W.FJ. MUELLER 7 , D.L. 
OLSON 3 , M.D. PARTLAN 4 , N.T. PORILE 1 . R. POTENZA 2 , G. 
RA1 3 . J. RASMUSSEN 3 , H.G. RITTER 3 , J. ROMANSKI 2 . J.L 
ROMERO 4 , G.V. RUSSO 2 , H. SANN 7 , RJ>. SCHARENBERG 1 , A. 
SCOTT 6 , Y. SHAO 6 , B. SRJVASTAVA 1 , TJ.M. SYMONS 3 , MX. 
TINCKNELX 1 , C. TUVE 2 , S. WANG*. P. WARREN 1 . D. 
WEERASUNDARA 6 , H.H. WIEMAN 3 , K.X WOLF 5 , ! Purduc IL 
2 fNFN. Catania . 3 Lawrcncc Berkeley Lab. . 4 UC - Qayis , 5 T?xas A&M. 
LL ° Kcnt Smc XL ' SSL DarmstatiL. 8 SUNY Am Bnrck. - By 
studying protons from the multifrag mentation of lAGeV Au + C 


collisions, we can distinguish two stages in the reaction process. A 
method for separating these two components will be presented. 

14:30 

EB 6 Cavitation and Penetration in Central Collisions Induced 
bv Light Ions . G. WANG, K. KWIATKOWSKI, and V. E. VIOLA, 
Indiana University ; and W. BAUER and P. DANIELEWICZ, 
Michigan State University . Two different BUU calculations are 
employed to examine the dynamics of central collisions induced 
by light-ion projectiles. For projectile energies above about one 
GeV incident on heavy target nuclei, a region of low density 
develops in the core of the nucleus at times of the order of 30 
fm/c. The simulations predict penetration of the target by the 
projectile momenium front above energies of 4-6 GeV, leaning to 
a saturation of deposition energy. These results are examined 
in the context of marked cnanges in reaction observables for 
such collisions. 


EB 7 One Ctgn Pro duction of * Alpha Panicles bv l2 C + 1I£ 

->6 4 He. E. NORBECX, Y.W. CHENG. L.3. YANG, U. of Iowa. 
F.D. INGRAM, Rock Valiev College. According to the Bom 
approximation, the cross section for 1 ^C + *^C *■ 6 * He is 

proportional to p(p‘]|(v r |V|vjf in which n = 14 is the number of 
independent parameters, and V is the potential acting between a 
clusters in the beam with the a clusters in the target. y f is the final 
state with each cluster expressed as a plane wave, y , is the initial state 
consisting of a cluster 12 C wave functions times plane waves for the 
relative motion. The phase space factor, p, is slowly varying. The 
matrix element, (v f (VtVi), is large if the momentum of one or more 
clusters are spectators which have the same momentum after, as before 
the reaction. If two a particles in the target and two in the beam are 
spectators, the threshold for the reaction is 88 Mev, and the sum of the 
energies of the two remaining a particles is 1/3 of the beam energy less 
about 14 Mev binding energy. Such peaks in the sum of the energies of 
two a particles have been observed for 12 C energies from 660 to 1860 
Mev. 
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14:54 *.* .;* . .• .. W: . - . . . " - ; " •“ * 

EB8 ..x ■ - •< ♦!-.> - v 

Quantum field description of high energy multiple sequential 
collisions in dense media . Jicai Pan, Deptr of Physics, McGill 
University, Montreal, Canada . High energy sequential collisions in 
dense media have, so far, been discussed in the framework of classi- 
cal kinetic theory. In this picture, quantum fluctuation is not taken 
into' account.. As a result, inteference between sequential collisions r 
is usually overestimated. We propose here generalized Lagrangians 
that consist of fields with different initial space-time coordinates to 
develop a full quantum description of sequential collisions in dense 
electromagnetic and strong interacting media. Unlike Glauber the- 
ory that is a quantum mechanical description of elastic and quasi- 
elastic hadron-nucleus scatterings, the present formulation provides 
a unified quantum field description of particle production, radia- 
tion, scattering, and particularly the interference in high energy se- 
quential collisions in dense matter that may be produced in high en- 
ergy nuclear collisions. The destructive interference of soft dilepton 
and photon productions, known as Landau- Permoranchuk-Migdal 
effect, in both hadronic gas and quark-gluon plasma is weaker than 
that obtained in classical kinetic theory. 


Voi. 40, No. 10(1995) 


1633 



Attachment 2 
Full Length Articles 



Accurate Universal 


Parameterization of Absorption Cross Sections 


R. K. Tripathi 

Southern Illinois University, Carbondale, Illinois 62901 
Hampton University, Hampton, VA 23668 

Francis A. Cucinotta and John W. Wilson 
NASA Langley Research Center, Hampton, VA 23681 


Abstract 

We present Q simple universal parameterization o. lolal reaction cross sections 
tor any syslem of colliding nuclei valid lor lire entire energy range from a few A 
MeV to a few A GeV. Tire universal picture presented here treats proton-nucleus 
collision as o special case of nucleus-nucleus collision, where the projectile has 
charge and mass number of one. The parameters are associated with the physics 
of the collision syslem. In general lerms Coulomb interaction modifies cross 
sections al lower energies and lire effects of Pauli blocking is imporlanl al higher 
energies. The agreement between Ihe calculated and experimental data is 
better than all earlier published results. 



1 .Introduction 


The Iranspor, alien ol energelic ions in bu, k mglier is of direct inieres, in 

several areas ,1, including shielding agains, ions originaiing from eiiher space 

terrestrial accelerators, cosmic ray propagation studies in galactic 

medium or radiobioiogicai effecls resting „ 0 m ihe war. place or Cinico, 

carcinogenesis, terrestrial radiation therapy. and radiobiological 

research; know, edge o, ihe beam composiiian and inieraclions is necessary ,o 

properly evaluale Ihe elfecls on human and animal tissues. Far Ihe proper 

radiation exposures both reliable transport codes and accurate 
input parameters are needed. 

One such imporian, inpu, is ihe ioia, reaciian cross sec, ion. deiined as ihe 
lolal mrnus Ihe elaslic cross seclions ior Iwo colliding ions; 


'el 


ViBW ^ i,S imPOdanCe ' he tetaI -Han has been exiensiveiy 

died both theoretically |l-14) and experlmenlally [15-20] for Ihe past (ive 
decades. A db, ailed its, o, references is give Ref,,. I3 . , 4J . Empilical 

prescriptions have been deveioped ,2-,. ,0. , ,. ,3, tar ihe faia, reaciian cross 

-Cions wo, k ing in various energy ranges and combine, ion of in, eroding ions The 

p-esen, mode, wor k s in o„ energy ranges for any combination of in, eroding ions 

including proion-nucieus collisions and is more accurate ihan eariier reported 
empirical models. 



2. Model Description 

The present model is on updaled and revised version of the empirical 
model developed al NASA Langley Research Center and reported earlier 
[tO). Most of the empirical models approximate total reaction cross section of 


Bradt-Pelers form: 


- 8)2 


where r 0 is energy independent and 6 is either energy-independent or energy 
dependent parameter, and A, and A, are the profile and large! mass 
numbers, respectively. This form of parameterization works nicely lor higher 
energies. However, for lower energies Coulomb interaction becomes important 
and modifies reaction cross sections significantly. In addition. Strong absorption 
models suggest energy dependence of the interaction radius. Incorpora ling these 
effects, and other effects discussed latte, in lire text, we propose the following 
form for the reaction cross section: 

- -£> 

We notice lhal Coulomb interaction, where r 0 = 1.1 lm. and Ecm is inMeV. 
modifies cross sections al lower energies and gels less imporlant as the energy 
increases (typically after several tens at A MeV|. In Eq. (3) B is the energy 
dependent Coulomb interaction barrier (right hand (actor in Eq. 3). and is given 
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( 4 ) 


B = 


1 . 4 4 ZpZj’ 

R 


where, 


R = r P + r T + 


1 . 2 (Ap /3 + 




(5) 


wilh ( i - PJ). 

r, = 1 • 29 (r t ) M ,6) 

There is energy dependence in Ihe recclion cross seclion at 
intermediate and higher energies mainly due to two effects - transparency and 
Pauli blocking. This is taken into account in 8 E . which is given by. 

Q = 1.85 S + 0 . 16S/EJ, " C-£ (7) 

+ 0.91 {A t ~ 2Z t ) Zpl { A-pAp ) 

where S is the mass asymmetry term and is given by, 

Ap n Ar /y ( 8 > 

" W 1 + tf'* 

and is related to Ihe volume overlap ol Ihe collision system. The last term on Ihe 
right hand side of Eq. (71 accounts lor Ihe isotope dependence of Ihe reaction 
cross seclion. The term C E is related to Ihe transparency and Pauli blocking and 

is given by. 


4 


( 9 ) 


C E = D (1-exp ( -£/40) ) -0 . 292exp ( -E/792) 
x cos (0 .229E 0 ' 453 ) 

Here D is related lo Ihe densily dependence of Ihe colliding system scaled wilh 
respect lo the density of C + C system, ie: 


D = 1.75 


P Ap + Pa, 


( 10 ) 


P A C * P A c 

The densily of a nucleus is calculated in Ihe hard sphere model 1241, and 
for a nucleus of mass number A, is given by. 


Pa, 


471 3 

' f I 


an 


3 1 


where Ihe radius of the nucleus r, is defined in Eq. (6| wilh rooi-mean-square 
radius, |r,), m „ oblained directly from experiment [25]. There is inleresling physics 
associated wilh constant D. This in effecl simulates Ihe modifications of the 
reaclion cross sections due lo Pauli blocking. This effect is new and has not been 
taken info qccounl in other empirical calculations. This helps presenl a universal 

picture of the reaction cross sections. 

At lower energies (below several lens of A MeV) where ihe overlap of 
interacting nuclei is small (and where Coulomb interaction modifies the reaction 
cross sections significantly) Ihe modificalions of Ihe cross sections due lo Pauli 
blocking are small, and gradually play an inaeasing role as Ihe energy increases. 


5 



since this leads to higher densities where Pauli blocking gets increasingly 
important. Interestingly enough for the prolon-nucleus case.since there is not 
much compression effect, a single constant value of D=2.05 gives very good 
results for all proton-nucleus collisions. For alpha - nucleus collisions, where there 
is a little compression, the best value of D is given by, 

D = 2 .77-8 . OxlO"M r +l .8xl0" 5 i4r (12) 

-0.8/ (l+exp(250 -E) /75) 

For lithium nuclei because of the 'halos', compression is less and hence the Pauli 
blocking effect is less important and a reduced value of D/3 gives better results 
for the reaction cross sections al the intermediate and higher energies. 

There are no adjustable parameters in the model except that for proton 
nucleus collisions this method of calculaling the Coulomb energy underestimates 
its value for the very light closed shell nuclei of alpha and carbon, and these 
should be increased by a factor of 27 and 3.5 respectively for a better fit. 


3.Results / Conclusions: 


Typical results obtained from the model are shown in Figures 1 through 5. 
Agreement with experimental dala is excellent and is better than all other 
empirical models reported earlier. This is particularly important in view of the fact 
that the agreement is excellent Ihroughoul the whole energy range - uplo a few 
A GeV. The mode! has been tested with all Ihe available data for projectiles 
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proton through krypton and targets alpha through bismuth tor the energy range 
from a few A MeV upto a few A GeV and is tound to give excellent resells (or all 
fhe systems throughout the energy range, tn view of the simplicity and accuracy 
of the model it is a welcome improvement tor transport calculations. 

„ will be interesting to see how the model compares with the new 
experimental data as and when these become available. 
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Figure Captions: 


Figures 1 : 


Figures 2 : 


Figures 3: 


Figures 4 : 


Figures 5 : 


Absorption cross sections for proton - beryllium collision as a function 
of proton kinetic energy. Solid line represents present model and the 
experimental dala are from [15]. 

Absorption cross sections for proton - aluminum collision as a function 
of proton kinetic energy. Solid line represents present model and the 
experimental dala are from [15]. 

Absorption cross sections for alpha - carbon collision as a function of 
incident ion kinetic energy. Solid line represents present model and 
the experimental data are from [21 - 22]. j 

Absorption cross sections for carbon - carbon collision as a function 
of incident ion kinetic energy. Solid line represents present model and 
the expeiimental dala are from [1,9,13,16]. 

Absorption cross sections for calcium - calcium collision as a function 
of incident ion kinetic energy. Solid line represents present model and 
Ihe experimental dala aie from [1,9,13,16]. 
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Abstract: 

We present a simple universal parameterization of total reaction cross sections 
for any system of colliding nuclei valid for the entire energy range from a few A 
MeV to a few A GeV. The universal picture presented here treats proton-nucleus 
collision as a special case of nucleus-nucleus collision, where the projectile has 
charge and mass number of one. The parameters are associated with Hie physics 
of the collision system. In general terms Coulomb Interaction modifies cross 
sections at lower energies and the effects of Pauli blocking is Important at higher 
energies, The agreement between the calculated and experimental data is 
better than all earlier published results. 


Introduction: 

The transportation of energetic ions in bulk matter is of direct interest in 
several areas [1] including shielding against ions originating from either space 
radiations or terrestrial accelerators, cosmic ray propagation studies in galactic 
medium or radiobiological effects resulting from the work place or clinical 
exposures. For carcinogenesis, terrestrial radiation therapy, and radiobiological 
research; knowledge of the beam composition and interactions is necessary to 
property evaluate the effects on human and animal tissues. For the proper 
assessment to radiation exposures both reliable transport codes and accurate 

input parameters are needed. 

One such important input is the total reaction cross section, defined as the 
total minus the elastic cross sections for two colliding ions: 


Or = °T ~ °el 

, n view of its importance tire total reaction cross section has been extensively 
studied both theoretically (1-141 and experimentally 115-24] for the past five 
decades. A detailed list ol reterences is given in Ref.[l. 13. 16], Empirical 
prescriptions have been developed [2-4. ,0. 1 1. 13[ for me total reaction crass 
sections working in various energy ranges and combination of interacting ions. The 
present model works in all energy ranges with uniform accuracy for any 
combinofion of'ln t eracling ions including proton-nucleus collisions and is more 
accurate than earlier reported empirical models [10] which were accurate above 
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100 A MeV but showed large errors upto 25 percent at lower energies. 


Model Description: 


Most of the empirical models approximate total reaction cross section of 
Bradt-Peters form: 

a R = vrSu¥ 3 + A}'* - S ) 2 (2) 

where r 0 is energy independent and 6 is either energy-independent or energy 
dependent parameter, and A p and A T are the projectile and target mass 
numbers, respectively. This form of parameterization work^ nicely for higher 
energies. However, for lower energies Coulomb interaction becomes Important 
and modifies reaction cross sections significantly. In addition, strong absorption 
models suggest energy dependence of the Interaction radius. Incorporating these 
effects, and other effects discussed latter in the text, we propose the following 
form for the reaction cross section: 


a* = 7 irt(At > /3 + A + fi E ) 2 (1 - 


Where r 0 = t.l fm, and Ecm is colliding system cenler of mass energy in 
MeV. The last term is Ihe Coulomb interaction term which modifies the cross 


section at lower energies and gets less important as the energy increases 


(typically after several tens of A MeV). In Eq. (3) B Is the energy dependent 
Coulomb Interaction barrier (right hand factor In Eq. 3). and is given by, 
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( 4 ) 


B = 


1 . 44ZpZ r 

R 


where Z„ (Z,) Is olomlc number of Ihe proiecllle (lo.gef) and radius far evaluating 
the Coulomb barrier height is, 


R = r P + r T + 


1.2 ui /3 + ^ /3 ) 


E 


1/3 


(5) 


cm 


where r, Is equivalent sphere radius and is related to Ihe r fm|j radius by , 

r = 1 29 r i (6) 

' I - 1 - • mts.i 

I 

with ( i = PJ)- 

There Is energy dependence In Ihe reaclion cross seclion al inlermediale 
and higher energies mainly due lo Iwo effects - transparency and Pauli blocking. 
This is taken into account in fl E , which is given by, 


fi E = 1.85S + 0.16S/E C ^ 3 " Ci 
+0.91(i4 r -2Z r )Z P / (i4jv4p) 

where S is. the mass asymmetry term and is given by. 


P) 


1/3 A !/ 3 


s = 


( 8 ) 


1/3 


+ A 


1/3 


and is related to the volume overlap of Ihe collision system. The last lerm on Ihe 
right hand side of Eq. (7) accounts for Ihe Isotope dependence of the reaction 
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cross section. The term C E is related to the transparency and Pauli blocking and 


is given by, 


C E = D (1-exp ( -E/40) ) -0 . 29 2exp ( -E/192) 
x cos ( 0 . 229E 0,453 ) 

Where the collision kinetic energy E is in units of A MeV. Here D Is related to the 
density dependence of the colliding system scaled with respect to the density of 
,2 C +' 2 C colliding system, ie: 

P/4 + P A 

D = i . 7 5— (10) 

P-4 C + P-*c ( 

The density of a nucleus is calculated in the hard sphere model. There is 
interesting physics associated with constant D. This in effect simulates the 
modifications of the reaction cross sections due to Pauli blocking. This effect is 
new and has not been taken into account in other empirical calculations. This 
helps present a universal picture of Ihe reaction cross sections. 

At lower energies (below several tens of A MeV) where the overlap of 
interacting nuclei Is small (and where Coulomb Interaction modifies Ihe reaction 
cross sections significantly) Ihe modifications of Ihe cross sections due to Pauli 
blocking are small and gradually play an increasing role as the energy increases, 
since this leads to higher densities where Pauli blocking gels increasingly 
Important. Interestingly enough for the proton-nucleus case, since there Is not 
much compression effect, a single constant value of D=2.05 gives very good 
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results for all proton-nucleus collisions. For alpha - nucleus collisions, where there 
is a little compression, the best value of D is given by, 

D = 2 .77 -8 . OxlO"M r +l . 8xlO“ 5 Ar (n) 

-0.8/ (1+exp (250-E) /75) 

For lithium nuclei because of the 'halos’ [21, compression Is less and hence the 
Pauli blocking effect is less important and a reduced value of D/3 gives better 
results for the reaction cross sections at the intermediate and higher energies. 

There are no adjustable parameters in the model except that for proton 
nucleus collisions this method of calculating the Coulomb interaction barrier 
underestimates its value for the very light closed shell nuclei of alpha and carbon, 
which are very tightly bound and hence compact. Consequently, for these two 
cases Coulomb barrier should be increased by a factor of 27 and 3.5 respectively 
for a better fit. 

Results/ Conclusions: 

Figures (1-45) show the plot of available results for proton-nucleus, alpha- 
nucleus, and nucleus-nucleus collisions. Figures (6,18) also show comparison with 
Ref. [10]. The data set used for figures (1-5) has been collected from Refs. [15, 23] 
and for figures (6-14) has been obtained from Refs. [16, 17, 22, 23]. There is 
extensive data sel available for C + C system (Fig. 18) and has been taken from 
Refs. [16, 1 7, 23, 24], For Ihe remaining figures data has been collected from the 


6 



compilation ol data sel from Reis. [9. 16-20], The ag.eemeni with experiment is 
excellent and is belief than all olher empirical models reported earlier. This is 
parliculaily imporlanl in view of the loci lhal Ihe agreement is excellenl 
lliroughoul Ihe whole energy range - up to a lew A GeV. We no, ice again lhal 
at lower energy end Ihe cross sections are modilied by Ihe Coulomb inleraclion 
and at Ihe inlermediale and high energy end Pauli blocking effects become 
increasingly imporlanl. II will be inleresling lo see how Ihe model compares will, 
the new experimental data as and when these become available. 


) 
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Figure Captions: 


Figures (1-5) : 


Figures (6-14) : 


Figures ( 1 5-45) : 


Reaction cross sections as a function of energy for proton- 
nucleus collisions. Solid line represents present model. 

Reaction cross sections as a function of energy for alpha- 
nucleus collisions. Solid line represents present model. The 
dashed line in Fig. 6 is from Ref. [10). 

Reaction cross sections as a function of energy for nucleus- 

nucleus collisions. Solid line represents present model. The 

I 

dashed line in Fig. 18 is from Ref. [18]. 
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Abstract 


A fully self-consistent, quantum statistical model is used to 
investigate the emission of light fragments from equilibrated nuclear 
systems with varying prolon-to-neutron ratios over a range of internal 
excitation energies. A very strong correlation between the species of 
emitted fragments and the composition of the equilibrated system is 
found. This finding suggests that light fragments emitted by excited 
nuclear systems may serve as an indicator of the equilibrated population. 
We also note that a portion of the emitted alpha particle spectrum is 


condensed. 
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Obtaining information about nuclear matter under unusual 
conditions of density and temperature (excitation energy) has been a goal 
of research in nuclear physics for some time 11-12). Because nuclear 
collisions may last for only a short period of time (~ 10' 22 sec), the 
formation of unusual matter is often transitory. Because most 
observations are endpoint in nature, information about transitory 
behavior is often lost. In this letter, we suggest that studies of the 
distributions of light particles emitted from an equilibrated nuclear 
systems may be useful for ascertaining the composition of the emitting 
excited system. Although there is some debate concerning the existence of 
equilibrated systems in heavy ion collisions, there is general agreement, 
for the low energies considered herein, that at least a portion of the 
colliding matter does become equilibrated. 

In the work described herein, we show that an unequivocal 
correlation exists between the abundances of the emitted light fragments 
and the overall composition of the equilibrated system. This correlation is 
strong enough that the light fragment abundances may be used as 
indicators of the transitory composition of the nuclear collision 
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participants. We also note that a measurable fraction of the emitted alpha 
particles are condensed. 

To begin, consider an assembly of A nucleons consisting of N 0 
neutrons and Z Q protons in a volume of V Q . Let the proton-to-neulron 
ratio be denoted by y (= Z 0 /N 0 ). Allow the system of A nucleons to evolve 
in thermodynamic equilibrium at temperature T and density p, and to 
emit light particles: neutrons (n), protons (p), deulerons (d), tritons (t), 

3 I Ie (h), and alphas (a). Besides their sets of intrinsic properties (spin, 
mass, etc.) which differ for each particle, we also classify them in terms of y 
so that they can be related to the total equilibrium population A. Specific 
values for each are (yin parentheses): n (0), t (0.5), d (1), a (1), and 3 He (2). 
For protons, y is infinite since we are dividing by zero. In the model, 
charge and baryon number are conserved by requiring that 

Z n.N. = N 

Yj z.N. = 7. (1) 

t r o ' 

where the N ( refer to the number of particles of species i. The system 
evolves in chemical equilibrium, which implies that for any species i, the 
chemical potential is 

it . = z.u + m .it +e. 

r i v p in x 


(2) 
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where z i and n. are the number of protons and neutrons in the i— species, 
and Ej is its binding energy. The distribution of fermions is given by [13] 


N. = 2v.V rt~ 1/2 x : 3 

1 d I 0 I 



(3) 


where g ; = 2 Sj + 1 is the spin degeneracy, p. is the chemical potential X s is 
the thermal wavelength for the E- particle of mass mj 


X. = 2 n ti |2n m. k R 7 j (4) 

where k n is Boltzmann's constant, and the species index is i (= p, n, h, t) as 
appropriate. Values for the Fermi-Dirac integral functions are tabulated in 
the literature [14]. For bosons [13, 15, 16], we have 

r r-1 


N.= 


exp -u./Jt 



-1 


+ g.V XT 3 

0 l 



(5) 


where the first term on the right side of Eq. (5) gives the number of 
condensed bosons, the second term is the number of noncondensed 
bosons, and the species index is i (= d, ex). The Bose-Einstein integral 
functions appearing in Eq. (5) haven been extensively studied elsewhere 
[13,15]. 

The evolution of the system is studied as follows: for a given 
temperature T, density p (= A/V n ), and proton /neutron ratio y, Eqs. (1) - (5) 
are solved self-consistently until a convergent solution is found. The 
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density is then incrementally increased, for the same T and y, until 
another self-consistent solution is found. The set of calculations is 
stopped at any density for which no self-consistent solution can be found. 
The calculations for increasing density are repeated for different 
temperatures, and then the entire process is again carried out for different 
y values. 

To illustrate, consider systems with A = 80 and y = 1.0, 0.5, and 0.25, 
respectively. These choices for the proton /neutron ratio span most 
systems of interest in nuclear and astrophysics studies. The total number 
of equilibrated nucleons (A = 80) is reasonable for collisions involving 
heavier nuclear systems where some but possibly not all of the nucleons 
may be in equilibrium. Temperature of 1, 2.5, 5, and 7.5 MeV are assumed. 
At higher temperatures, the alphas and other composite species will 
undergo significant breakup. 

The results of the model calculations are displayed in figures 1 to 3. 
For y = 1 (fig. 1) the most abundant light fragment species is the alpha 
particle which also has a proton/neutron ratio of y = 1. Note that the 
distributions of alpha particles are very similar for all excitation energies 
(temperatures) considered. For lower densities and high temperatures 
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(5 and 7.5 MeV), deuterons (also with y = 1) appear. The abundances of 
other species with y * 1 are rare. 

Figure 2 displays model predictions for a proton-to-neutron ratio of 
0.5. Here, the most abundant fragment produced is the triton ( 3 II), which 
also has y = 0.5. The next most abundant fragment species is the neutron 
(y = 0), followed by the alpha particle (y = 1) and deuteron (y = 1). Again, as 
was the case with Fig. 1, the most abundant fragment species is that which 
has the same proton-neutron ratio as the total equilibrated population. 

Next, we consider the extreme case for a proton-neutron ratio of 
y = 0.25. Among the emitted species considered herein, none have this 
value of y. Based on the results displayed in figs., 1 and 2, however, we 
would expect that the most likely emitted fragment species would be one 
with a low y value, such as the neutron (y = 0). From fig. 3, we observe 
that this is the case. Neutrons are the most abundant species - followed by 
the triton (y = 0.5). Other species which appear (barely) are those with y = 1. 

Among the light nuclear fragments considered here, the alpha 
particle is the most-tightly bound. Therefore, from binding energy 
considerations, one would expect alpha production to dominate — 
irrespective of the assumed equilibrium conditions. Instead, our 
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calculations indicate that particle emissions from equilibrated systems are 
dictated by the actual composition (y values) of the source regions rather 
than by binding energy considerations. This finding may be useful for 
monitoring the transitory equilibrated population in heavy ion collisions. 

This study also provides indications of another interesting 
phenomena regarding the population of emitted alpha particle — namely 
that a measurable fraction of the emitted alphas are condensed. This 
effect, as expected, is a maximum for the y = 1 equilibrium population 
which favors alpha particle production. In fact, for every temperature and 
density considered when y = 1, most of the emitted alphas are condensed. 
These results are shown in fig. 4, where the number of condensed alpha 
particles is displayed as a function of density, temperature, and y value. 
Note that condensed alpha particles are also present for y * 1. For y = 0.5, 
the trends in alpha particle emission as a function of density and 
temperature are similar to those for y = 1, where an increase in the 
number of condensed alphas with increasing system density is observed at 
each temperature. For y = 0.25, however, we note that very few condensed 
alphas are produced - irrespective of system density and/or temperature. 
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Since we are unable, a priori, to determine the actual composition 
of the population in thermal equilibrium for any nuclear collision, 
experimental detection of condensed alphas might be enhanced by 
selecting the most favorable conditions. From fig. 4, this appears to be 
lower temperatures (T = 1, 2.5 MeV) and lower densities (p < 0.5 normal 
nuclear matter). It is satisfying to note that these conditions are similar to 
those found at the surface of a finite nucleus where, at T = 0 MeV, alpha 
particles are naturally present. Some years ago, an analysis of nuclear 
matter at T = 0 MeV by Mueller and Clark [17] also predicted the presence 
of alpha particles at densities lower than that of normal nuclear matter. 

In this letter, we have demonstrated that the distributions of light 
particles emitted from thermalized, equilibrated nuclear systems are 
mainly dictated by the proton-lo-neutron ratio of the equilibrated source, 
rather than by binding energy. This finding may be useful for monitoring 
the presence of transitory equilibrium in heavy ion collisions. Although 
more realistic calculations should be performed to confirm these results, it 
is anticipated that the general conclusions presented here would not be 


substantially altered. 
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Figure Captions 

Figure 1. Distributions of emitted particles for A = 80 and y= 1 

(N q = Z Q = 40) as a function of temperature and density ratio 

(p o = 0.17 fm' 3 ). Curves for the more abundant fragment species are the 

only ones displayed. 

Figure 2. Distributions of emitted particles for A = 80 and y = 0.5 

(N Q = 53, Z o = 27) as a function of temperature and density ratio 

(p o = 0.17 fm' 3 ). Curves for the more abundant fragment species are the 

only ones displayed. 

Figure 3. Distributions of emitted particles for A = 80 and y = 0.25 
(N q = 64, Z Q = 16) as a function of temperature and density ratio 
(p o = 0.17 fm" 3 ). Curves for the more abundant fragment species are the 
only ones displayed. 

Figure 4. Number of condensed alpha particles as a function of 
temperature and density ratio (p () = 0.17 fm' 3 ) for different compositions 
(y = 1, 0.5, 0.25) of equilibrated matter. 
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Abstract 

A recent set of light ion experiments are analyzed using the GRNTRN transport code 
and the NUCFRG2 fragmentation database generator code. Although the NUCFRG2 
code reasonably represents the fragmentation of heavy ions, the effects of light ion 
fragmentation requires a more detailed nuclear model including shell structure and short 
range correlations appearing as tightly bound clusters in the light ion nucleus. The 
NUCFRG2 code is augmented with a quasielastic alpha knockout model and 
semiempirical adjustments in the fragmentation process allowing reasonable agreement 
with the experiments to he obtained. A final resolution of the appropriate cross sections 
must await the full development of a coupled channel reaction model in which shell 
structure and clustering can be accurately evaluated. 


Introduction 


The need tor accurate transport methods and corresponding aiomic/nuciear database 
tor protection against radiauons in space was demonstrated in a recent review of issues in 
space radiation protection ( 1 ). Although the earliest efforts in code validation experiments 
were placed on the relatively light ion beam of neon because of the potential importance of 
that beam to radiation therapy ( 2 - 5 ), the first beam studied specifically for space radiation 
protection was the iron beam which is the single most important species in long term 
space exposures ( 6 , 7 ). Among the specific issues in the iron beam experiments are the 
production cross sections for light fragments which were in doubt and in fact the reason 
lor development ot the NUCFRG code ( 8 ). Recent experiments at the GSI accelerator in 
preparation for medical therapy has recently provided data on the transport in water of the 
light ions of carbon, nitrogen and oxygen ( 9 - 11 ). Aside from the fact thai the most 
abundant ions in the space environment with charge of three or more is the CNO group 
produced in abundance in star interiors, these are also interesting ions for which the 
NUCFRG code (based on a liquid drop model) is expected to be least appropriate 
( 12 - 14 ). Hence, the current data set is an important test for the fragmentation database 
used iij space studies. 

In the present report, we give a brief review of the current state of the transport 
formalism and application to the analysis of the GSI beam data. The initial results from 
the NUCFRG2 database are found to be in error by 30-40 percent. It is surmised that the 
fault lies in the liquid drop assumption. The nuclear interaction viewed as a combination 
of inelasuc. quasielastic. and nonelastic events leads us to assume the main quasielastic 
event in which specific panicles or clusters are removed from a reasonably stable core 
cannot be accomadated within a liquid drop formalism ( 13 , 14 ) while the highly 
nonelasuc events in which the mass removed from the projectile is fully dissociated is 
well represented by the NUCFRG2 code. One can see these dynamic differences in the 
12 C/ 13 C ratio in the fragmentation of 40 Ar (ratio of 1 ) m comparison with the value for 


fragmentation of l6 0 (ratio nmi 

0f 2 3) as measured b y Tull (15) and Olson et al. (16) 

respectively. The clusier effec, for >«0 projeciiles is la ken herein from the work of 
Cuctnoun and Dubey „3. H, and some addntonal ad hoc ad,us.men.s are made .0 the 

df0P m ° del ^ -ents until a mom complete cluster code ,s 

available. The resultant database is an improvement over the NUCFRG2 code but L, and 

Be fragments remain untested by the GS. expenmen, s. In the absence of such data the 

.urther development of the Custer mode, is pur best hope of resolvmg these cross 
sections. 


Transport Methods 

lh= domain of a spa.ially uniform beam, the Bol.amann transport epuat.on 

assuming a straight-ahead approx, mat, on is useful in evaluating many field puantities 

high-energy ton transport. Even for narrow and directed ion h,., m 

uirecied ion beams, angular and spectral 

corrections have been applied successfully as multiplicative la, 
type./ al x-.wilh energy £ (AMeV) is given as (17) 


factors (3-5). The ion flux of 


p JfA, l0E] 

Sj(E)P:{E) 0 n°' E j) 


T l Ej dE' A i p i( £ ") r 

7 ic Sj{E)Pj{E)h dE '' a jk ( < E ' - E '')9k[x + Rj(E)- Rj{E'),E"] 


u: 


where S,(E } is the stopptng power. P j(E) is rhe nuciear aiienuaiion coefficren, given as 


P j< E ) = '*p[-J 0 i <Tj(E')Aj<IEys j (E')l 


( 2 ) 


Ej -V[*+ *,■<£)] « energy a ,he boundary lor an ion j a , wiih residual 

energy £, and Rj(E) is ihe residual range, £-‘[/? y (£)] = £. We may rewrite ihe 
solution in terms of the Green's function as 


*1' U - £) = X lo G » (*- E - E o MO. E v )dE„ 


C3) 


where G ]m (x, E.E 0 ) satisfies an integral equation similar to Eq. (1) as 

ijlE)Pj(E) f 1 £ Sj(E) P ’j(E) (4) 

Je ■■ dE " a )t < £ ' ■ E ") a km [ * + Ej (E) - Rj (£•).£•' ,E„ j 


where 


C/«i(0,£.£ £I ) = 8j m 5(E-E U ) 


( 5 ) 


We use the Neuman expansion as a penurbative series 


G Jm (a% E,E 0 )-Y cj' J ( a, £. £„ ) 


i=0 


( 6 ) 


where the leading term is 


Arr,.*!*) - x , r r . 
"" ' " j 5, (£)£,(£) ^ 5 ( E J ~ E o ) 


(7) 


and the higher order terms are given by 
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The first iterate of Eq. (8) is given as 


G<;>(*,£.£ 0 ) = 2 \‘<JE 


E, A.Pj(E') 




MEAZdEJd 
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( 9 ) 


Skm d(E' t -E 0 ) 


where 




( 10 ) 


If we make the usual assumption that the interaction is dominated by perpherial processes 
then 


C7 /Hi (£ ',£••) = Oj„ ,(£••) 8(E ' -£") 


(ID 


for which- the second term of the Neuman series becomes 


g|‘J(a.£,£ ( ,) = — ^ 

J,n{ u) Sj(E) 


v, 




V j 


P,(E) 


o„„(n 


W) 


( 12 ) 


tor values of E suet) that 


-^4 «,„(£„)- .v] < R, (E)< *„,(£„) - .v 

v i Vi 


( 13 ) 


where V, is the range scaling parame.er (v, . zj / Aj). h was shown 

by Chun e. al. (It, rhai the second term can be approx, mated by a linear function of 
energy as 


G Jm ( X 'E-E 0 ) ^ f G ; J ( a.\ £, £ 0 nun ) + (x. £. £ 0 ^ )J 

G i>!S X ' E ' E ou\zx )-G\ 1 ,1 (*, £, E 0 nun ) 


£ - F 

o mar tL . o min 


r_ .5r> max + E t> min 


( 14 ) 


when; £„ ^ and E 0 mul are associated whh the allowed range in relation ( 13). Eq. ( 1 4) 
allows a simple numerical eyaluation of the ,on flux spectra 15. IS). The h.gher order 
terms we evaluate using the nonpenurbative method discussed elsewhere (5). 


comparison With Experiment 

Experiments were performed at the GS1 accelerator usmg beams of I2 C . 14 N . and 

' 60 “ e " erSieS ° f 674(±2 » MeV/ “ **» the transmitted flux of charge 5 to 8 were 
measured behind a water large, of vanahle thickness (9-111 The measured transmitted 

,nux of the same charge (open circles) is shown in fig. 1 along with the solution for the 
primary beam flux (dashed curve) and the calculated ilux of all ions of charge equal to the 
imual beam (solid curve). It appears that the total absorption and neutron removal cross 
sections of the NUCFRG2 are reasonably correct. The measured flux with a single 
charge removed (open circles) is shown in Tig. 2. The NUCFRG2 code (filled circles) 
tends to overestimate the single charge removal cross section for I2 C and )<.() projectiles 
and underestimates the cross section for » N projectiles. This unsyslemauc behav , or u 

indicative of structure dependent effects and perhaps results from the fact that the carbon 
and oxygen nuclei consist of integral numbers of highly stable alpha panicles and nitrogen 
does not. Results for the revised NUCRFG cross sections are also shown as the firs, 
collision tom (dashed, md the complete soluuon (sohd) which ts in reasonable agreemen, 


— “ «■-« « <»■ - 

w “ “ *“■ I" m 3 ...» a. «ucnc» 

“ ™ c, " 0! *. »*..,» J „„ «0 „ „ 

— - 

fragments being produced ... ■ 

knockout cross section leaving lhe , 2 C cor£ in ^ ^ 

:: ~r ,so " d ^ bnnss sood ~ ~ - « — -; 

*n e carbon fragments produced on the first rnllicin • » 

in the fi Bure Th P i ■ , n “ Sh ° Wn 25 ^ dashed curve 

L on y triple charge removal data is for the ,6 0 i k 

% 4 for which the NUCFRG" 7 ° ^ “ Sh ° Wn in 

INUCFRG2 is in reasonable agreement. 

As a result of the present comparison the NUCFRm H , 

satisfactory agreement .it,, the tncasuretnents .rHecJ^" ^ ^ 

a 1 , 1 . , , . mjn ei al - ( 17 ) as shown in fie s 

'JT ' nC,US,0n ‘ WS dUSl " **» » ' m Porlam f„ ling lhL . ^ helwcen 

experiment ant, theory. aOhtttona, ad hoc ad,us, meats .ere made the mC L- 

“ be “ r ?■ - OS, data. The ftna, NUCFRGo c “ 

«h those measured a, GSI in table , The use of 2' “ Sl ’° Wn " C ° mP “ S0 " 

solid c , ' curve for the first colliston fl ux and the 

aM ihc h,ehcr ^ ^ 

’ thL Ll Jnd production cross section, , r 

sections arc not represented in these l a .p C , 
measurements and are let, „„„ C Jlesl 

calculations wi„ he heipfu, ^ ^ ^ ^ 

hopefuhy he avai,ab,c the near future' Most im ‘ ^ “ C, ‘ re “ , ‘ S 


targets. There is expected to be a large ^knockout cross section for other 4n nuclei such 
as 20 Ne, 24 Mg, and 28 Si which are important contributors to galactic cosmic ray 
exposures ( 14 ). Also, the knockout of other light clusters will become important heavy 
ion fragmentation for all nuclei which have large spectroscopic constants for clusters 
outside closed subshells in the ground state of the projectile or target. There is strong 
energy dependence from the nuclear form factors and the effects of pion production as 
clearly shown in the few hundred MeV to one. GeV region in fig. 6. Other structure 
dependent effects are expected to show strong energy variations. Fortunately the energy 

dependence is less severe m light targets and low energy which is helpful in developing 
medical therapy beams. 


Concluding Remarks 

The value of having transport experiments to guide semiempirical models used to 

generate nuclear databases for estimation of shielding properties is aptly demonstrated in 

the present paper. The resulting revisions in the NUCFRG2 code will increase its 

uselulness in future studies. It is clear from the present study that a final database 

generator will require cluster models for the light ions. Partial results of such models was 

instrumental in correcting some of the deficiencies in the revised NUCFRG2 code 

presented herein. Clearly future versions should use exclusively cluster models for the 
light ions.' 
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Table 1. Augmented NUCFRG2 charge chtmgtng cross sections (mbn) 
and experimental values in a water target (9) 
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1. INTRODUCTION 

Obtaining information about nuclear matter under unusual conditions of density 
and temperature (or excitation energy) has been a goal of research in nuclear physics 
for nearly two decades [1-12]. In particular, there has been an intrinsic interest in 
Bose- Einstein condensation of particles with a concomitant phase transition in nuclear 
collisions. Studies of multifragmentation reactions may provide information about 
liquid-vapor and other types of nuclear phase transitions, and ultimately information 
about the nuclear matter equation of state. Because nuclear collisions may last for 
only a short period of time (~1(P 22 sec), the formation of unusual nuclear mutter 
may be highly transitory. Since most experimental observations are endpoint, in nature 
(e.g., cross sections or yields), information about transitory behavior is often lost or 
overlooked. In previous work [13] we demonstrated the possibility that Bose-Einstein 
condensation of nuclei is just such a transient phenomenon, which has thus far been 
undetected. Over a decade ago, Bose condensation of pious was predicted to occur 
in the expansion phases of nuclear collisions and other forms of matter [4 6]. To 
date, none have been observed. There is a significant difference, however, between 
Bose-Einstein condensation of nucha and (hat of pious. To form nuclei, matter must 
have some minimum density because several nucleons must be in close physical and 
momentum proximity. On the other hand, pion production is generally a one-nucleon 
phenomenon and therefore does not depend upon a minimum density of nuclear matter. 
Hence, Bose condensation of pious should be observable in end-point measurements but 



condensation of nuclei will not. Its observation must be related to pro-ultimate stages 
of the reaction (e.g., compression). Bose-Einstein condensation has been a well-known 
phenomenon in condensed matter physics and has recently become a vigorous area of 
research in atomic physics [14, 15]. There is no reason for this fundamental phenomenon 
to be absent from unclear physics. 

We begin the paper by briefly reviewing the (jnantum statistical model used to 
investigate the Bose-Einstein condensates [13]. Theoretical results for symmetric 
( Z () = N 0 ) and asymmetric* (Z Q ^ No) nuclear systems are presented (Z () (N a ) refers 
to the number of protons (neutrons) initially present). Then some interesting aspects 
of light fragment emissions from equilibrium nuclear systems are discussed. Finally, 
possible experimental signatures of these Bose condensed nuclei are described. 


2. THEORY 

Since our initial focus was on the possible existence of Bose condensate nuclei, a 
simple model based upon quantum statistics was assumed [13]. It ignores interactions 
between the particles, which should be a reasonable assumption for the low densities 
and temperatures being considered. More detailed and reliable predictions should 
incorporate the effects of particle interactions. In principle, we know how to do this 
[11]; such complexities, however, are left for future work. The present, simple model 
will suffice to demonstrate the possible existence of the phenomenon and to investigate 
qualitatively some of its features. 

2.1 Quantum Statistical Model 

To begin, consider an assembly of A nucleons consisting of N () neutrons and Z () 
protons in a volume V () . The smallness of the finite volumes could mask a clear signal 
of any phase transition. However, by exercising care, experimental evidence of Bose 
condensation may be observable. Assume that, the system of A nucleons is evolving in 
thermodynamic equilibrium at a temperature T and density p to form the light ions: 
12 C, U B, 1() B, 9 Be, 7 U, c Li, 4 He(a), 3 He(h), 3 Il(t), 2 1I(<1), together with protons (p) 
and neutrons ( 71 ). Although there is some debate concerning the existence of fully 
equilibrated systems in heavy ion collisions, there is general agreement that at least 
a portion of the colliding matter does become equilibrated, and can be described by 



equilibrium thermodynamics. In the model, charge and baryoii number are conserved 
by requiring that 


£>i N i = N„ 

i 

^2ziNi = Z 0t (1) 

i 

where the N t refer to the number of particles of species i. 'I he system evolves in 
chemical equilibrium, which implies that the chemical potential is 

Hi = Zillp 1 UiUn -h (2) 

where Zj and 7 /.j are the number of protons and neutrons in the i species, and f j is 
its binding energy. For convenience, the particles are treated as point particles moving 
in a reduced volume given by 

V pt = V 0 -Y^N i V ii (3) 

i 

where V, is the eigenvoluine of the i Ul particle. The number density is 

Po — A I V n , (-*) 


and the density of point particles is 

p v t = A/V pl . (b) 

The distribution of fermions, from quantum statistics, is given by [1C] 

Ni = 2giV v i*~ { l 1 \T i Ffpim/ksT), 

where (j l = 2 S t \ 1 is the spin degeneracy and A; is the thermal wavelength for the i th 
particle of mass m t 

A i = 2nh{2nmik D rr 1 ^ (7) 

where k B is Boltzmann’s constant and the species index is i = (p, u, h, t, 7 Li, 9 Be, n B) 
as appropriate. Values for the Fermi- Dirac integrals 

f°° 1/2 

— (lx 

1 + C x ~ 1 ' 

are well documented in the literature [17]. For bosons quantum statistics yields [1G] 


Ffd( u ) - f 

J () 


( 8 ) 


Ni = </j[exp(-/i/A:0T) ~ J ] 1 + OiVpt^i *F BE (-Pi/k B n 


(9) 



where the first term on the right side of Eq. (9) gives the number of condensed bosons, 
the second term is the number of noncondensed bosons, the species index is z(= <1, «,®Li, 
1 () B, 12 C), and the Bose-Einstein integrals [18] are 



The evolution of the system is studied as follows: for a given temperature 7’, 
density and proton-to-nentron ratio 7 ( = Z„/N 0 ), Eqs. (G) and (9) are solved 
self-consistently, while simultaneously satisfying E<is. ( 1 ) and (2), until a convergent 
solution is found. The density is then incrementally increased, for the same T and 7 , 
until another self-consistent solution is found. The set of calculations is stopped at 
any density for which no self-consistent solution can be found. This occurs when the 
magnitude of the chemical potential for one of the bosonic species is at a minimum. 
The calculations for increasing density are repeated for different temperatures, and 
then the entire procedure is again carried out for another value of 7 . 


3. RESULTS 

We have studied systems with 7 = 1 . 0 , 0.5, and 0.25, and with A = 80, 100,320,610, 
and 1280. These choices for 7 span most systems of interest in nuclear and astrophysics 
studies. The smallest number of equilibrated nucleons (A = 80) is reasonable for 
collisions involving heavier nuclear systems where some but possibly not all of the 
nucleons may be in equilibrium. The largest value of A{A — 1280) is a computationally 
reasonable ensemble of nucleons for representing stellar interiors or infinite nuclear 
matter. Temperatures of 2.5 and 5 MeV were assumed. At higher temperatures, the 
excitation energy would be near or above the average binding energy per nucleon and 
significant breakup of the composite species being formed could occur. 

3.1 Results for Symmetric ( 7 = 1 ) Systems 

Representative results of model calculations for symmetric systems with A — 160 
and 1280 at T = 2.5 MeV are displayed in figures 1 to 3. Results for A = 80 at 
T = 2.5 MeV and 5 MeV were presented elsewhere [13). Figure 1 displays particle 
distributions as a function of the density ratio Ppt/pmn where p um — 0.17 fm ~ 3 is the 




Fig. 1. Distribution of most abundant particles for symmetric (y= 1) systems at T = 2.5 MeV. 
Other species have abundances which are smaller than the lowest curves shown here. 



Fig. 2. Percent of condensed bosons for symmetric (y = 1) systems at T = 2.5 MeV. 
Other bosons have values smaller than those displayed. 



Fig. 3. Chemical potentials for the most abundant bosons for symmetric (y = 1) systems 
at T = 2.5 MeV. Other bosons have values smaller than those displayed. 



usual nuclear matter value. Note that bosons are more abundant than fermions. The 
most abundant ion species are 12 C and 4 He, both with 7=1. Note also that the 
distributions scale with the total mass number A such that for each density ratio, the 
number of particles for a particular species when A = 1280 is eight times larger than 
when A = 160. J his same scaling is present, for the other A values and for T = 5 MeV. 
Figure 2 displays the percentage of condensed bosons. From these results, it is apparent 
that t here is a finite probability of producing detect able numbers of condensed alpha 
pai tides and carbon nuclei at densities which typicall}' occur in heavy ion collisions 
(p £ 2 Pmn)- Note that the percentage of condensed bosons increases with increasing 
density, but the percentages decrease with increasing A. Hence there is no simple 
scaling with mass number. Finally, figure 3 displays the chemical potentials for the 
most abundant bosons ( 12 C and ^He). All other boson species have values which are 
smaller (more negative) than those shown in the figure. Here we note that the chemical 
potentials are nearly identical for both A values. This result holds for all calculated A 
values and is not unexpected since Eq. (2) has no A dependence. Similar features are 
observed for the T — 5 MeV calculations. 

3.2 Results for Asymmetric (77^ 1) Systems 

Results of model calculations for asymmetric systems with 7 = 0.5, A = 160 and 
1280, ami 1 = 2.5 MeV are displayed in figures 4 to 6. From figure 4, we note that 
fermions are more abundant than bosons, which is the opposite trend to that which 
occurs for symmetric systems. However, as was the case for symmetric systems, the; 
particle distributions do scale with total mass number A. Naively, one would expect 
that the most abundant species should be those which are most tightly bound ( 12 C or 
lie). Instead, these results indicate that the most frequent species that are produced 
are those which have 7 values close to the value of the overall equilibruted system. In 
this case, the most abundant species are ^11 (7 = 0.5) and neutrons (7 = 0). The total 
number of bosons produced is (juile small. 

Figure 5 displays the percentages of condensed bosons. These percentages are 
comparable to those displayed in figure 2 for symmetric systems. However, because of 
the small number of bosons produced, the number of condensed bosons is probably too 
small to delect in asymmetric systems. 

Figure 6 displays chemical potentials for the bosons in the 7 = 0.5 system. Note 
again that there is little or no A dependence. Note also that the magnitudes for this 
system (7 = 0.5) are much larger than for the symmetric case (7=1) displayed in 
figure 3. 




Fig. 4. Distribution of most abundant particles for asymmetric (y = 0.5) systems at T = 2.5 MeV. 
Other species have abundances which are smaller than the lowest curves shown here. 



Fig. 5. Percent of condensed bosons for asymmetric (y = 0.5) systems at T = 2.5 MeV. 
Other bosons have values smaller than those displayed. 



Fig. 6. Chemical potentials for the most abundant bosons for asymmetric (y = 0.5) systems 
at T = 2.5 MeV. Other bosons have values smaller than those displayed. 



3.3 Summary of Results 


The two main findings described in the previous sections were: (1) Bose-Einstcin 
condensation of nuclei, such as 4 Ile or 12 C, is a real possibility in equilibrated nuclear 
systems composed of equal or nearly equal numbers of protons and neutrons; and (2) the 
distributions of light particles emitted from thermalized, equilibrated nuclear systems 
are mainly dictated by the proton- to-neutron ratio of the equilibrated source, rather 
than by binding energy. Therefore, the composition of the light fragments emitted in 
the collision may be an indicator of the composition of the equilibrium population of 
nucleons in the source. 


4. EXPERIMENTAL DETECTION OF BOSE-EINSTEIN 

CONDENSATION 

Because we are unable, a priori, to determine the actual composition of the 
population in thermal equilibrium for any nuclear collision, experimental detection 
of condensed nuclei might be enhanced by selecting the most favorable production 
conditions. From the previous section, this appears to be symmetric collision systems 
( Z 0 — N (} = A/2), low temperatures (T — 2.5 MeV or lower) and low densities 
{p < 1.0 pinn). It is satisfying to note that these conditions are similar to those found 
near the surface of a finite nucleus where, at T = 0 MeV, alpha particles are naturally 
present. Some years ago, an analysis of nuclear matter at T = 0 MeV by Mueller and 
Clark [19] also predicted the presence of alpha particles at densities lower than that of 
normal nuclear matter. 

The critical question is how to observe; Bust; condensate nuclei once they are 
produced. For any distribution of reaction products, total yields or cross sections 
are most commonly measured. A measurement of this type, however, would neither 
confirm nor deny the existence of condensed nuclei. Figure 1 indicates that the 
total yields of bosons (sums of condensed and noncondensed) vary with temperature 
and density. Figure; 2 indicates that the percentages of condensed bosons also vary. 
Therefore, experimental observations must be able to distinguish between condensed 
and noncondensed particles. The most promising method appears to be measurements 
of the momentum distributions of each boson species. In the collision eenter-of- 
inass system, condensed bosons will be produced with zero/near-zero momenta. The 
noncondensed bosons will have much larger momentum values in that coordinate 
system. Therefore, it appears that the clearest, most unambiguous signal that 



Bose-Einstein condensation has occurred is to detect bosons with near-zero momenta in 
the center-of-mass system. In the laboratory, the condensed bosons would be detected 
at 0° to the beam with a momentum equal to that of the collision system center of 

mass. 


5. CONCLUDING REMARKS 

In this paper we have used a fully self-consistent quantum statistical model to 
investigate the phenomenon of Bose- Einstein condensation of nuclei in symmetric and 
asymmetric ensembles of nucleons in thermal and chemical equilibrium. It is found that 
Bose condensation of nuclei is favored in symmetric or nearly symmetric, equilibrated 
nuclear systems, at low temperatures (T = 2.5 MeV or below) and low densities 
(densities at or below normal nuclear matter). The most direct evidence of their 
existence is to experimentally detect bosons with little or no momenta in the system’s 
center-of-mass frame. 

We also found that the distributions of light particles emitted from thermalized, 
equilibrated nuclear systems are mainly dictated by the proton- to-neutron ratio of the 
equilibrated source, rather than by binding energy. This finding may be useful for 
monitoring the presence of transitory equilibrium in heavy ion collisions. 

Finally, we emphasize that the findings presented in this work, although definitive, 
deserve further study using more sophisticated models. Nevertheless, the specific 
experimental signature of Bose- Einstein condensation of nuclei proposed in our work 
will not be substantially altered by the use of improved theoretical models. In view of 
the importance of the phenomenon to fundamental physics, experimental verification 
of its existence or nonexistence should be a high priority. 
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ABSTRACT 

The fragmenting of high energy, heavy ions (HZE particles) by hydrogen targets is an important, 
physical process in several areas of space radiation research. In this work quantum mechanical 
optical model methods for estimating cross sections for HZE particle fragmentation by hydrogen 
targets are presented. The cross sections are calculated using a modified abrasion-ablation collision 
formalism adapted from a nucleus-nucleus collision model. Elemental and isotopic production cross 
sections are estimated and compared with reported measurements for the breakup of neon, sulphur, 
and iron, nuclei at incident energies between 400 and 910 MeV/nucleon. Good agreement between 
theory and experiment is obtained. 

INTRODUCTION 

There is a need for reliable methods of accurately estimating cross sections for high-energy heavy 
ion breakup by hydrogen for a variety of space radiation applications. In astrophysics, interstellar 
hydrogen comprises the major type of material encountered by galactic cosmic rays (GCR) as they 
travel through the universe. Hence, accurate cross sections are crucial for understanding cosmic-ray 
propagation and source abundances. Hydrogen is also a major constituent of human tissue and 
appears to be the most effective GCR shield material per unit mass for long duration, manned space 
missions /l, 2/. Therefore, accurate cross sections are needed for proper risk assessment of critical 
organ exposures of astronauts. Typically, cross sections used in many of these studies have been 
obtained from semiempirical parameterizations /3, 4/ which have various fitting parameters. In this 
paper we present fundamental, quantum-mechanical, optical model methods for estimating these cross 
sections using a knockout (abrasion) - ablation collision formalism obtained from nuclear scattering 
theory /5, 6/. The model has no arbitrary fitting parameters. Predictions of element and isotope 
production cross sections are in good agreement with recently reported measurements /7/. 
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THEORY 

In an abrasion-ablation model, the projectile nuclei, moving at relativistic speeds, collide with 
stationary target nuclei. In the abrasion step (particle knockout), those portions of the nuclear 
volumes that overlap are sheared away by the collision. The remaining projectile piece, called a 
prefragment, continues its trajectory with essentially its precollision velocity. Because of the 
dynamics of the abrasion process, the prefragment is highly excited and subsequently decays by the 
emission of gamma radiation or nuclear particles. This step is the ablation stage. The resultant 
isotope is the nuclear fragment whose cross section is measured. 

Although a picture of overlapping nuclear volumes being sheared off may be reasonable for heavier 
nuclei colliding with each other, it is not reasonable for a single nucleon striking another nucleus. 
Instead, a more reasonable physical picture involves individual collisions between the projectile 
constituents and the target proton. Some struck projectile nucleons exit the fragmenting nucleus 
without further interaction, and some interact one or more times with the remaining constituents 
before departing. These interactions are called frictional-spectator interactions (FSI). The remaining 
nucleus (prefragment), in an excited state because of the energy deposited during the collision, then 
de-excites by particle- or gamma-emission processes. This picture is easily described by an abrasion- 
ablation-FSI model where the abrasion stage is described by a quantum-mechanical, optical model 
knockout formalism, and the ablation stage is modeled with cascade-evaporation techniques. There is 
no excess surface area energy. Instead, the prefragment excitation energy is assumed to be provided 
by FSI contributions from the abraded nucleons. 

In the optical potential knockout (abrasion) formalism /6 /, the cross section for producing a 
prefragment of charge Z PF and mass A PF , because of collision with a hydrogen target, is given by 
equations (1) - (3) of ref. /5/ with/^ = 1 and P r (^ T ) = 6 (£ r ) . 

Prefragment excitation energies are estimated from the FSI energy contribution which is calculated 
with a modified form of the model of Rasmussen /8/. In the modified model, the average energy 
deposited in the prefragment for each FSI is 

< E psi > m 8.5R MeV (1) 


where the nuclear radius is R = 1.29 R 

rim' 

Therefore, the cross section for a prefragment species (Z PF , A PF ) which has undergone frictional 
spectator interactions is obtained from eq. (5) of ref. /5/ and, the final hadronic cross section for 
production of the type i isotope is obtained from eq. (6) of ref. /5/. Finally, the elemental production 
cross sections are obtained by summing all isotopes of a given element according to 

a (z.) = £ a (z,a) 

"“ c l v A- nucy i' ij 


(2) 
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Fig. 1. Isotope production cross section for -^S at 400 MeV/nucleon fragmenting in 
hydrogen targets. For clarity only even - Z (charge number) isotopes are displayed. 
Also displayed are recently reported experimental measurements 111 . 
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Fig. 2. Element production cross section for ^Ne at 910 MeV/nucleon and for ^Fe at 
400 MeV/nucleon fragmenting in hydrogen targets. Also displayed are recently reported 
experimental measurements /7/. 
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RESULTS 

Figure 1 shows isotope production cross sections for 32 S beams at 400 MeV/nucleon fragmenting in 
hydrogen targets. The figure also shows recent experimental measurements /7/. For brevity, only the 
even - Z (charge number) fragment isotopes are shown. Element production cross sections are 
displayed in Figure 2 for 22 Ne beams at 910 MeV/nucleon, and for 56 Fe beams at 400 MeV/nucleon. 
The figure also shows recent experimental measurements /7/. For these collisions, the agreement 

between theory and experiment is good, especially considering that no arbitrary fitting parameters are 
in the theory. 

CONCLUSION 

A simple optical potential knockout-ablation fragmentation model has been presented for use in 
studies of high-energy heavy ion breakup by proton. The model has no arbitrarily adjusted 
parameters. Model predictions are in good agreement with reported laboratory measurements of 

fragment production cross sections for the fragmenting of sulfur, neon, and iron beams by hydrogen 
targets. 
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ABSTRACT 

The fragmenting of high energy, heavy ions by hydrogen targets is an important physical process in 
several areas of space radiation protection research. Quantum mechanical, optical model methods for 
calculating cross sections for particle fragmentation by hydrogen have been developed from a modified 
abrasion-ablation collision formalism. The abrasion stage is treated as a knockout process which leaves 
the residual prefragment in an excited state. In the ablation stage the prefragment deexcites to produce the 
final fragment. The prefragment excitation energies are estimated from a combination of liquid drop 
model considerations and frictional-spectator interaction processes. Estimates of elemental and isotopic 
production cross sections are in good agreement with published cross section measurements. 

INTRODUCTION 

Reliable methods of accurately and quickly estimating cross sections for the breakup of high-energy 
heavy ions (HZE particles) by hydrogen targets are needed for a variety of space radiation protection 
applications including risk assessment and shielding of astronauts from galactic cosmic rays (GCR). 
Many hydrogen target cross sections used in these studies are obtained from the semi empirical 
parameterizations of Silberberg et al (1976) and Webber et al (1990), which have various fitting 
parameters. In previous work (Townsend 1994, 1996) we have taken a more fundamental approach using 
a knockout-ablation collision formalism based upon quantum-mechanical, optical model methods 
obtained from nuclear scattering theory. The model had no arbitrary fitting parameters. In this work we 
improve the model by incorporating previously-neglected excitation energy contributions from liquid drop 
model considerations. The improved model is used to estimate element and isotope production cross 
section for a variety of projectiles ions colliding with hydrogen targets. 

THEORY 

In the knockout-ablation fragmentation model, the projectile nuclei, moving at relativistic speeds, collide 
with stationary target protons. In the knockout step individual collisions between the projectile 
constituents and target proton occur. Some of the struck projectile nucleons exit the fragmenting nucleus 




31 9991 d 1 


•99 t>08 esn ana yyaNbdS ddsu dd si: so 96, ?t nnr 



without further interaction, and some interact one or more times with the remaining constituents before 
departing. These interactions are called frictional-spectator interactions (FSI). The remaining nucleus, 
called a prefragment or spectator is in an excited state because of the energy deposited during the 
collision. The prefragment subsequently decays by particle - or gamma - emission processes to produce 
the final fragment species whose cross section is measured. 

In the optical potential knockout formalism (Townsend 1994) the cross section for producing a 
prefragment of charge Z PF and mass A PP , because of collision with a hydrogen target, is given by 


K»D-nwrp»)r <» 

where 

T{ b)«exp{-o(e) [2* B(e)J 3/J J dl a +*<. + y) ex p|V t2B(e)] } (2) 


In Eqs. (1) and (2), b is the impact parameter vector, e is the two-nucleon kinetic energy in their center-of- 
mass frame, z„ is the target proton position in the projectile nucleus rest frame, y is the projectile- 
nucleon-target proton relative separation vector, p, is the projectile nucleus ground state number density 
distribution, a ( e ) is the constituent - averaged nucleon-nucleon cross section, and B(e) is diffractive 
nucleon-nucleon scattering slope parameter. Equation (1) assumes a hypergeometrical charge dispersion 
model to describe the distribution of knocked-out nucleons z out of Z projectile protons and n out of N 
projectile neutrons are knocked out, Note that N + Z= A p and A PF = A p - n - z where A p is the projectile 
nuclear mass number. 

Prefragment excitation energies are computed from E„ e - E IP + E FSl . E FSI denotes the FSI energy 
contribution which is calculated using a modified form of the model of Oliveira et al (1979). In the 
modified model, the average energy deposited for each FSI is <E fS i > * 5R MeV where the nuclear radius 
is R - 1 .29 Rn™. The liquid drop model energy E LD is estimated by considering the differences in volume 
energy, surface energy and coulomb energy between the pretragment and the incident projectile nucleus. 
For * 6 Fe it is approximately 2 MeV/knockout. For 32 $ it is about 4 MeV/knockout. 

Depending upon the magnitude of its excitation energy, the prefragment decays by particle and/or photon 
emission. The probability a y (q) that a prefragment species j, which has undergone q FSI, decays to 
produce a final species i is obtained using the EVA Monte Carlo cascade-evaporation computer code 
(Morrissey 1979). Therefore, the final cross section for producing fragment species i is 

o„ k (Z,,4)=X ”'L CL uW a ‘»'(?J' A i’ q ) (3) 

J f-0 

where the summation over j accounts for contributions from different prefragment isotopes j, and the 
summation over q accounts for the effects of different FSI excitation energies. Finally, the elemental 
production cross sections are obtained by summing over all isotopes of a given element according to 
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RESULTS AND CONCLUSION 


Table 1 lists element production cross sections for i6 Fe beams at 40QA MeV fragmenting in hydrogen 
targets. Results for the previous model (Townsend 1996) and the current work are displayed along with 
the experimental measurements of Guzik et al (1994). Note that the current model predictions are in 
much better agreement with the experimental data. Table 2 displays predicted isotope production cross 
sections for 32 S beams at 600A MeV fragmenting in hydrogen targets. The agreement between theory and 
experiment is reasonably good, especially considering the lack of arbitrary fitting parameters in the 
theory. 
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Element Produced 


Mn 

Cr 

V 

Ti 

Sc 

Ca 

K 

At 

Cl 

S 

P 

Si 

Al 


Experiment 


117.47 ± 70.54 
108.54+ 6.21 
72.98 ± 4.25 
65.25 ± 3.46 
40.14 ± 2.92 
29.12 ± 2.40 

17.74 ± 2.12 
13.43 ± 2.00 
7.36 ± 1.99 

5.74 ± 2.41 
3.25 ± 2.76 
3.86 ± 2.08 
7 51 + 2.73 


Cross Section (in millibarns) 
Previous Work This Work 


166 

115 

85 

83 

57 

35 

36 
26 
19 
16 
9 

6.8 

2.5 


195 

120 

84 

82 

59 

33 

28 

18 

13 

9.3 

4.6 

2.5 

1.5 


„ O *• „ c Fnr 600 AMeV 32 S Beams Colliding With Hydrogen. For 
™' 2 ' Experimental Data Are From 

Guzik et al (1994). 


Isotope 

Produced 


29 


Si 
Si 
Si 
Si 
’Mg 
25 Mg 
'Mg 
Mg 
’Ne 


26 

27 

26 


24 

23 

22 


2l Ne 

2 »Ne 

l9 Ne 

i« 0 


17 


16 

15 


O 

o 


Cross Section 
Experiment 


(in millibarns) 

This Work 


O 


18.63 ± 2.14 
40.01 + 3.11 
38.10 ± 3.20 
4.98 ± 1.04 
13.43 ± 2.21 
27.58 ± 3.14 
25.90 i 3.05 
2.76 1 0.96 
7.39 ± 1.95 
14.63 ± 2.86 
11.16 ± 2.07 

I. 27+ U6 
3.37 ± 2.35 
6.78 ± 1.83 

II. 81 ± 2.60 
2.09+ 0.81 


21.7 

55.7 
44.1 
10.6 

9.7 

22.7 
18.5 

4.9 

3.1 

9.7 

5.9 
1.0 

1.1 

3.7 

7.1 
2.6 
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Abstract 


A semiempirical abrasion- ablation model has been successful in gen- 
erating a large nuclear database for the study of high charge and energy 
(HZE) ion beams, radiation physics, and galactic cosmic ray shield- 
ing. The version reported herein has coulomb trajectory corrections, 
improved transmission factors, improved surface energy corrections, and 
light fragment emission was added . The cross sections that are gener- 
ated are compared with measured HZE fragmentation data from various 
experimental groups. A research program for improvement of the data- 
base generator is discussed. 


Introduction 

An adequate and reliable nuclear database that 
assesses the quality of heavy ion beams for various 
technological efforts is needed. For example, the 
nuclear fragmentation properties of shielding mate- 
rials can after the protection of astronauts by an or- 
der of magnitude through the selection of appropri- 
ate shield materials (refs. 1 and 2). The radiation 
quality of heavy ions, which is related to the ability 
to cause biological injury, is an essential parameter 
in high altitude commercial aviation, radiotherapy, 
microelectronic signal processing, and information 
storage. Understanding single event upset damage 
to microelectronic systems is becoming more impor- 
tant as more aircraft and spacecraft control functions 
are handled by microprocessors. This damage is of 
special concern for miniature spacecraft in which re- 
duced telemetry requires intensive onboard process- 
ing by low T power microprocessors with large memory. 
Such small scale devices are very sensitive to single 
event upsets and evaluation of onboard shield worth 
is critical to an adequate design. The specification of 
the nuclear fragmentation cross sections is critical in 
all these applications. 

Over the years the theoretical description of 
nuclear fragmentation in heavy ion collisions has 
been described with abrasion-ablation models (refs. 3 
to 11 and recent work by Cucinotta, Townsend, and 
Wilson of Langley Research Center). In these mod- 
els, fragmentation occurs in two stages. In the fast 
abrasion stage, the projectile and the target over- 
lap and matter is sheared away from both nuclei. 
The remnants of the colliding nuclei are called the 
prefragments (projectile or target) and are assumed 
to be left in a state of excitation. The ablation 
stage is the description of the decay of the pre- 
fragment nuclei. The emphasis of these models is 
typically the prediction of inclusive mass yields of 
the final fragments that are observed. In the semi- 
empirical descriptions of these reactions (refs. 3 to 6), 
the overlap volume of projectile and target is esti- 


mated by using a classical approach. The excitation 
energy of the prefragments is estimated by using a 
surface distortion model with correction terms and 
energy transfer across the interface of the interac- 
tion zone. These models provide reasonable overall 
agreement with measured data; however, they lack 
a description of nuclear structure effects and a de- 
scription of the nuclear diffuseness related to skin 
thickness. A fundamental and more complex prob- 
lem is the degree to which the distribution of levels 
of prefragment nuclei must be considered to provide 
the correct description (ref. 5). 

Hiifner, Schafer, and Schiirmann used the 
Glauber model in a first attempt at formulat- 
ing a quantum mechanical abrasion-ablation model 
(ref. 7). In this model, closure is made on the fi- 
nal states of the target in describing the projectile 
fragmentation and the unobserved nucleons abraded 
from the projectile. Energy conservation is also ig- 
nored and a final closure approximation is assumed 
for the prefragment states that occur following the 
removal of a fixed number of nucleons. The advanta- 
geous factorization properties in the Glauber model 
of the nuclear amplitude then allow closed-form ex- 
pressions for the abrasion cross sections to be found. 
The Glauber model of the abrasion cross section can 
then be shown to correspond closely to the semi- 
empirical models when the abrasion cross section for 
a given product relates to the volume of the projec- 
tile and the target removed in their overlap. A study 
of the closed-form expression for the nuclear absorp- 
tion cross section in the eikonal form of the first-order 
optical potential model led to a recasting of the abra- 
sion model as an optical model bv using the binomial 
distribution (refs. 8, 9, 12, and 13). A comparison of 
the abrasion cross sections (ref. 7) with the optical 
models (refs. 8, 9, and 12) reveals that the two differ 
only by the assumptions of coherence and closure in 
the projectile intermediate states. The optical mod- 
els being preferred for the sum rule on the abrasion 
cross sections to satisfy unitarity. 


The Glauber model or the optical model of abra- 
sion began to employ sophisticated evaporation or 
cascade/evaporation codes to describe the ablation 
stage. These codes rely on a correct average exci- 
tation energy to be used to start the evaporation 
process. A major shortcoming in the physical de- 
scription results from the use of closure on the pre- 
fragment final states in the Glauber model (ref. 7). 
In the Glauber model, all information on the distri- 
bution of the actual levels excited from abrasion in 
the prefragment state is lost and is replaced by an 
average state that is described uniquely by the abra- 
sion cross section and average excitation energy. For 
light prefragments (mass number (A) is less than 16) 
where nuclear structure effects are large and reso- 
nance levels separated by several MeV, the use of an 
average prefragment state is highly questionable. 

The use of a cascade model in the ablation stage 
is also noted. In the abrasion stage, nucleons in 
the projectile are knocked into the continuum by 
the target. These escaping nucleons will multiple 
scatter inside the projectile and cause further nucleon 
knockout. This process may be described as final 
state interaction (FSI) with the Moller operator and 
is expected to be highly dependent on the trajectory 
of the initial cascading nucleon (ref. 14). In the 
optical model (ref. 12) and the semiempirical model 
(refs. 5 and 6), a final state interaction correction 
to the prefragment excitation energy has been used 
to mimic the cascade effect of projectile knockouts. 
It is unclear whether the cascade effects described 
by the Monte-Carlo codes for describing ablation 
are distinct from the FSI corrections that describe 
prefragment nuclei. It was further found in the 
semiempirical model that statistical fluctuations in 
the FSI correction must be considered (ref. 5). 

In a more recent formulation of the abrasion- 
ablation model that uses the Glauber amplitude, 
three major improvements are made (refs. 10 and 11 
and recent work by Cucinotta, Townsend, and 
Wilson of Langley Research Center). First, energy 
conservation is treated in describing nuclear abra- 
sion. Second, the treatment of the excitation of spe- 
cific levels in the prefragments is considered for the 
first time. Here the prefragment excitation is consid- 
ered as a core excitation during the knockout stage of 
abrasion. By using fractional parentage coefficients 
to couple nucleons to the core (prefragment) in the 
projectile ground state, we expect that for many nu- 
clei the complex configurations of the nuclear ground 
state that result in virtual states of relative excita- 
tion are such that the core excitation is dominated by 
diagonal transitions to excited states of the prefrag- 
ment. As the number, of nucleons lost in abrasion 


becomes large, the use of parentage coefficients to 
form the prefragment level spectrum will become in- 
tractable. For many nucleon knockouts, off-diagonal 
coupling is expected to become more dominant and 
statistical methods will become necessary to deter- 
mine the distribution of prefragment levels. The 
preequilibrium models developed by Feshbach et al. 
(ref. 15), Griffin (ref. 16), and Tamura, Udagawa, and 
Lenske (ref. 17) may be amenable to the description 
of the heating of prefragment nuclei in heavy ion col- 
lisions. The third development in the reformulation 
of the abrasion- ablation model is the description of 
cluster knockout in nuclear abrasion. This descrip- 
tion allows the treatment of nuclear structure effects 
in nuclear abrasion, which is important for many 
projectile and target nuclei of interest (C, O, and 
Ne). The reformulation also considers the momen- 
tum distribution for nucleon production from abra- 
sion (ref. 18) as well as ablation (ref. 19). Although 
this more systematic approach will ultimately meet 
the need for a high quality nuclear database, current 
interaction studies rely on the semiempirical model. 

The genesis of the semiempirical model is in the 
abrasion- ablation model of Bowman, Swuatecki, and 
Tsang (ref. 4) as discussed in references 5 and 6. 
In the development of the model, the transmission 
factors of a projectile and target were averaged and 
included to account for the mean free path in nu- 
clear matter (ref. 6). Then, a semiempirical higher 
order correction was given to the surface deforma- 
tion energy of the abrasion products (ref. 6) and the 
energy transfer across the interaction zone bound- 
ary was treated as a two- valued distribution (ref. 5). 
The final charge distribution of the fragmentation 
products was approximated by Rudstam’s formalism 
(ref. 6). The available experimental nuclear fragmen- 
tation data were very limited at the time of the first 
reporting of the model (ref. 6). The model agreed 
with experimental data to the extent that the exper- 
imentalists agreed among themselves. 

Since the inception of the NUCFRG model, ad- 
ditional experimental fragmentation cross sections 
(refs. 20 to 23) and thick target fluence data (refs. 24 
and 25) have become available for validation of 
the database. Incremented improvements have been 
made by including the coulomb dissociation con- 
tribution (ref. 26), an energy-dependent nuclear 
mean free path (ref. 27) based on the analysis of 
Dymarz and Kohmura (ref. 28), and the nuclear 
radii extracted from experimental charge distribu- 
tions (ref. 27). These improvements were included 
in the publicly released version of the HZEFRG1 
code (ref. 27). In addition to these improve- 
ments, the version of the code described herein 
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also includes the following. A coulomb trajectory 
correction that includes the effects of an energy 
downshift has been added (ref. 29). The trans- 
mission factor is evaluated at the maximum over- 
lap in the interaction zone. The spectator nucle- 
ons in the interaction zone are now assumed to be 
poorly bound to the spectators of the abraded frag- 
ment outside the interaction zone and undergo pre- 
equilibrium emission. A unitarity correction is made 
for targets with A > 63. Finally, a correction to 
the semiempirical excess surface distortion energy is 
made for light projectiles. 


by coulomb excitation, which is discussed by Norbury 
et al. (ref. 30). At smaller distances, the overlap of 
the nuclear densities strongly interact and mass is 
removed from the projectile and the target. 

Abrasion Process Description 

The strength of the interaction varies over the 
interaction zone. The projected interaction potential 
on the impact plane is given as 

J V{R+z)dz^ap p C P {R) p T C r (R) ( 6 ) 


In the report presented herein, the model data- 
base is compared with available experimental data. 
Weaknesses within the model are thereby uncovered 
and a research program for database improvement is 
outlined. 


Semiempirical Fragmentation Model 

The equation of motion for a low energy ion in a 
nuclear coulomb field is given by energy conservation 
as 

1 .9 t 2 ZpZpe 2 


*** ~ + 2 ^ 


+ 


(l) 


where Ftot is the total energy in the center of mass, 
r is the relative distance between the charge centers 
with time derivative r,/i is the reduced mass, £ is 
the angular momentum, Zp and Zp are the atomic 
numbers of the projectile nucleus and target nucleus, 
respectively, and e is the electric charge. (That is, 
e 2 = 2 Rya 0 , where Ry is the Rydberg constant and 
a 0 is the Bohr radius.) The angular momentum is 
given as 

( 2 ) 


ft = 2pE tot b 2 


The distance of closest approach is given by equa- 
tion ( 1 ) for r = 0 as 


c _ E tot b 2 t ZpZpe 1 

•^tot — 9 1 


( 3 ) 


which is written as 


b 2 = r(r - Tm) 


( 4 ) 


where R is a position vector in the impact plane, z 
is the longitudinal position component, a is the two- 
body cross section that includes Pauli blocking, p p 
and p T sure the projectile and target mass density, re- 
spectively, and Cp(R) and Cp(R) are the projectile 
and target chord, respectively, at position R along 
z. (See ref. 31.) The amount of nuclear material re- 
moved from the projectile in the collision at a given 
impact separation is the volume of the overlap re- 
gion times an attenuation factor that is evaluated at 
the maximum product of the chords in equation ( 6 ). 
The formula for the number of participating projec- 
tile constituents in the interaction zone is 

A ab r = F Ap[l — exp(— C j7 A)] (7) 

where Ap is projectile mass number, A is mean free 
path, and Cp is the value of Cp(R) that maximizes 
equation ( 6 ) and is given for r T > r p as 

{ 2y]rl — x 2 x > 0 1 

\ ( 8 ) 

2 v r T ~ r 2 x < 0 j 

where 

x = (rj> + r 2 - r|0/(2r) (9) 

and r p and r T are the projectile and target 
radius, respectively, and are related to the root- 
mean-square charge radius (rnns) of electron scat- 
tering ss 1.29^7*2^ — 0.84 2 ^ (ref. 31). When 
r p > r T , the chord Cp is given as 


where 


r 777 — 


ZpZpe 1 

Etot 


( 5 ) 


Note that r m is the distance of closest approach for 
zero impact parameter. 


At a given impact parameter, there is a distance 
of closest approach r for which the interaction takes 
place. When r is large, the interaction is dominated 



x > 0 
x < 0 


in which 

x — (vp -hr 2 — rp)/( 2 r) 


( 10 ) 

( 11 ) 


The quantity F in equation (7) is the fraction of 
the projectile in the interaction zone as given in the 
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appendix (refs. 5, 6, 32, and 33). (The b in ref. 5 is 
replaced by the r in the equations presented herein.) 
The number of projectile spectator constituents in 
the interaction zone is given as 

Aspc = ApF exp(-CV/A) (12) 

The spectator constituents are assumed to be only 
loosely bound to the projectile constituents outside 
the interaction zone. The nuclear mean free path is 
taken as 

A = 16.6/£°' 26 (13) 

where E is the projectile energy in MeV/nucleon. 
(See refs. 27 and 28.) The change ratio of the re- 
moved nucleons is assumed to be that of the initial 
projectile nucleus, an assumption that ignores polar- 
ization effects. 

Surface Distortion and Collisions! 
Excitation 

The projectile constituents outside the interac- 
tion zone (spectators) retain the same relation among 
themselves after the collision as before the collision. 
(This retention of relationship is a sudden approxi- 
mation that is strictly applicable at higher energies.) 
The mass removed by the interaction has altered the 
overall stability of the spectators. This instability is 
related to the reduced binding energy when the nu- 
clear surface is other than its minimum energy spher- 
ical configuration. 

The excess surface area is given as 

AS = 4t wr 2 P [1 + P - ( 1 - F) 2 > z ] ( 14) 

The functions F and P are defined in the appendix. 
(See ref. 32.J For small surface distortions AS in 
units of /m , the excitation energy E§ in units of 
MeV is approximated by 

E$ = 0.95 AS (15) 

At the impact separations r<r ? + r r , the projectile 
spectator group is left far from equilibrium and the 
0.95 MeV/fm 2 coefficient requires correction that is 
taken herein as a semiempirical parameter / given by 

/ = 1 + 5F + [1500 — 320(Ap — 12)] F z (16) 

where the quantity in the square brackets is limited 
to values between 0 and 1500. The cubic coefficient 
F 3 provides a correction for light projectiles that are 
unstable because of large surface distortions. The 
semiempirical surface excitation energy is then 

E s = E' s f (17) 


Note that the correction factor approaches 1 as A a b r 
becomes small. As before, we assume that fragments 
with a mass number of 5 are unbound, 90 percent 
of fragments with a mass number of 8 are unbound, 
and 50 percent of fragments with a mass number of 
9 are unbound. 

A second source of excitation energy is the trans- 
fer of kinetic energy of relative motion across the in- 
tersecting boundaries of the two ions. The rate of 
energy loss of a nucleon when passing through nu- 
clear matter is taken as 13 MeV/fm (ref. 34). The 
energy deposit is assumed to be symmetrically dis- 
persed about the azimuth so that at the interface 
6.5 MeV/fm per nucleon is the average rate of en- 
ergy transferred into excitation energy. This energy 
is transferred in single particle collisions. In half the 
events, the energy is transferred to excitation energy 
of the projectile and in the remaining events, the 
projectile excitation energy remains unchanged. The 
estimate of this contribution is made by using the 
length of the longest chord Ci in the projectile sur- 
face interface. This chord length is the maximum 
distance traveled by any target constituent through 
the projectile interior and is given by 

{ 2 {vp + 2 rrx — r 2 — rf) 1 ! 2 ( r > rp) 1 

r 

2 r P (r < r T ) ) 

(18) 

The number of other target constituents in the par- 
ticipant and spectator interface may be found by esti- 
mating the maximum chord transverse to the projec- 
tile velocity that spans the projectile surface interface 
Ct which is given by 

Ct = 2 (r 2 p - b 2 P ) 1/2 (19) 

where 

bp = (rp -hr 2 - rf)/(2r) (20) 

The total excitation energy transferred across the 
participant and spectator interface is then 

E' x = 13 Q + - 13 C t {C t - 1-5) (21) 

o 

where the second term contributes only if Ct > 
1.5 fm. The effective longitudinal chord length for 
these remaining nucleons is assumed to be one third 
the maximum chord length. 

Nuclear Ablation 

The decay of highly excited nuclear states is dom- 
inated by heavy particle emission. In the present 
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model, a nucleon is assumed to be removed for every 
10 MeV of excitation energy and is determined by 

^abl = (&S E x )/lO A spc (22) 

where A spc are the loosely bound projectile specta- 
tors in the interaction zone that are emitted prior to 
the equilibrium deexcitation process. (See ref. 4.) In 
accordance with the previously discussed direction- 
ality of the energy transfer, the E x is double valued 



' K 


II 

1° 

( p ‘ = 9 


where Pj is the corresponding probability of occur- 
rence for each value of E x in the collision. 

Nuclear Abrasion- Ablation Model 

The number of nucleons removed through the 
abrasion-ablation process is given as a function of 
impact parameter b as 

AA = A a , hl (b) + A M (b) (24) 

The impact parameter is related to the impact sepa- 
ration r by equation (4) for a coulomb trajectory. 

A second correction to the trajectory calculation 
comes from the transfer of kinetic energy into bind- 
ing energy during the release of particles from the 
projectile. (Obviously, energy is also lost in releasing 
particles from the target, which we do not yet cal- 
culate.) The total kinetic energy in passing through 
the reaction zone is reduced to 

E f = E t - 10 A A (25) 

which assumes that 10 MeV is the average binding 
energy*. The kinetic energy* used in the closest ap- 
proach calculation is the average of the initial and 
the final energies and is given as 

Exot = \{E { + E f ) = E t - 1(10 A A) (26) 



Figure 1. Nucleon removal number as a function of impact 
parameter in carbon-copper collisions. 

The nuclear fragmentation cross sections dis- 
cussed herein axe approximated as the abrasion- 
ablation model of Bowman, Swiatecki, and Tsang 
(ref. 4). The cross section for removal of A A nu- 
cleons is estimated as 

cr(A A) = 7T&2 — nbl (27) 

where b 2 is the impact parameter for which the 
volume of intersection of the projectile contains A a b r 
nucleons and the resulting excitation energies release 
A a bl additional nucleons at the rate of 1 nucleon for 
every 10 MeV of excitation so that 

A a br ( b 2) + A ab l (^> 2 ) = AA - 1 (28) 

and similarly for b\ 

Aabr(f>l) + A abl (&i) = AA+ 1 - (29) 

The charge distributions of the final projectile frag- 
ments are strongly affected by nuclear stability. We 
expect that the charge distribution given by Rudstam 
(ref. 35) for a given <j(AA) to be reasonably correct 
as 


Obviously, Etot as given by equation (26) is very 
crude and substantial improvements can be made. 
The values of A A for carbon projectiles on a cop- 
per target and for copper projectiles on a carbon 
target are shown in figure 1 for high energies. A 
real collision would be given by a statistical distribu- 
tion between the limits shown by these two curves. 
The average event will be calculated as if the two ex- 
tremes occurred with equal probability as noted in 
equation (23). 


t(A f ,Z f ) = Fi exp \-R\Z f - SA f + TAy\ zn 


a(AA) 

(30) 


where the values of R = 11.8/Aj?, D = 0.45, 5 = 
0.486, and T = 3.8 x 10~ 4 axe taken from Rudstam 
and F\ is a normalizing factor so that 


'Y1 CT ( A F’ Z f) = Aj4) (31) 
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The formula of Rudstam for a(AA) was not used 
because his AA dependence is too simple and is not 
useful for heavy targets. For fragments with a mass 
of 9, the cross sections were reduced by a factor of 
2.6, and for fragments with a mass of 7, cross sections 
were increased by a factor of 1.25 and distributed 
between 7 Li and 7 Be with factors of 0.52 and 0.48, 
respectively. 

The charge of the removed nucleons AZ is calcu- 
lated according to charge conservation 

A Z=Zp-Z F (32) 

where Zp is the charge of the fragment and Zp is 
the charge of the projectile. The charge is divided 
according to the following rules among the nucleons 
and the light nuclei produced by the interaction. The 
abraded nucleons are those removed from the por- 
tion of projectile in the overlap region. Therefore, 
the abraded nucleon charge is assumed to be propor- 
tional to the charge fraction of the projectile nucleus 
and is given as 


^abr = ZpA^/Ap (33) 

This assumption, of course, ignores the charge sep- 
aration due to the giant dipole resonance model of 
Morrissey et al. (ref. 33). The charge release in the 
ablation is then given as 


Z a.b\ = AZ - Z abr (34) 

which conserves the remaining charge. Similarly, 
mass is conserved to obtain A^bi- 

The alpha particle is unusually tightly bound in 
comparison to other nucleon arrangements. Because 
of this unusually tight binding, the helium produc- 
tion is maximized in the ablation process as 


where the term i ranges over the mass numbers 2, 3, 
and 4 for ablated particles. 

The calculation is performed for A A = 1 to 
A A = Ap — 1 for which the cross section associated 
with A.4 > Ap — 0.5 is missed. This missed region 
corresponds to the central collisions for which it is 
assumed that the projectile disintegrates into single 
nucleons if rp < rj then 

N p = Z P (38) 

and 


N n — Ap — Zp (39) 

Otherwise, this missed region is ignored. The ener- 
getic target fragments are being ignored as well as the 
mesonic components. The peripheral collisions with 
A A < 0.5 are also missing. The most important pro- 
cess in these collisions with large impact separations 
will be the coulomb dissociation process (ref. 30). 

Electromagnetic Dissociation 

The total electromagnetic cross section for one 
nucleon removal that results from electric dipole (El) 
and electric quadrupole (E2) interaction is written as 


0em — &E\ + a E2 

= f [ N E1 ( E )°E\ i E ) + N E2 i E ) °E2 (-£)] dE 

(40) 

where the virtual photon spectra of energy E pro- 
duced by the target nucleus are given by 


n ei{ e) = L 2 -z^1, 


£KqKi - -Z 2 3 2 ( K 2 - K 2 ) 


(41) 


for the dipole field and by 


N a — [Int(Z a bi/2), Int(>l a bi/ 4)] m i n i muni (35) 

where Int(x) denotes the integer part of x. The 
other light isotopes are likewise maximized from the 
remaining ablated mass and charge numbers in the 
order of decreasing binding energy. The number of 
protons produced is given by charge conservation as 

Np = Z abl - Zi^i ( 36 ) 

i 

Similarly, mass conservation requires the number of 
neutrons produced to be 

Nn = A ahl -N p -J2AN i (37) 


"«<*>- i l z2a w 


2(1 - / 3 2 )K 2 


+ £(2 - 0 2 ) 2 K 0 K X - l - i 2 0 2 (K 2 - K 2 ) 


(42) 


where a is the fine structure constant for the 
quadrupole field. (See ref. 36.) The terms api(E) 
and cte2(E) ar ^ the corresponding photonuclear re- 
action cross sections for the fragmenting projectile 
nucleus. The terms Ko and K\ in the expression for 
Npi and Np 2 are modified Bessel functions of the 
second kind and are also functions of the parameter 
£ with 


* = 


27 iEb„ 


y(3hc 


(43) 
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where E is the virtual photon energy, b^n is the 
minimum impact parameter below which the colli- 
sion dynamics are dominated by nuclear interactions 
(rather than electromagnetic (EM) interactions), /? is 
the speed of the taxget (measured from the projectile 
rest frame) as a fraction of the speed of light c, Ti is 
Planck’s constant, and 7 is the Lorentz factor from 
special relativity that is given by 7 = (1 — /? 2 ) -1 / 2 . 
The minimum impact parameter is given by 


In equations (48) and {49), Np is the number 
of neutrons, Zp is the number of protons, and Ap 
is the mass number of the projectile nucleus. The 
fractional exhaustion of the energy weighted sum rule 
in equation (49) is (ref. 34) 

{ 0.9 {A P > 100) } 

0.6 (40 < A P < 100) \ (50) 

0.3 (40 < A P ) J 


^min “ (1 + x d) b c + 


7rao 


where = 0.25 and 


(44) 


<20 


ZpZpe 2 
mo/? 2 c 2 


(45) 


allows for deviation of the trajectory from a straight 
line (ref. 37 ). The critical impact parameter for single 
nucleon removal is 

be = 1.34 [>lp /3 + - 0.75 (aJ/ S + Ap /3 )] (46) 


where b c is in units of fm and Ap and Ap are 
the projectile and the target nucleon mass numbers, 
respectively. 

The photonuclear cross sections <T£i(E) and 
<j E2 {E) are Lorentzian shaped and somewhat sharply 
peaked in energy. Therefore, the photon spectral 
functions can be taken outside the integral of equa- 
tion (40) to yield an approximate form given by 
(ref. 36) 


^em ~ ^Ei(EgDk) J °E\ ( E ) dE 

+ N E2{ E GQr) E GQR J <J E2i E ) ~2 ( 47 ) 

where jEgdr an d E GQ R are the energies at the 
peaks of the giant dipole resonance (GDR) and gi- 
ant quadrupole resonance (GQR) photonuclear cross 
sections, respectively. These integrals of photonu- 
clear cross sections over energy are evaluated with 
the following sum rules: 


In equation (47) £?GDR £gqr are the energies 
at the peaks of the GDR and GQR photonuclear 
cross sections. For the dipole term it is 


he 


( 1 4 * 6 + 3 U \ 

1 + U 6 ) 

27T 

87 

V 1 + £ + U J 


which is expressed in units of MeV with 


- 1/2 

(51) 


and 


u - IL A~ l/Z 
u - Q' Ap 

Ro = r 0 A]i 3 


(52) 


(53) 


where e = 0.0768, Q* = 17 MeV, J = 36.8 MeV, 
ro = 1.18 fm, and m* is 7/10 of the nucleon mass. 
(See ref. 34.) For the quadrupole term, it is simply 
given by 

E GQR = ^3 (54) 


which is expressed in units of MeV. 


Finally, the single proton or single neutron re- 
moval cross sections are obtained from a em (eq. (47)) 
with proton and neutron branching ratios g p and g n , 
respectively, as 


a(i) = (t = P or n) (55) 


The proton branching ratio has been parameterized 
by Westfall et al. as 


g p = min 


1.95 exp(-0.075Z P ) 

Ap 


(56) 


J o E \(E) dE = 60 


NpZp 

Ap 


which is expressed in units of MeV-mb and 


/ I TH 

0E2{E)e 2 = °- 22 f Z P A T 


which is expressed in units of /ib/MeV. (See ref. 


w r here Zp is the number of protons, and the mini- 

(48) 

mum value of the two quantities in square brackets 
is taken. (See ref. 34.) This parameterization is satis- 
factory for heavier nuclei {Zp > 14). For light nuclei, 
however, the following branching ratios are used: 


(49) 

( 0.5 

(Z P < 6 ) "j 



9p = 1 0.6 

(6 < Zp < 8 ) > 

(57) 

36.) 

(0.7 

(8 < Zp < 14) J 
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For neutrons, the branching ratio is given by 

9n = 1 — Qp (58) 


Database Evaluation 

Measurements have been made for carbon ion 
beams on carbon target at the four energies of 
250 A MeV (ref. 38), 600 >1 MeV (ref. 22), 1.05 A GeV 
(ref. 39), and 2.1 A GeV (ref. 39) and are shown in 
figure 2 with results from NUCFRG2. These frag- 
mentation cross sections are among the best known 
(ref. 40). The effects of the coulomb trajectory are 
clearly apparent in the energy dependence of the 
lighter mass fragment cross sections of Li and Be 
below 100 A MeV. These coulomb effects will be 
even more important for projectiles and targets with 
greater charge. Figure 3 shows the NUCFRG2 model 
at very low energy (11.7 A MeV) for 16 0 projectiles 
onto an 92 Mo target where coulomb trajectory cor- 
rections are very important (ref. 40). The cross sec- 
tions of the resulting charge removal seem well repre- 
sented by NUCFRG2 even at such low energies. The 
addition of exchange poles to the model would bring 
the cross sections of A Z = 0 into agreement, which 
can be judged by the proton exchange pole contribu- 
tion for A Z = — 1 as shown in figure 3. Clearly, the 
model gives a far better result than expected. 

There are three projectile and target combina- 
tions for which two groups, 1.55 A GeV (ref. 20) and 
1.88 A GeV (ref. 34) have measured cross sections at 
nearly the same energy. On the basis of NUCFRG2, 
very small cross section differences are expected at 
these energies. (See fig. 2.) The cross sections from 
NUCFRG2 tend to agree more closely with the ex- 
periments of Westfall et al. (ref. 34) and are 10 to 
50 percent higher than those measured by Cummings 
et al. (ref. 20) as shown in figures 4 to 6. However, 
for charges between 10 to 13, NUCFRG2 agrees more 
closely with the data of Cummings et al. 

To better quantify the comparison of results 
shown in figures 4 to 6, a chi square analysis is used. 
A comparison of NUCFRG2 model with the exper- 
iments of Westfall et al. (ref. 34) and Cummings 
et al. (ref. 20) is shown in table 1 for iron projec- 
tiles on three targets. Shown in table 1 are the total 
chi square value and the average chi square contribu- 
tion per degree of freedom n. Clearly, the data for 
producing .41 fragments in the Westfall et al. exper- 
iments show large systematic errors and is the dom- 
inant contribution to the chi square value. Except 
for the Al datum, the model shows good agreement 
with the data of Westfall et al. for carbon and copper 



Figure 2. Charge removal cross sections for carbon ions on a 
carbon target as a function of projectile energy. 



Figure 3. Charge- removal cross sections for 11.7 A MeV 16 O 
projectiles onto 92 Mo targets. 

targets. The greater discrepancy for the lead targets 
surely results from simplified nuclear matter distri- 
bution in NUCFRG2. A diffuse model instead of 
the uniform spheres of the NUCFRG2 computation 
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o Experimental data at 1.55 A GeV (ref. 20) 
+ Experimental data at 1.88 A GeV (ref. 34) 



Figure 4. Charge removal cross sections for Fe ions on a 
carbon target compared with measurements by Westfall 
et al. and Cummings et al. 


o Experimental data at 1.55 A GeV (ref. 20) 
+ Experimental data at 1.88 A GeV (ref. 34) 



Figure 5. Charge removal cross sections for Fe ions on a 
copper target compared with measurements by Westfall 
et al. and Cummings et al. 

is recommended. This growing discrepancy with in- 
creasing target mass leads to a lack of unitarity given 
by the condition 

Ap <?abs = ^ (^9) 


o Experimental data at 1.55 A GeV (ref. 20) 
+ Experimental data at 1.88 A GeV (ref. 34) 



Figure 6. Charge removal cross sections for Fe ions on a lead 

target compared with measurements by Westfall et al. and 

Cummings et al. 

In generating nuclear data for transport studies, a 
correction is applied to targets with a charge num- 
ber greater than 29 to ensure mass conservation. 
Comparing the NUCFRG2 model with the data of 
Cummings et al. shows similar trends with target 
mass, but the overall agreement with the data of 
Cummings et al. is inferior to agreement with the 
data of Westfall et al. as we have noted in the dis- 
cussion of figures 4 to 6. 

The chi square analysis has been used to compare 
how well one experimental group compares with the 
results of another group. These results are shown in 
table 2. From table 2, it is clear that the NUCFRG2 
model better represents the two sets of experimental 
data than either experimental data set represents the 
other. The systematic errors led to a chi square per 
datum of 10 to 15, which might be used as a measure 
of goodness of fit of the NUCFRG2 model. 

The fragmentation cross section for several pro- 
jectile species on carbon targets at several energies 
are shown in table 3 (refs. 19 to 22, 34, 38, and 41). In 
the table, cross sections for the fragment charge, the 
mass for isotopic measurements, the statistical uncer- 
tainty, the results of NUCFRG2, and the chi square 
per datum are shown. If the error in the experimental 
data were only statistical, then a chi square per da- 
tum value of 1 to 2 would be appropriate and a data 
set with near zero would be viewed with suspicion. 
Clearly, large values of chi square per datum indicate 
possible systematic errors or errors in the NUCFRG2 
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model and it is difficult to make a clear judgment 
in all cases. Although the 12 C projectile isotonic 
breakup cross sections fit the data perhaps too well, 
the isotopic distributions of NUCFRG2 are distinctly 
different from the experiments in most cases. 

Unlike the 12 C projectiles, the isotonic and the 
isotopic fragmentations of 16 0 are outside the statis- 
tical uncertainty of the experiments. The 23 Na and 
24 Mg fragmentations show a strong even-odd effect 
on fragment charge, which are not well represented 
in the results of the NUCFRG2 model. The more 
massive projectiles of 40 Ar and 56 Fe are better rep- 
resented by the NUCFRG2 model for the isotope dis- 
tributions. The isotonic fragmentation cross sections 
for 56 Fe on carbon targets are represented well by the 
NUCFRG2 model. 

The isotopic distributions of light fragments 
shown in table 4 are improperly distributed in the 
breakup of 12 C projectiles on numerous targets. In 
particular, the isotopes of B are produced in dispro- 
portionate numbers for all targets. We also note 
that the cross sections from the NUCFRG2 model 
for heavier targets are systematically low and proba- 
bly result from the assumption of uniform spheres in 
the nucleus (ref. 31). This error is corrected in the 
model presented herein by applying a renormaliza- 
tion factor Nr to individual fragment cross sections 

Nr = A p <7 abs / ^2 A j ( 60 ) 

j 

where p denotes the projectile and j the fragment. 
This unitarity factor is used in the final database 
but is not included in tables 1 to 7. 

The light fragment distributions through carbon 
isotopes are examined best by comparing the 16 0 
breakup cross sections. We again suspect signifi- 
cant nuclear structure effects that are not well rep- 
resented by the Rudstam formalism. Heavier targets 
show larger differences with the NUCFRG2 model 
presented herein because of effects of the diffuseness 
at the nuclear surface. The poor representation of 
the isotopic distribution of light fragments persists 
in the 42 Ar fragmentation in KC1 shown in table 6. 

The remaining 56 Fe fragmentation data in table 7 
are the measurements of Westfall et al. (ref. 34) an d 
Cummings et al. (ref. 20). Generally the NUCFRG2 
model agrees with the two experiments (at least to 
the degree that they agree with each other) except for 
the few spurious data points in the data of Westfall 
et al. The model is most accurate for light targets. 

The variation of chi square per n over the avail- 
able experimental data is summarized in table 8. Re- 


call that the estimate of systematic experimental er- 
ror obtained by comparing the Westfall et al. iron 
data with that of Cummings et al. gives the experi- 
mental chi square per n of 4 to 12. The corresponding 
model chi square per n for NUCFRG2 for these ex- 
perimental data sets is from 2 to 7, which shows that 
NUCFRG2 tends to split the difference between the 
two experiments. The iron fragmentation for targets 
below suffer show comparable model chi square per n 
between 2.3 and 7.4. The model chi square per n for 
targets above suffer indicate systematic model errors, 
which need to be resolved. Clearly, the current ex- 
perimental data are adequate as a measure of model 
improvement. However, current experimental accu- 
racy is inadequate to evaluate the resulting nuclear 
database to the accuracy required for shield design. 

More recently, fragmentation cross sections were 
measured at the Lawrence Berkeley Laboratory for 
600 A MeV iron beams on 2-cm-thick polyethylene 
(CH2 ) x targets (ref. 41). The results are shown 
in figure 7. A systematic error was introduced by 
an electronic trigger inefficiency to fragments lighter 
than Ne. The chi square per datum for fragments 
heavier than Ne is 3 x 10~ 4 , while for the complete 
data set the chi square per datum is 6 x 10 -4 . 


O NUCFRG2 



Figure 7. Charge removal cross sections for Fe ions at 600 A 

MeV in polyethylene. 

The comparison of the NUCFRG2 model with the 
measured fragmentation cross sections mainly tests 
the dependence of a(AA) and the applicability of 
the formalism of Rudstam for the charge distribution. 
A more sensitive test of the model representation of 
the ablation process is to compare the numbers and 
the types of particles produced. The multiplicities 
of the charged reaction products were measured by 
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Figure 8. Relative probabilities for multiplicity of various charge states for ^Fe ion beams in nuclear emulsion. 


Dudkin et al. (ref. 42) and are shown in figure 8 for 
comparison with calculations with the fireball model 
(ref. 43) and the version of NUCFRG2 presented 
herein. The primary difference between the two mod- 
els is the semiempirical correction to the excitation 
energy required by the NUCFRG2 model to fit the 
atmospheric air shower data (ref. 42) and the imple- 
mentation of the deexcitation process. 

As shown by Dudkin et al., the fireball model 
shows even qualitative differences in the frequency 
distribution of multiplicities in nuclear emulsion 
when compared with the experimental results 


(ref. 42). In general, the fireball model over- 
estimates the events of low multiplicity and cor- 
respondingly underestimates the high multiplicity 
events. The NUCFRG2 model gives a much im- 
proved distribution of events, although the predicted 
number of high multiplicity events appears greater 
than is seen experimentally. However, the resolution 
of the multiplicity when many high-energy secon- 
daries are produced at one apex is not good. Clearly, 
if some of the observed several prong events were in 
fact of higher multiplicity, then good agreement be- 
tween the NUCFRG2 model and the experiments is 
conceivable. 
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Concluding Remarks 

The analysis of experimental data of iron beam in- 
teractions in polyethylene and aluminum targets has 
resulted in an improved semiempirical model for nu- 
clear database generation. However, this analysis is 
still unable to unambiguously resolve differences be- 
tween various experimental groups that use the same 
projectile and target combinations. Further improve- 
ments to the basic model have been made as a result 
of the presented evaluation. These improvements are 
correcting the charge distributions in the formulation 
of Rudstam for light nuclei where shell structure ef- 
fects cause important deviations from simple system- 
atic behavior. Even with these improvements, there 
is need for improved target mass dependence where 
the skin diffuseness is suspected to require further 
correction. 

There are two approaches to improve the gen- 
eration code of the NUCFRG2 nuclear database. 
The first is incremented improvements to the semi- 
empirical model. The second approach is more radi- 
cal and involves the further development of quantum 
based methods. Ultimately, model development is 
still limited by the systematic errors in the experi- 
mental fragmentation data and the paucity of exper- 
imental data. 

The semiempirical model is first limited by the 
assumed uniform nuclear matter distribution. This 
assumption is the main source of nonconservation 
of mass and charge for massive targets that is tem- 
porarily corrected by forcing unitarity; however, a 
fully correct description must replace the uniform 
sphere model with realistic nuclear density distribu- 
tions. Although ad hoc corrections have been made 
for structure effects in the low mass fragments (mass 
number less than 10), errors remain in the isotope 
distributions lighter than Ne. Improved methods 
for representing the distribution of excitation energy 
and corrections to Rudstam’s distribution will be re- 
quired. The distribution of mass and charge in the 
final ablation products depends on the excitation 
energy. Although the multiplicities are reasonable 
and are greatly improved over the fireball model, the 
agreement with atmospheric air shower data might 
be improved. 

The microscopic description of nuclear fragmen- 
tation proceeds from a multiple-scattering theory 
(MST), a description that uses the Glauber or 
eikonal approximation of the multiple-scattering se- 
ries, which results from a systematic reduction of rel- 
ativistic MST or a nonrelativistic MST. A relativistic 
or nonrelativistic model may be cast as a distorted- 


wave series by using an average optical potential for 
elastic transitions. 

A very difficult task still remains in treating mul- 
tiple inelastic transitions that occur in heavy ion frag- 
mentation. The difficulty arises first because of the 
many irreducible diagrams that contribute to any 
integral equation that would be formulated. Typi- 
cally, integral equations of six or more dimensions 
will occur for inelastic transitions. This occurence 
of six or more dimensions is in comparison with the 
3-dimensional integral equations that appear in the 
optical model formulation of elastic heavy ion scat- 
tering. Also, a perturbative approach becomes diffi- 
cult because of the large number of terms required for 
heavy ion scattering and the complexity of summing 
over intermediate state variables. 

A more practical approach is to use the forward 
scattering assumptions of the eikonal model, which is 
expected to be valid at high energies. This approach 
allows for closed-form expressions to be derived for 
the multistep processes. Here, the eikonal approx- 
imation can be applied in both the nonrelativistic 
MST or a relativistic MST. A relativistic eikonal 
model could allow the effects of negative-energy 
states on nuclear fragmentation to be considered. A 
study of the one- and two-step contributions to frag- 
mentation could be made in the relativistic or non- 
relativistic models to provide a test of the validity 
of the eikonal approximation and the neglecting of 
off-shell effects. In all of the approaches mentioned, 
a treatment of meson and antinucleon production in 
the MST should be considered. 

An alternative to the use of an MST or the 
Glauber model in formulating the fragmentation 
process is to derive transport equations from an 
MST. In the MST, it is difficult to include the 
cascade of projectile knockouts through the pro- 
jectile prefragments. In nucleon-nucleus scattering, 
the number of cascade particles is small at low to 
medium energies that are typical of the relative en- 
ergy expected between knockouts and prefragments. 
However, the multiple-scattering cascade terms will 
depend strongly on the initial trajectory of the 
knockouts, which will lead to a heavy computational 
burden. A derivation of the simple final state in- 
teraction (FSI) corrections may help to simplify an 
MST approach. The transport models are a more 
tractable approach to the cascade problem. These 
models are usually used to study the knockout spec- 
trum of nucleons. In some cases, the spectrum of 
residual energies after the cascade is used to predict 
the final fragment mass yields. The reliance upon 
classical methods that ignore quantum effects and 
the use of Monte-Carlo simulations that require large 
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computational times may limit the usefulness of this 
approach. 

In describing the abrasion- ablation process, the 
dynamical model of choice should be used to formu- 
late the creation of the prefragment state in an arbi- 
trary configuration. The description of the prefrag- 
ment configuration requires the variables mass and 
charge number, the distribution of excitation ener- 
gies, and the spin to be complete. This description 
should distinguish whether the creation of the pre- 
fragment occured after multistep nucleon removal or 
cluster removal and predict the correct yields and the 
energy spectrum of light fragments in the reaction. 
Such a description would require a large amount of 
nuclear structure input. For example, the study of 
cluster abrasion requires a more detailed description 
of the nuclear ground state than the independent par- 
ticle model. The development of the formalism for 


overlap function of the nuclear ground state for com- 
peting cluster configurations is required to aid the 
description of the abrasion process. Clearly, nuclear 
shell effects will be more correctly described in such 
a formalism. Much information on the level densities 
and decay modes of nuclei is well known, although 
the treatment of extremely high excitation energies 
is less understood. Also less understood is whether 
there is a one-to-one correspondence between the de- 
cay of excited prefragments produced in heavy ion 
fragmentation and the equilibrium decay treated by 
statistical methods. Further studies in these areas 
should be made. 


NASA Langley Research Center 
Hampton, VA 23681-0001 
July 7, 1995 
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Appendix 

Abrasion Model Formula 

Two functions F and P used in the formalism given in the text are described in this appendix. For rp > rp, 
where rp is target radius and rp is projectile radius, we have 


P = 0.125(/xi/) 1 / 2 fi-2 


1-0 


0.125 


0.5(pi/) 1/2 f i-2j +1 


1-0 


F = 0.75(1 - i/) 1 / 2 - 0.125(3(1 - v) 1 ' 2 - 1] (!—£)' 


(Al) 

(A2) 

(A3) 
(A4) 

li - (l/i/) - 1 = rj/rp (A5) 

(See ref. 32.) Equations (Al) and (A2) are valid when the collision is peripheral (i.e., the two nuclear volumes 
do not completely overlap). In this case, the impact separation r is restricted so that 


and 


with 


and 


V = rp/(rp + r T ) 
0 = b/ (r P + rp) 


rp — Tp < r < rp + rp 


(A6) 


If the collision is central, then the projectile nucleus volume completely overlaps the target nucleus volume 
{r <rp - rp), and all the projectile nucleons are abraded. In this case, equations (Al) and (A2) are replaced 
by 

P = -l (A7) 

and 

F = 1 (A8) 

and there is no ablation of the projectile, because it was destroyed by the abrasion. 

For the case where rp > rp and the collision is peripheral, equations (Al) and (A2) become 


and 


P = 0.125(/iJ/) 1 ^ 2 (l — 2^ 

-[(1/*0(1-V) 1/2 - 



F = 0.75(1 - i/) 1/2 (h_£) 2 - 0.125 j 3 ^ - ^ )1/2 

_ n-(i-A* 2 ) 3 / 2 ][i-(i-ii) 2 i 1 ^ 2 j ( i-/? ) 


where the impact separation is restricted so that 


(A9) 


(A10) 


14 


r p -r T <r <r p +r T 


(All) 



(A12) 


(See ref. 33.) For a central collision (r < r p - r T ) with r p > r T , equations (A9) and (A10) become 


and 


P = 


-d-^) 


2 N 1/2 



F = 


i-(i-^ 2 ) 3/2 ] 



(A13) 
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Table 1 . Chi Square Analysis of Iron Fragmentation Model 



NTJCFRG2 and Westfall et al. 

NUCFRG2 and Cummings et al. 

System 

x 2 

X 2 without Al data 

x 2 /« 

X 2 ln without Al data 

x 2 

X 2 /n 

Fe + C 

50.2 

16.0 

5 

1.8 

48.3 

3.7 

Fe + Cu 

200.6 

22.9 

20 

2.5 

78.4 

6.0 

Fe + Pb 

177.4 

56.2 

18 

6.2 

83.1 

6.7 


Table 2. Chi Square Analysis of Iron Fragmentation Experiments 



Westfall et al. and Cummings et al. 

Cummings et al. and Westfall et al. 

System 

x 2 

% 2 without Al data 

X 2 /" 

X 2 /n without Al data 

x 2 



X 2 /n without Al data 

Fe + C 

85.6 


8.6 

4.8 

54.6 

33.6 

5.5 

3.7 

Fe + Cu 

424.4 


42.4 

12.0 

160.3 

69.4 

16.0 

7.7 

Fe + Pb 

348.8 

79.5 

34.9 

8.8 

143.1 

55.8 

14.3 

6.2 


18 





























Fragment 

Projectile charge 
n B 

326 A MeV a 4 

3 

,2 C 

250 A MeV b 6 

6 
5 

5 

4 
4 
4 
3 
3 

3 

12 C 

403 A MeV a 5 

4 

12 C 

418 A MeV a 5 

4 

12 C 

561 A MeV a 5 

4 

12 C 

600 A MeV c 6 

6 

5 
5 
4 
4 
4 

12 C 

693 A MeV a 5 

4 


Fragment 
mass number 

Experiment cross 
section, mb 

Experiment 
uncertainty, mb 

NUCFRG2 
cross section, mb 

Chi square per 
degree of 
freedom (n) 


105.9 

1.59 

99.31 



30.1 

.45 

56.84 






6.508 

11 

55.97 

4.06 

56.19 


10 

5.33 

.81 

.37 


11 

65.61 

2.55 

56.33 


10 

47.50 

2.42 

57.46 



5.88 

9.70 

3.47 


9 

10.44 

.85 

14.21 


7 

22.64 

1.49 

20.37 


8 

1.33 

1.00 

.11 


7 

17.19 

3.00 

20.07 


6 

26.35 

2.10 

30.80 






8.504 


106.0 

1.59 

1 14.62 



29.6 

.89 

40.79 






.929 


111.2 

1.67 

114.50 



32.1 

.96 

40.66 






.949 


108.7 

1.63 

113.07 



30.3 

.91 

40.16 






1.295 

11 

53.6 

.80 

54.1 


10 

2.1 

.11 

.3 


11 

70.7 

1.06 

54.3 


10 

38.6 

.58 

55.3 


10 

5.6 

.28 

3.3 


9 

9.6 

.29 

13.52 


7 

15.5 

.47 

19.53 






3.092 


110.1 

1.65 

112.41 



34.9 

1.05 

39.76 






0.321 













Table 3. Continued 


Projectile 

Fragment 

charge 

Fragment 
mass number 

Experiment cross 
section, mb 

Experiment 
uncertainty, mb 

NUCFRG2 
cross section, mb 

Chi square per 
degree of 
freedom ( n ) 

12 c 

915 AMeV a 

5 


109.4 

1.64 

111.1 



4 


33.6 

1.01 

39.0 


12 C 

1016 AMeV a 

5 


113.2 

1.70 

110.8 

0.387 


4 


36.8 

1.10 

39.0 


12 C 

1050 AMeV d 

6 

ii 

44.70 

2.80 

52.7 

.072 


6 

10 

4.44 

.24 

.4 



5 

11 

48.60 

2.40 

52.9 



5 

10 

27.90 

2.20 

54.4 



4 

10 

5.34 

.29 

3.2 



4 

9 

10.70 

.50 

13.1 



4 

7 

18.60 

.90 

19.0 



3 

8 

2.40 

.18 

.1 



3 

7 

21.50 

1.10 

20.5 



3 

6 

27.10 

2.20 

28.5 



2 

6 

1.83 

.19 

.4 


12 C 

1572 A MeV a 

5 


103.9 

1.56 

109.2 

10.836 


4 


35.6 

1.07 

38.1 


12 C 

2 100 A Me V d 

6 

11 

46.50 

2.30 

51.3 

.211 


6 

10 

4.11 

.22 

.3 



5 

11 

53.80 

2.70 

51.6 



5 

10 

35.10 

3.40 

53.7 



4 

10 

5.81 

.29 

3.2 



4 

9 

10.63 

.53 

12.7 



4 

7 

18.61 

.93 

18.5 



3 

8 

2.18 

.15 

.1 



3 

7 

21.50 

1.10 

20.0 



3 

6 

30.00 

2.40 

27.9 



2 

6 

2.21 

.22 

.4 


,4 N 

516 AMeV a 

6 


169.2 

2.54 

124.0 

9.456 


5 


63.1 

1.89 

85.8 


. 

4 


27.9 

.84 

33.8 

7.828 


20 






Projectile 


Fragment 

charge 


Fragment 
mass number 


Experiment cross 
section, mb 


16 , 


O 


441 A MeV a 


16 q 

491 A MeV a 


16 , 


0 


669 A MeV a 


16 , 


o 


903 A MeV a 


16 , 


o 


1563 AMeV a 


16 i 


o 


2100 AMeV d 


7 

6 

5 

4 


7 

6 

5 

4 


7 

6 

5 

4 


7 

6 

5 

4 


7 

6 

5 

4 


8 

8 

7 

7 

7 

7 

6 

6 

6 

6 

6 

5 

5 

5 


15 

14 

15 
14 

13 
12 

14 
13 
12 
11 
10 
13 
12 
11 


162.9 

160.2 

60.7 

13.6 


146.4 

146.2 

54.7 

13.4 


158.5 

159.6 
56.5 
17.3 


154.4 

152.9 

52.4 

20.3 


125.3 

123.2 

46.6 

18.2 


42.90 

1.67 

54.20 

41.80 

8.06 

.73 

4.71 

27.70 

65.10 

18.46 

2.51 

.44 

2.44 

26.0 


Experiment 
uncertainty, mb 

NUCFRG2 
cross section, mb 

Chi square per 
degree of 
freedom ( n ) 

2.44 

133.8 


2.40 

111.3 


1.82 

75.7 


.65 

32.9 

9.718 

2.20 

133.3 


2.19 

110.6 


1.64 

75.5 


.67 

32.8 

7.491 

2.38 

131.7 


2.39 

109.3 


1.69 

74.7 


.87 

32.3 

10.661 

2.32 

130.4 


2.29 

108.2 


1.57 

73.9 


1.02 

32.1 

8.331 

1.88 

128.4 


1.85 

106.1 


1.40 

72.8 


.91 

31.7 

4.500 

2.30 

57.6 


.12 

.6 


2.90 

58.0 


3.30 

61.7 


.42 

6.9 


.07 

.3 


.31 

10.8 


1.40 

47.7 


5.20 

40.3 


.92 

6.0 


.16 

.2 


.05 

.4 


2.15 

4.3 


1.30 

32.7 













Table 3. Continued 


Projectile 


2100 A MeV 


20 Ne 

468 A MeV a 


20 Ne 

599 A MeV a 


20 Ne 

608 A MeV a 


20 Ne 

1057 A Me V £ 


2J Na 

461 A MeV a 


Chi square per 

Fragment Fragment Experiment cross Experiment NUCFRG2 degree of 
charge mass number section, mb uncertainty, mb cross section, mb freedom ( n ) 














Chi square per 

Fragment Fragment Experiment cross Experiment NUCFRG2 degree of 
charge mass number section, mb uncertainty, mb cross section, mb freedom (n) 


"Mg 

309 A MeV a 


"Mg 

48! A MeV a 


"Mg 

739 A MeV a 


"Mg 

1455 A MeV a 


ZI A\ 

582 A MeV a 
















Projectile 

Fragment 

charge 

Fragment 
mass number 

Experiment cross 
section, mb 

Experiment 
uncertainty, mb 

NUCFRG2 
cross section, mb 

40 Ar 

1650 A Me V e 

17 

36 

49.50 

16.60 

47.71 


17 

35 

51.00 

19.50 

27.34 


17 

34 

12.00 

6.45 

4.93 


16 

38 

4.35 

1.50 

.15 


16 

37 

11.70 

2.70 

1.19 


16 

36 

12.40 

3.00 

4.48 


16 

35 

24.00 

5.40 

21.83 


16 

34 

49.50 

10.80 

38.49 


16 

33 

31.50 

7.05 

27.31 


16 

32 

10.60 

3.00 

5.74 


16 

31 

.54 

.42 

.99 


15 

36 

.615 

.195 

.06 


15 

35 

2.10 

.33 

.58 


15 

34 

5.85 

.825 

2.67 


15 

33 

18.00 

1.80 

14.36 


15 

32 

27.00 

2.70 

32.65 


15 

31 

21.00 

2.10 

27.26 


15 

30 

3.90 

.63 

6.83 


15 

29 

.315 

.195 

1.05 


14 

34 

.01 

.07 

.02 


14 

33 

1.32 

.24 

.30 


14 

32 

3.00 

.24 

1.71 


14 

31 

11.00 

1.80 

9.40 


14 

30 

37.50 

3.15 

27.76 


14 

29 

25.50 

3.00 

26.96 


14 

28 

13.00 

1.95 

8.17 


14 

27 

.69 

.285 

1.11 


13 

32 

.129 

.111 

.01 


13 

31 

.705 

.21 

.15 


13 

30 

3.00 

.405 

1.10 


13 

29 

10.40 

2.25 

6.08 


13 

28 

19.50 

2.40 

23.48 


13 

27 

25.50 

2.85 

26.25 


13 

26 

7.20 

1.28 

9.85 


13 

25 

.315 

.165 

1.17 


12 

30 

.165 

.126 

.01 


12 

29 

.66 

.21 

.01 


12 

28 

2.10 

.495 

.71 


12 

27 

6.75 

1.02 

3.86 


12 

26 

24.00 

2.85 

19.77 


12 

25 

22.50 

3.45 

25.33 


12 

24 

14.20 

1.80 

11.46 


12 

23 

.96 

.315 

1.25 


11 

27 

.36 

.135 

.03 


Chi square per 
degree of 
freedom (n) 



























Chi square per 


Fragment 

Fragment 

Experiment cross 

Experiment 

NUCFRG2 

degree of 

Projectile 

charge 

mass number 

section, mb 

uncertainty, mb 

cross section, mb 

freedom (n) 

40 At 

1650 A Me V c 

5 

8 

0.21 

0.255 

0.01 

5.313 

56 Fe 

330 A MeV a 

25 


244.3 

7.33 

196.8 



24 


182.7 

5.48 

132.3 



23 


121.1 

3.63 

107.7 



22 


110.7 

3.32 

92.2 



21 


89.2 

2.68 

81.4 



20 


79.6 

2.39 

73.4 



19 


51.6 

2.58 

67.2 



18 


44.4 

2.22 

62.0 



17 


38.2 

1.91 

57.5 



16 


42.4 

2.12 

53.5 

5.471 

56 Fe 

434 A MeV a 

25 


223.8 

3.36 

194.8 



24 


175.1 

2.63 

131.2 



23 


116.1 

3.48 

106.7 



22 


116.0 

3.48 

91.4 



21 


79.8 

2.39 

80.8 



20 


73.9 

2.22 

73.0 



19 


52.3 

2.62 

66.8 



18 


48.2 

2.41 

61.7 



17 


39.5 

1.98 

57.3 



16 


44.3 

2.22 

53.4 



15 


24.89 

2.49 

49.9 



14 


45.7 

4.57 

47.0 

4.355 

56 Fe 

520 A MeV a 

25 


206.2 

3.09 

193.6 



24 


163.4 

2.45 

130.2 



23 


115.2 

3.46 

105.9 



22 


112.3 

3.37 

91.2 



21 


76.7 

2.30 

80.4 



20 


71.3 

2.14 

72.7 



19 


51.8 

2.59 

66.7 



18 


52.8 

2.64 

61.5 



17 


40.6 

2.03 

57.3 



16 


41.7 

2.09 

53.3 



15 


29.5 

1.48 

49.9 



14 


48.6 

4.86 

46.9 

2.960 








Chi square per 
degree of 
freedom {n) 






Chi square per 

Fragment Fragment Experiment cross Experiment NUCFRG2 degree of 
Projectile charge mass number section, mb uncertainty, mb cross section, mb freedom ( n ) 











Table 3. Continued 



3D Fe 

724 A MeV a 


;>0 Fe 

944 A MeV a 


30 Fe 

1086 A MeV a 


3t> Fe 

1409 A MeV a 


I 








Projectile 


56 Fe 

1409 A MeV a 


Chi square per 

Fragment Fragment Experiment cross Experiment NUCFRG2 degree of 
charge mass number section, mb uncertainty, mb cross section, mb freedom (n) 


3D Fe 

1512 A Me V a 


1570 A Me V r 


140.73 

105.33 

79.32 

75.17 

57.29 

63.37 

43.62 

47.65 

41.45 
46.47 

39.45 
50.99 
41.23 

45.45 
35.83 
44.79 

















Target 


'Be 


1050 A Me V a 


'Be 


2100 A MeV‘ 


27 


A1 


1050 AMeV a 


Fragment 

charge 

Fragment 
mass number 

Experiment cross 
section, mb 

Experiment 
uncertainty, mb 

NUCFRC 
cross section 

6 

11 

44.70 

2.60 

52.88 

6 

10 

4.02 

.23 

.35 

5 

11 

50.70 

3.20 

52.99 

5 

10 

28.80 

2.30 

55.01 

4 

10 

5.08 

.39 

3.32 

4 

9 

11.60 

.76 

13.14 

4 

7 

17.80 

.90 

19.03 

3 

9 

.75 

.08 

.03 

3 

8 

2.36 

.14 

.10 

3 

7 

23.40 

1.20 

20.61 

3 

6 

24.80 

2.00 

29.03 

2 

6 

2.09 

.17 

.37 

6 

| 

1 

11 

46.70 

2.30 

51.51 

6 


4.20 

.21 

.34 

5 


53.20 

2.90 

51.64 

5 


31.10 

2.60 

53.81 

4 


5.97 

.31 

3.25 

4 

9 

10.98 

.55 

12.85 

4 

7 

18.91 

.95 

18.66 

3 

9 

.92 

.08 

.03 

3 

8 

2.52 

.16 

.10 

3 

7 

22.80 

1.10 

20.23 

3 

6 

33.10 

2.70 

28.37 

2 

6 

2.54 

.25 

.36 

6 

11 

57.80 

3.90 

59.27 

6 

10 

5.06 

.37 

.38 

5 

11 

64.50 

5.30 

60.09 

5 

10 

30.40 

3.50 

59.66 

4 

10 

6.49 

.48 

3.60 

4 

9 

13.90 

.90 

14.34 

4 

7 

19.90 

1.10 

20.80 

3 

9 

.82 

.16 

.04 

3 

8 

2.87 

.27 

.11 

3 

7 

28.50 

1.40 

22.53 

3 

6 

24.90 

2.90 

31.62 

2 

6 

2.00 

.29 

.40 


I Chi square per| 
degree of 
freedom (rt) 


10.026 


11.964 


13.464 











Table 4. Continued 


Target 

Fragment 

charge 

Fragment 
mass number 

Experiment cross 
section, mb 

Experiment 
uncertainty, mb 

NUCFRG2 
cross section, mb 

Chi square per 
degree of 
freedom (n) 

27 Al 







2100 A MeV a 

6 

11 

59.50 

3.10 

58.69 



6 

10 

4.99 

.34 

.37 



5 

11 

65.20 

4.80 

59.84 



5 

10 

36.40 

4.80 

58.27 



4 

10 

7.02 

.40 

3.52 



4 

9 

12.74 

.71 

14.02 



4 

7 

25.80 

1.30 

20.43 



3 

9 

.88 

.12 

.04 



3 

8 

2.79 

.23 

.11 



3 

7 

27.30 

1.40 

22.14 



3 

6 

36.30 

2.90 

30.64 



2 

6 

2.82 

.27 

.39 








13.846 

63 Cu 







1050 A Me V a 

6 

11 

78.10 

8.10 

70.91 



6 

10 

7.53 

.70 

.43 



5 

11 

80.10 

7.90 

73.88 



5 

10 

36.40 

9.90 

67.08 



4 

10 

7.69 

.61 

4.05 



4 

9 

14.30 

1.20 

16.29 



4 

7 

25.00 

1.90 

23.08 



3 

9 

1.05 

.38 

.04 



3 

8 

3.99 

.70 

.12 



3 

7 

32.60 

1.90 

25.00 



3 

6 

33.10 

6.00 

35.91 



2 

6 

3.01 

.68 

.46 








23.774 

63 Cu 







2100 AMeV a 

6 

11 

81.40 

6.30 

72.59 



6 

10 

5.38 

.55 

.42 



5 

11 

84.40 

9.00 

77.16 



5 

10 

43.70 

9.80 

66.28 



4 

10 

8.57 

.70 

4.00 



4 

9 

16.10 

1.30 

15.88 



4 

7 

33.70 

2.30 

22.85 



3 

9 

1.38 

.36 

.04 



3 

8 

3.89 

.47 

.12 



3 

7 

31.90 

2.30 

24.75 



3 

6 

47.30 

4.50 

35.30 



2 

6 

3.21 

.47 

.45 




. 




23.436 


34 









Chi square per 


Fragment 

Fragment 

Experiment cross 

Experiment 

NUCFRG2 

degree of 

Target 

charge 

mass number 

section, mb 

uncertainty, mb 

cross section, mb 

freedom (n) 

107 Ag 

1050 A Me' V a 

6 

ii 

98.00 

13.00 

83.18 



6 

10 

7.70 

1.00 

.47 



5 

11 

110.00 

15.00 

89.42 



5 

10 

43.00 

12.00 

74.01 



4 

10 

8.40 

1.20 

4.47 



4 

9 

23.70 

2.70 

17.35 



4 

7 

21.60 

2.70 

25.05 



3 

9 

1.15 

.49 

.04 



3 

8 

2.80 

1.20 

.14 



3 

7 

42.10 

3.40 

27.14 



3 

6 

38.10 

7.60 

38.30 



2 

6 

3.60 

1.40 

.49 

19.443 

107 Ag 

2100 AMeV a 

6 

11 

101.90 

9.60 

89.88 



6 

10 

7.03 

.88 

.47 



5 

11 

109.00 

13.00 

100.17 



5 

10 

65.00 

17.00 

73.03 



4 

10 

8.81 

.91 

4.41 



4 

9 

18.60 

1.70 

17.10 



4 

7 

41.20 

3.30 

24.63 



3 

9 

1.20 

.33 

.04 



3 

8 

3.27 

.53 

.13 



3 

7 

40.30 

3.30 

26.68 



3 

6 

46.10 

5.60 

37.56 



2 

6 

3.50 

1.10 

.48 

19.202 

2 °8 pb 

1 050 A MeV a 

6 

11 

128.00 

22.00 

106.46 



6 

10 

10.90 

1.70 

.54 



5 

11 

149.00 

25.00 

120.13 



5 

10 

51.00 

18.00 

83.37 



4 

10 

10.90 

1.80 

5.04 



4 

9 

22.20 

3.70 

19.68 



4 

7 

37.80 

4.70 

28.02 



3 

9 

1.76 

.81 

.05 



3 

8 

4.90 

1.60 

.16 



3 

7 

45.20 

4.80 

30.36 



3 

6 

51.00 

13.00 

42.94 



2 

6 

7.30 

2.70 

.55 

42.293 


















Target 

Fragment 

charge 

Fragment 
mass number 

9 Be 



2100 A MeV a 

8 

15 


8 

14 


7 

15 


7 

14 


7 

13 


7 

12 


6 

14 


6 

13 


6 

12 


6 

11 


6 

10 


5 

13 


5 

12 


5 

11 


5 

10 


4 

10 


4 

9 


4 

7 


3 

7 


3 

6 

27 A1 



2100 AMeV a 

7 

15 


6 

14 


6 

13 


5 

12 


5 

11 


4 

9 


3 

7 

63 Cu 



2100 A Me' V a 

8 

15 


8 

14 


7 

15 


7 

14 


7 

13 


7 

12 


6 

14 


6 

13 


6 

12 


6 

11 


6 

10 


Experiment cross 
section, mb 

Experiment 
uncertainty, mb 

NUCFRG2 
cross section, mb 

Chi square per 
degree of 
freedom (fi) 

43.00 

2.20 

57.60 


1.60 

.10 

.60 


54.10 

2.70 

57.80 


49.50 

4.00 

63.00 


8.01 

.40 

7.04 


.66 

.05 

.33 


5.21 

.30 

11.10 


28.60 

1.40 

48.20 


60.80 

4.90 

40.91 


21.00 

1.10 

6.08 


2.81 

.17 

.21 


.50 

.04 

.41 


2.75 

.15 

4.36 


27.50 

1.40 

32.95 


19.20 

1.50 

33.02 


3.92 

.27 

1.99 


9.79 

.50 

9.91 


22.00 

1.10 

16.83 


27.00 

1.40 

18.24 


33.50 

2.70 

26.60 

5.722 

66.00 

4.30 

69.48 


6.29 

.46 

12.59 


31.40 

2.00 

55.56 


3.61 

.24 

5.03 


31.00 

1.60 

37.64 


11.22 

.68 

11.34 


34.80 

1.80 

21.12 

3.021 

74.00 

7.80 

85.47 


2.14 

.42 

.78 


98.20 

9.80 

92.10 


72.00 

14.00 

82.02 


14.70 

1.60 

9.22 


.42 

.18 

.43 


7.76 

.92 

14.45 


35.80 

3.70 

63.11 


92.00 

14.00 

53.65 


27.00 

2.60 

7.91 


4.45 

.52 

.27 










Table 5. Continued 


Target 


2100 A MeV 


charge 


Fragment 

Experiment cross 

Experiment 

NUCFRG2 

mass number 

section, mb 

uncertainty, mb 

cross section, mb 



0.82 

2.98 

35.90 

35.20 
.30 

6.51 

12.30 

32.00 

3.63 

38.70 

61.20 


0.54 

5.72 

42.91 

42.64 

.21 

2.57 

12.94 

21.79 

.11 

23.60 

33.87 


99.00 

2.20 

121.00 

68.00 

18.60 

1.11 

7.50 

39.40 

104.00 

37.80 

4.20 

.65 

4.04 

43.60 

26.60 
5.65 

13.80 

36.40 

39.80 

49.40 


108.69 

.85 

123.72 

89.42 

10.16 

.47 

15.75 

69.55 

58.14 

8.70 

.30 

.60 

6.20 

47.16 
47.09 

2.84 

14.16 
23.80 
25.79 
38.06 


135.00 

2.80 

202.00 

71.00 

17.00 
12.30 

45.40 

126.00 

36.90 

7.20 


22.00 

1.50 

26.00 

22.00 

3.20 

2.20 
8.30 

25.00 

5.70 

1.40 


162.79 
.96 

199.79 
100.65 

11.48 

17.73 

78.53 

66.34 

9.90 

.34 








Target 

Fragment 

charge 

Fragment 
mass number 

Experiment cross 
section, mb 

Experiment 
uncertainty, mb 

NUCFRG2 
cross section, mb 

Chi square per 
degree of 
freedom (n) 

208 pb 







2100 AMeV a 

5 

13 

0.70 

0.44 

0.68 



5 

12 

3.98 

.75 

7.07 



5 

11 

52.90 

5.90 

53.67 



5 

10 

35.70 

11.00 

52.95 



4 

10 

6.80 

1.10 

3.20 



4 

9 

15.30 

2.10 

16.36 



4 

7 

43.30 

6.40 

27.84 



3 

7 

39.70 

4.30 

30.17 



3 

6 

56.00 

13.00 

42.04 



17.217 













Target 


KC 

1 650 A Me V 


Fragment 
mass number 

Experiment cross 
section, mb 

Experiment 
uncertainty, mb 

NUCFRG2 
cross section, mb 

34 

17.00 

10.00 

6.05 

35 

38.00 

15.00 

33.26 

36 

6.80 

3.80 

57.01 

37 

42.00 

31.00 

41.18 

39 

56.00 

29.00 

87.56 

31 

1.90 

1.50 

1.23 

32 

20.00 

6.30 

7.08 

33 

41.00 

8.50 

33.54 

34 

50.00 

11.00 

47.23 

35 

32.00 

8.60 

26.56 

36 

29.00 

6.00 

5.35 

37 

19.00 

4.10 

1.40 

38 

9.90 

2.70 

.18 

29 

.42 

.33 

1.33 

30 

7.40 

1.70 

8.53 

31 

21.00 

3.80 

33.93 

32 

25.00 

3.00 

40.25 

33 

22.00 

5.70 

17.64 

34 

2.10 

1.50 

3.27 

35 

2.80 

.93 

.71 

36 

.20 

.07 

.07 

27 

1.40 

.60 

1.42 

28 

15.00 

2.80 

10.41 

29 

38.00 

5.20 

33.97 

30 

43.00 

5.20 

34.70 

31 

14.00 

9.00 

11.70 

32 

1.50 

1.80 

2.11 

33 

.54 

.89 

.37 

34 

.16 

.10 

.03 

25 

.86 

.44 

1.53 

26 

8.00 

1.30 

12.57 

27 

37.00 

4.30 

33.54 

28 

18.00 

4.70 

29.91 

29 

22.00 

5.30 

7.66 

30 

1.50 

.49 

1.38 

31 

.58 

.32 

.19 

23 

.40 

.18 

1.62 

24 

21.00 

3.50 

14.83 

25 

26.00 

3.90 

33.08 

26 

29.00 

3.20 

25.24 

27 

10.00 

4.70 

4.93 

28 

.42 

.80 

.90 









Table 6. Concluded 


Target 


1650 A Me V 


Fragment 

charge 


Fragment 
mass number 

Experiment cross 
section, mb 

Experiment 
uncertainty, mb 

NUCFRG2 
cross section, mb 

Chi square per 
degree of 
freedom (n) 

29 

0.68 

0.34 

0.09 


21 

.76 

.80 

1.84 


22 

12.00 

2.20 

17.37 


23 

30.00 

5.90 

31.67 


24 

19.00 

4.30 

20.60 


25 

5.20 

2.90 

3.17 


26 

1.70 

1.40 

.57 


19 

.46 

.24 

2.16 


20 

14.00 

3.30 

19.69 


21 

25.00 

4.50 

30.40 


22 

16.00 

2.70 

16.52 


23 

5.50 

1.70 

1.98 


24 

1.40 

.51 

.35 


18 

5.80 

2.00 

21.65 


19 

16.00 

3.30 

29.02 


20 

9.10 

1.90 

12.66 


21 

4.60 

1.70 

1.30 


22 

2.20 

.93 

.20 


15 

2.00 

1.10 

3.02 


16 

24.00 

6.10 

23.14 


17 

17.00 

3.80 

27.18 


18 

8.30 

2.50 

9.26 


19 

6.10 

1.40 

.86 


20 

.99 

.30 

.11 


13 

1.20 

.76 

3.57 


14 

12.00 

5.30 

24.04 


15 

27.00 

7.00 

26.01 


16 

7.70 

3.10 

6.44 


17 

4.80 

.73 

.56 


10 

.48 

.55 

.16 


11 

1.50 

.67 

4.16 


12 

14.00 

6.00 

24.66 


13 

13.00 

4.20 

24.42 


14 

4.00 

2.70 

4.23 


15 

2.80 

.98 

.35 


9 

1.20 

.88 

1.89 


10 

5.60 

1.80 

24.83 


11 

11.00 

3.70 

22.55 


12 

1.60 

1.70 

2.63 


13 

1.90 

2.00 

.21 

14.886 


From reference 2 1 












Table 7. Chi Square Analysis of Fragmentation Cross Sections of Iron Projectile Ions on 

Targets Other Than Carbon 

[See page 45 for footnotes] 

Chi square per 

Fragment Fragment Experiment cross Experiment NUCFRG2 degree of 
charge mass number section, mb uncertainty, mb cross section, mb freedom ( n ) 


1 880 A Me V 


141.00 

98.00 

88.00 

75.00 

67.00 

64.00 

56.00 

55.00 

38.00 

56.00 

57.00 

57.00 

50.00 


187.66 

127.29 

104.06 

89.86 

79.78 
72.48 
66.76 
61.97 

57.78 
54.16 
50.95 
47.82 
45.25 


156.00 

11 1.00 
88.00 

83.00 

77.00 

68.00 

65.00 

54.00 

54.00 

63.00 

57.00 

75.00 

50.00 


186.66 

126.30 

103.18 

88.97 
79.10 

71.98 
66.03 

61.34 

57.35 
53.52 
50.33 
47.23 
44.68 


174.04 

127.60 

91.05 

84.12 

73.41 

68.92 

52.89 

52.72 

45.24 

52.27 

43.47 

58.21 


210.00 

137.75 

113.32 

97.92 

87.22 

79.57 

73.26 

68.26 
64.01 

60.22 
56.78 

53.58 






Table 7. Continued 








Chi square per 


Fragment 

Fragment 

Experiment cross 

Experiment 

NUCFRG2 

degree of 

Target 

charge 

mass number 

section, mb 

uncertainty, mb 

cross section, mb 

freedom (n) 

27 A1 

1569 AMeV b 

13 


45.37 

1.82 

50.85 



12 


51.76 

1.94 

48.17 



11 


45.23 

1.81 

46.10 



10 


49.11 

1.88 

43.75 

2.354 

32 S 

1880 AMeV a 

25 


250.00 

22.00 

219.02 



24 


128.00 

16.00 

140.59 



23 


86.00 

16.00 

115.73 



22 


64.00 

10.00 

100.47 



21 


91.00 

13.00 

89.77 



20 


97.00 

14.00 

81.95 



19 


55.00 

21.00 

75.66 



18 


74.00 

13.00 

70.46 



17 


66.00 

14.00 

66.35 



16 


74.00 

12.00 

62.24 



15 


50.00 

8.00 

58.95 



14 


106.00 

14.00 

55.66 



13 


78.00 

18.00 

52.98 

7.376 

63 Cu 

1569 AMeV b 

25 


238.96 

6.78 

265.32 



24 


147.44 

3.73 

158.74 



23 


98.89 

3.00 

132.79 



22 


98.45 

2.97 

116.63 



21 


73.64 

2.57 

105.32 



20 


80.32 

2.67 

97.14 



19 


59.98 

2.31 

90.85 



18 


61.18 

2.32 

85.56 



17 


49.41 

2.09 

81.17 



16 


59.58 

2.27 

77.52 



15 


49.82 

2.08 

74.04 



14 


72.20 

2.48 

70.87 



13 


51.47 

2.10 

68.46 



12 


61.03 

2.27 

66.09 



11 


50.17 

2.06 

64.08 



10 


54.55 

2.14 

62.03 

4.867 

63 Cu 

1880 A Me V a 

25 


219.00 

20.00 

268.30 



24 


149.00 

16.00 

158.27 



23 


121.00 

15.00 

132.58 



22 


101.00 

14.00 

1 16.20 



21 


100.00 

15.00 

105.34 



20 


98.00 

14.00 

97.25 



13 









Table 7. Continued 


Target 


Chi square per 

Fragment Fragment Experiment cross Experiment NUCFRG2 degree of 
charge mass number section, mb uncertainty, mb cross section, mb freedom ( n ) 


1 880 A Me V 


U/ Ag 

1880 A Me V 1 


88.00 

95.00 

86.00 

56.00 

88.00 

72.00 

179.00 


280.00 

218.00 

117.00 

124.00 

104.00 

1 18.00 

79.00 

84.00 

79.00 

96.00 

64.00 

158.00 

112.00 


56.00 

457.00 

206.00 

150.00 

152.00 

129.00 

107.00 

111.00 

100.00 

101.00 

109.00 

133.00 

81.00 


14.00 

15.00 

13.00 

11.00 

15.00 

11.00 

27.00 


23.00 

21.00 

15.00 

16.00 

13.00 

14.00 

11.00 

14.00 

14.00 

13.00 

13.00 

20.00 

19.00 


82.00 

34.00 

22.00 

19.00 

19.00 

18.00 

17.00 
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ATOMIC ELECTRON CORRELATION AND PARTICLE AND ANTI-PARTICLE 
INDUCED SINGLE- AND DOUBLE-IONIZATION 
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Physics Department, Southern Illinois University, Carbondale, Illinois, 62901, U.S.A. 
^Current Address: Shawnee Community College, Ullin, Illinois, 62992, U.S.A. 


1. INTRODUCTION 

In a series of experiments, Andersen et al. [1,2] have found that the ratio, R, of double- 
ionization of He to single-ionization of He, by incident proton and anti-proton differs substantially 
in the incident energy range of 0.3 to 40 MeV. In this energy range one expects the impulse 
approximation to be valid and that the ionization probabilities are well described by Lewis- 
Merzbacher (noted hereforth as LM) Theory [3] which predicts this ratio to be the same in both 
cases. According to this theory, the interaction V (LM) causing the transition is the potential 
between the incident projectile and two bound electrons and is given by (minus and plus signs are, 
respectively, for incident proton and anti-proton) 


V(LM) = ±— ^ — 



(1) 


where R, r,, and r 2 are coordinates of the projectile, electron one and electron two, respectively. 
In Bom approximation, the scattering amplitude f LM (0) is simply proportional to the square of the 
matrix-element of (1) 


/ LU ( 0 )~ | ///exp (/(*,«, -*,>>, 


)}i|f/r, ,r 2 ) ViLM^Jr^dr^dR \ 2 


( 2 ) 


In (2) i| ifs v r 2 ) and i|/ 0 (r,r 2 ) are final and initial electronic wave functions, respectively. The 
incident and outgoing projectile wave functions are represented by plane waves of wave numbers 
kji o and k x n v respectively, n 0 and /?, being appropriate unit vectors. Since the scattering 
amplitude is an absolute square of the matrix-element of V(LM), incident proton and anti-proton 
should yield the same cross section for single— as well as double-ionization. 



Calculations of Reading and Ford [4-6] using forced impulse approximation use (1) as the 
basic interaction causing the transition but include correlation in electronic wave functions. Then 
calculated R agrees to some extent with the data for incident proton but falls short 

observation for incident anti-proton. 

Das. in his thesis [7] has pointed out that the two electrons in the final channel are far apart 
both for single- and double-ionization and hence do not interact anymore. Thus, the interaction 
causing the transition is not given by (1) but the interaction between two electrons must be added 
to (1). In section 2 we derive the expression for scattering amplitudes starting from the 
Schroedinger equation of the system and then do systematic approximation to examine this point. 
We find that electronic interact, on must be added to (1). The consequences for the inclusion of 
this term in the calculations are that single- and double-ionization of He by p and p should, in 
principle, be different. This has been observed in the recent experiments of Hvelplund et al. [8] 
and Andersen et al. [9], In section 3, we apply the theory to both single- and double-ionization 
at a few incident energies and find reasonable agreement with the data. Section 4 summarizes our 


conclusion. 

2. THEORY 

The Schroedinger equation for He-atom and an incident projectile in the center of mass 
system of incident projectile and He-nucleus is given by 


{v]+v\)-(*Ze 2 /R)-Ze 2 /r { -Ze 2 /r 2 +e 2 lr n * 

2p * 2m 






-E] iKr p r 2 ,/?)=0 


(3) 

In (3) n and m are reduced masses of projectile and electrons, respectively, r,, r 2 , and R are, 
respectively, coordinates of electron one, two and projectile. Minus and plus sign refer, 
respectively, to incident proton and anti-proton. One may expand t in the following complete 

orthonormal set: 

Wv r V R >T.n F n( R Wn( r v r 2> (4) 

where i|/ n (r,, r-,) is defined by the following Hermitian equation. 


2 



( 5 ) 


[ - il<7? *v]y— W. 

2m r. r. 


with (v|; n ,v|/ m ) = 6 n ..» 

Equations for F n (R) may now be obtained by taking scaler products of (3) with t|t 


(*R+ k n- 


Ze 


il)F n (R)=~i^r l Jr 2 ty' n (r v r 2 ) [ = 
R h 2 


Ze 2 

\R-r } 


Ze 


I *-r. 


“] E„,„ 


12 


(6) 

One may formally now write down the expression for differential cross sect.on for 
transition form an initial state to a state n . 

do _ = 4ttV K | J f ^ r id r 2 JRFjR )^ n ( r v r 2 ) (7) 

dQ h 2 k o 

HJr,,r 2 ,R)T.,„.. F ^JV^ 


with 


HJ'VS) = T 


Ze 2 

Ir.-m 


Ze 2 

k 2 -R| 



(8) 


In (7) F no (R) is the solution of the homogeneous pari of (6) with outgoing asymptotic 
condition. (7) is exact and clearly the e 2 /r, : should be included in the calculation. 

In Born approximation, incident and outgoing projectile are represented by plane waves 
which is equivalent to neglecting Ze 2 /R term on the left side of (6). We make further the 
following approximation 

E,„ WwWC wvJ (9a) 


where F 0 (R) is the incident projectile wave function and is correlated ground state 
electronic wave function in the incident channel. In case F 0 (R) is approximated by plane wave 
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the incident channel with wave number differential cross section is given by 

Ze 2 Ze‘ 


2 e 2 


do_ = 47tV^| | | \ dr jr 2 JRexp[i(k/t 0 -kfi ). + 

ctl h 2 k a • " I* r i' |if r i' 


]^( r r r 2)l‘ 


12 


( 10 ) 

The important thing in (8) and (10) is to note that the interaction responsible for the 
transition is not (1) but e 2 /, 12 term must be added to it. This conclusion does not depend on 

approximations (9a) and (9b) but is correct for the most general case given by (8). It is the 
consequence of the necessary asymptotic condition that in the final channel the asymptotic 
behavior of the wave function is given by a product of scattering amplitude, outgoing plane or 
coulomb wave function, intrinsic wave functions of ejected electron and He + for single- 
ionization (or intrinsic wave functions of two ejected electrons and wave functions and He ++ 
for double-ionization). 

Another important consequence of (8) is to note that one expects anti-proton induced 
single- and double-ionization probabilities to be greater than those by proton, in case matrix- 
elements have similar phases. The structure of (10) also implies that single- and double- 
ionization cross sections by anti-proton and proton are, in principle, different from each other, 
the extent of which depends on magnitude of matrix-elements of e lr n relative to those of 

other two interaction terms. 

The Lewis-Merzbacher approximation is obtained by (a) neglecting the e h 12 term 

in the interaction, (b) using products of hydrogenic type of wave functions for 
ty n (r v r 2 ) and ty g (r v r 2 \ and (c) representing incident and outgoing p + or p wave function 

by plane waves. 

3. APPLICATIONS 

We present here calculations at a few incident energies for single- and double- 
ionization by p + and p and the ratio R. The details of the calculation along with results for the 
entire energy range of 0.4 MeV to 4.0 MeV are presented in ref. [10]. 

3.1 Single-Ionization 

In Fig. 1 we have plotted calculations done in (LM) approximation, along with the data 
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of Rudd et at. [9] for incident proton. In LM approximation, the initial ground 

of Rudd 11 of wo hydrogemc lype wave 

function i^Oy.) is uncorrelated and is taken to o p 

functions, each with an effective charge of 1.6875. The final suite electronic wave funcnon ,s 
a product of a bound hydrogenic wave function with charge 2.0 and a ^ 

wave function in Sommerfeld representation with an effective charge • 

using an effective charge is that the ejected electron, being far away, sees a nuclear e g 
shielded by the other electron. Although the calculated results are close to the ata. 

some significant differences between them at a few tens of keV. and correlated 

In Fin 2 we have plotted calculations that include the e r 12 term in 
electronic ground state wave function, for which we have adopted the following variational 

wave function due to Hylleraas [12]. 

wi h (ID 

v^(r,,r 2 ) =Ne ~(y/a o )0 x +r 2 )( 1 +cr w 

01, N and a Q are the normalization constant and, Bohr radius, respectively V and c are 
determined from variational principle and found robe respectively, l.» and 0.14 • The 
variational ground state energy of -78.28 eV is very close to the observe va ue o . 

The final electronic wave function is taken to be the same as the one in the LM approximation. 
Whereas, there is no discernable difference between these calculations and those done in the 
L M approximation at high energies, the calculated results at low energies differ. Calculations 
done in this approximation are in better agreement with the dam compared to those done in the 

LM model. 

in the approximation that include the term and use (11) as ground state wave 
function, the s,ngle-,on,za.,on cross sections for incident proton and an„-proton are different 
which is the case experimentally [8, 9. 15]. The difference is significant both at a very ow 
and high energies. In Table 1 we have compared our calculations with the observed data or a 

few energies. There is a general agreement. 
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Single Ionization 



io" id 1 ur io 3 io 4 10 5 


Energy (keV) 


Fig. 1 Single-ionization cross sections calculated in LM approximation (solid dots) are 
compared with the data of Rudd et al. [ 1 1 1. 


Single Ionization 
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Energy (keV) 


Fig. 2 Single-ionization cross sections calculated using ( 10) and (1 1) (solid dots) are 
compared with the data of Rudd et al. [ 1 1 1. 
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Table 1 . Calculated single-ionization cross sections by p (column 2) and p (column 4) are 
compared with the data of ref. [15] for p + (column 3) and of refs. [8, 9] for p (column 5). 


Energy 

keV 

rjfn + ') in 10’ lo cm‘ 

o(d ) in 10' 16 cm“ 

th 

H 

th 

mm 

10 

0.082 


. 0.106 


20 

0.30 

0.20 

0.32 

0.41 

100 

0.85 

0.85 

0.88 

0.65 

500 

0.40 


0.48 

0.38 

1000 

0.24 

0.23 

0.26 


50000 

0.0065 


0.0293 



In the passing one may note that matrix-elements of e~/r l2 calculated by replacing 
ejected electron wave function by plane wave do not differ significantly from those calculated 
with Sommerfeld Coulomb function but the calculations simplify significantly. 

3.2 Double-Ionization 

In the LM approximation, double-ionization cross section for the case of incident 
proton and anti-proton should be the same but the measurements indicate them to be different 
[8, 9] and hence we are not presenting any calculation of double- ionization in the LM 
approximation. 

We present here calculations using (10) for both cases. Because of the e /r 12 term in 
the interaction, we expect the double-ionization probabilities for proton and anti-proton to be 
different. The initial ground stage i|r 4 (r,,r 2 ) is again represented by (11) and the final 

electronic state y„(r x s 2 ) is taken to be a product of two Sommerfeld Coulomb function with 

charge Z=2, except for the matrix-elements of e*/r 12 , for which a product of two outgoing 
plane wave function is used. The latter is justified on the ground that actual calculations done 
in the case of single-ionization with Coulomb and plane wave show no significant difference in 
numerical values of matrix-elements of this term. 
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Fig. 3 Calculated double-ionization cross sections for incident proton (solid dots) and 
anti-proton (solid diamonds) and experimental data for incident proton [13,14] 
and anti-proton [8,9) are plotted as a function of energy. 
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Fig. 4 Calculated ratio for incident proton (solid dots) and anti-proton (solid diamonds) 
are compared to respective experimental observations for proton (open squares 
1 1,2 1 and triangles [13] and anti-proton (open circles [1,2]). 
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Data for proton induced double-ionization taken by Shah and Gilbody [13] and Puckett 
and Martin [14], along with our calculations, are shown in Fig. 3. The agreement is 
satisfactory. In the same figure we also present our calculations for anti-proton induced 
double-ionization and compare them with the measurements of Hvelplund et al. [8, 9] and the 
agreement is good. Calculated cross sections for incident p are larger than those for proton, 

as anticipated. 

3.3 Ratio 

In Fig. 4, we have plotted the ratio of double- to single-ionization cross section, R, for 
incident proton and anti-proton for a few tens of keV and a few tens of MeV. Our calculations 
both for incident proton and anti-proton are in agreement with the data of Shah and Gilbody 
[1 1] and Andersen et al. [1,2]. The key factor responsible for the difference in two cases is 
the inclusion of e : /r 12 term in the interaction which makes double-ionization probabilities for 
the anti-proton substantially larger than those for the proton. 

4. CONCLUSION 

One may conclude that (a) the difference in measured R values for incident proton and 
anti-proton can in general be accounted for within the standard non-relativistic atomic physics, 
(b) In principle, one is to include two-electronic interaction in the calculation for both single- 
and double-ionization, (c) Single-ionization probabilities for incident proton and anti-proton 
are in principle different. This difference is pronounced at high energy and (d) double- 
ionization probabilities by anti-proton are substantially larger than those by protons from a few 
hundred keV to a few tens of MeV incident energies. 

Lastly at low energies one expects the cross sections for the p case to exhibit some 
structures originating from the quasi bound states in the attrafctive anti-proton-He potential 
[16]. These should be absent for incident proton because of the repulsive nature of proton-He 
potential. There are some indications of this in the data of [8, 9]. 
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On the ionization of He by proton and anti-proton: II. Double ionization and ratio 
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Abstract We have calculated double ionization cross section of He by proton and anti-proton in 
the incident energy range of 50 keV to 20.0 MeV using a coupled channel approach restricted to 
two channels, incident and final. The projectiles have been represented by plane waves. The 
ground state wave function of He is taken to be a variational correlated one, and the final state 
wave functions of two ejected electron is taken to be Coulomb functions for calculating matrix 
elements of the electron-projectile interaction, and plane waves for calculating matrix elements of 
electron-electron interaction. The calculation satisfactorily accounts for the data on double 
ionization cross section and on the ratio of double to single ionization cross sections. 
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I. Introduction 

In a series of experiments Andersen et al (1986, 1987) have found that the ratio R, of double 
to single ionization of He by proton and anti-proton differ significantly in the incident energy 
range of a few tens of keV to a few tens of MeV. Subsequently, measured single as well as 
double ionization cross sections of He by proton (Puckett and Martin 1970, Shah and Gilbody 
1985, 1982a, b; Shah et al. 1989) are found to differ significantly from those by anti-proton 
(Andersen et al 1989, 1990 and Hvelplund et al 1994). 

These findings are important because, within the framework of theories that represent 
projectiles by plane waves in the incident and outgoing channels and considering only electron- 
projectile interaction as perturbation, these ionization cross sections for incident proton and anti- 
proton and the ratio should be the same. 

A number of theoretical approaches have been used to understand these phenomena. The 
relevant theories for single ionization are summarized in Das and Malik (1996) (This paper is 
termed as I). 

Meng et al’s (1993) calculations of single and double ionization cross sections of He for 500 
keV incident proton and anti-proton using classical trajectory Monte Carlo method (Olson 1987) 
that includes electron-electron interaction approximately (noted henceforth as dCTMC) do not 
adequately explain the data Reading and Ford (1987, 1988 and 1994) have calculated the 
double-ionization cross sections by incident proton and anti-proton only at three energies namely 
at 300, 1000 and 6000 keV and estimated them at 500, 3000 and 10000 keV using a time 
dependent approach termed as forced impulse approximation or FIM. Their calculated results for 
the 300, 500 and 1000 keV incident proton energies are in good agreement with the data of Shah 
and Gilbody (1994) taken at 326, 500 and 1000 keV/amu. However, their calculated cross 
sections for the 300 and 500 keV incident anti-proton differ about 10 to 20% from the data taken 
at 270, 326 and 503 keV incident energies (Andersen et al 1986, 1987 and Hvelplund et al 1994). 
No calculation has been reported for lower incident energies in the FIM. For a proper 
understanding of the underlying reaction mechanism associated with the ionization process, it is 
important to understand the low energy data at least in the 100 to 300 keV region. 

Das (1994) in his thesis has calculated double ionization cross sections for incident p + and p‘ in 
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the energy range of 100 keV to 20 MeV using coupled channel approach but restricting it to two 
channels only. A few of these calculations at low and high energies have so far been reported 
(Das and Malik 1995) In this paper we present calculated results for the double ionization and 
the ratio in the energy range of a few tens of keV to 20 MeV and compare them with the 
observed data for both incident proton and anti-proton. The theoretical ratios have been obtained 
by using the calculated single ionization cross sections within the frame work of the same theory 
in I. The outline of the theory, which is discussed in details in I, is presented in the following 
section. Results and discussion, conclusions, and acknowledgment are presented in subsequent 
sections 

2. Theory’ 

In the coupled channel theory the equation governing the motion of proton or anti-proton 
incident on He, is given by (Mott and Massey 1 965 and I) 

= -/ dr x dr ilj',J r v r 2) H J R Sx T 2 ) X,„ <F n 'W,,, <r v r 2 ) (1) 

In (1) p is the reduced mass of the scattered particle with respect to He nucleus and R, r, and r 2 
are, respectively, radius vectors of the projectile, electron 1 and electron 2 in the center of mass 
system of incident projectile and He nucleus F n (R) is the wave function of the projectile in 
channel n. Minus and plus sign refer, respectively, to incident proton and anti-proton. k„ is the 
wave number of the nth channel and related to total energy of the system E and the non- 
interacting electron energies E^ by fi 2 k n :! =2 l u(E-E nm ). ijr^r,, r 2 ) in the above expression is a 
product of two hydrogen like wave functions in the states n and m including continuum states and 
is eigenfunctions of the following equation: 
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( 2 ) 


- v; »' vl) • - — - — .r,) =£„*_(»•, ^ 2 ) 

2 M 


r \ l 2 


with 


(4/ ,tlf ■ )=b /6 

V T /WJ’ T ^ //j / >j 


(3) 


The interaction energy in ( 1 ) is given by 


H J R Sv r i> + 


Ze 


r r R\ 


Ze 2 

\r 2 ~R\ 



(4) 


For single as well as double ionization, both electrons in the final channel are far apart from 
each other and represented by hydrogen-like eigenfunction of (2). 

One may obtain the expression for differential cross section by imposing the asymptotic 
condition 


HkJi-T\lnR) 

HmFjW =/„( 0 . 4 >)^-^ 

1*1— 


(5) 


where r|=,uZe 2 /}ik^. 

In (5) f n (0,4>) is the scattering amplitude at angles (0,40 defined with respect to the direction 
of incident momentum of the projectile The expression for differential cross section is given by 
(Mott and Massey 1965) 


do 


= l/„(0,<P)l 2 


_ 4tt 2 |! 2 k n 


If / / dr x drJR Ko( R W nm (r v r i) H J r v r 2 R )L„ < r i- r 2 ) 


( 6 ) 


Total cross section, o n , is then calculated by integrating over the solid angle dQ. F no (R) in (6) 
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is the solution of the homogeneous part of(l) with outgoing boundary condition (5). 

The expression (6) is exact but represents an infinite set of coupled equations and 
approximations are to be made for numerical evaluation. Since the incident projectile energy of 
10 keV to 20 MeV is considerably higher than the double ionization energy of He, it is reasonable 
to consider only two channels and use Born approximation in which the projectile in the incident 
and outgoing channels is described by plane wave. In this two-channel approximation, the 
summation over n' is replaced by one term denoted by ‘o’. Hence 




( 7 ) 


In (7) i|/ g (r 1 , r 2 ) is fully correlated wave function of the ground state of He. The differential 
cross section in this two channel plane wave Born approximation, denoted as 2cPWBA, is given 
by 


da _ 47i 2 p 2 L- • 
dQ ’ 


f dr x dr 2 dR exp (ik/t 0 -ik/t n )^ a (r v r 2 ) H mt {r v r 2 R)ty g (r v r 2 ) 
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In (8) and are, respectively, unit vectors in the directions of incident and outgoing wave 
numbers k„ and Iq,. 

It is important to note that the e 2 /r,, term is to be included in H^, in principle This conclusion 
does not depend on the two-channel approximation and the approximation (7) but is correct for 
the most general expression (6) and is a consequence of the boundary condition that two electrons 
in the final channel do not interact asymptotically and each is represented by a one-body wave 
function of hydrogenic type in the final channel. This term does not cause any transition in the 
absence of projectile-electron interaction because in that case, the projectile wave functions are 
zero and hence the matrix elements of the term e 2 /r 12 are zero. In the single ionization case, the 
contribution of the matrix elements of the e 2 /r 12 term to the total cross section is not significant for 
incident energies greater than a few hundred keV. 

The He-ground state wave function i|i g is taken to be the following (Hylleraas 1930) 
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^ g (r v r 2 )=Nexp(-y(r i +r 2 )/a o )( 1 +ct 12 ) 


(9) 


where N is the normalization constant and a 0 , the Bohr radius Parameters y and c have been 
determined from the Raleigh-Ritz variation principle to be 1.69 and 0.142, respectively. The 
variational ground state energy of -78.28 eV is very close to the observed value of -78.62 eV 
The electronic wave function in the final channel, i|r n (r,, r 2 ), is taken to be a product of two 
continuum Coulomb wave functions having two units of charge in the Sommerfeld representation 
(Bethe 1930; Lewis and Merzbacher 1958; Sommerfeld 1931) for calculating matrix elements of 
electron-projectile interaction and a product of two outgoing plane waves for calculating matrix 
elements of e 2 /r 12 . The latter approximation is motivated by the fact that in the case of single 
ionization, matrix elements of the e 2 / r 12 term calculated with the plane wave and Coulomb wave 
functions for the ejected electron yield about the same numerical value (Das and Malik 1996) in 
the entire energy range. 


3. Results and Discussion 

(a) Proton induced double ionization and the ratio 

In Table 1 and Fig. 1 we have presented our calculated results for double ionization cross 
section in the incident proton energy range of about 50 keV to 20 MeV and compared them with 
the measurements of Shah and Gilbody (1985), Shah et al (1989) and Puckett and Martin (1970). 
The measurements extend from about 16 keV to 2400 keV, although the data below 50 keV have 
not been plotted The calculation can account for the energy dependence and magnitude of 
observed cross section from 80 keV to 2400 keV very well. Below 80 keV the calculated values 
for the cross section seem to be slightly higher than the data which, however, have large errors. 
The observed cross section seem to be rather flat from 40 to 100 keV incident energy but not the 
calculated values 

Meng et al (1993) have calculated double ionization cross section in the dCTMC method at 
two energies. Their calculated values are 7.3 x 10‘ 19 cm 2 and 4.1 x 10‘ 19 cm 2 for 300 and 500 keV 


6 



incident energies, respectively, and are considerably larger than our calculated values and 
observed data. Calculations and estimates of Ford and Reading for the 300, $00, 1000 and 3000 
keV incident energies are about 10% lower than our results and slightly lower than the obsen/ed 
values but their calculations for the 6.0 and 10.0 MeV agree with our results. 

In Table 3 and Fig. 3 we have presented our calculated and the observed ratios of double to 
single ionization cross sections for incident proton using the values of single ionization calculation 
done in I in the 2cDWA, in the energy range of 50 keV to 20 MeV. The calculation reasonably 
accounts for the measured ratio in the entire energy range. 

(b) Anti-proton induced double ionization and the ratio 

We present the calculated double ionization cross section by incident anti-proton along with 
the observed data in the incident energy range of 40 keV to 50 MeV in Table 3 and have plotted 
them in Fig 3 The calculation does account for the data in this large energy range and the 
observed fact that double ionization cross section by incident proton systematically lie lower than 
that by incident anti-proton in the energy range of 100 to 500 keV. The calculation predicts this 
trend to continue right up to 50 MeV This is the consequence of the fact the three terms in (4) 
have the same sign for the anti-proton case but not for the proton case. 

The ratio of double to single ionization cross section for incident anti-proton is calculated 
using the results of the single ionization cross section calculated in I in the 2cDWA. The 
calculated ratio is presented and compared to experimental observation in Table j and Fig. j. The 
data cover a very large energy range from a few tens of keV to 20 MeV The theoretical 
calculation can satisfactorily account for the data essentially in the entire energy range. 

In I, the possibility of resonances or the formation of virtual (p - He) system has been raised. 
Such a physical situation would also influence double ionization cross section of He by anti- 
proton and hence the ratio. This process is estimated in I to occur at a few tens of keV. Kimura 
et al (1994) have incorporated such processes in their calculation of the ratio of incident energies 
below 50 keV and found them to be important. This process is absent for incident proton for 
which, on the other hand, the charge transfer channel becomes available. Kimura et al s 
calculations seem to explain the observed difference in R, for the two cases at energies below 50 
keV Incorporation of such phenomena is beyond the scope of our investigation which is. 
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therefore, restricted to incident energy grater than about 80 keV. At very low incident energy, 
Fermi-Teller effect (1947) might influence significantly the double-ionization cross section of He 

by anti-proton 

Satisfactory explanation of the observed data for double and single ionizations and their ratios 
for the incident energies above 80 keV indicates that the coupled channel method is a suitable one 
for describing ionization process. Because the perturbed potential (4) is different for double- 
ionization cross section for the two cases, theoretically, one predicts that the cross section for 
incident anti-proton will be higher than that by incident proton at incident energies greater than 
those used in current experiments. It would be interesting to verify this. 

Theoretically, ionization cross sections by electron impact should have the same perturbation 
(4) for anti-proton, but has the additional complication of incorporating the Pauli principle which 
is important at lower incident energy (Malik and Trefftz 1961, Trefftz 1963). For He target, the 
effect of the Pauli principle on the cross section may not be critical for incident electron energies 
greater than 1 0 keV and the situation is then similar to that for incident anti-proton, except for the 
difference in masses. Hence, the energy dependence of the ratio should be qualitatively the same 
in both cases which is indeed the case, experimentally. 

Within the framework of the 2cPWBA, alpha-induced ionization cross section should be 
similar to that by incident proton, except that the masses of two projectiles are different. This 
difference would cause the key features of the energy dependence of the cross sections such as 
maxima to occur at a different energy in the two cases but the general pattern should be the same, 
which is the case experimentally 

4. Conclusion 

The observed single and double cross sections and their ratio for incident proton and anti- 
proton can be well accounted for in the energy range of a few tens of keV to a few tens of MeV 
by the coupled channel method. It is important however to incorporate proper asymptotic 
condition in the incident as well as final channel. 
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Table Caption 


Table 1 


Table 2 


Table 3 


Comparison between the calculated double ionization cross section of He by 
proton, o (this work) and the selected observed ones, o (Shad et al) and o (puckett 
et al). Shah et al refers to (Shah and Gilbody 1985, 1982a, 1982b) and Puckett et 
al refers to (Puckett and Martin 1970). 

Comparison between the calculated double ionization cross section of He by anti- 
proton, o (this work) and the selected observed data noted as o (Aarhus), o 
(Aarhus) refers to data in (Andersen et al 1987a, 1987b) and (Hvelplund et al 
1994). 

Comparison between the calculated ratio of double to single ionization of He by 
proton and anti-proton noted in columns 2 and 5, respectively, and the 
corresponding data for incident proton and anti-proton noted in columns j and 6, 
respectively Aarhus group refers to the selected date of (Andersen et al 1987a, 
1987b, 1989a, 1989b) and (Hvelplund et al 1994). Shah et al refers to the selected 
data of (Shah and Gilbody 1985) and (Shah et al 1989). 


12 



Figure Caption 


Fig l Calculated double ionization cross section of He by proton noted as this work (p ) 

is plotted as a function of proton incident energy and compared to the selected 
data References to the data are given in the caption of Table 1 

Fig 2 Calculated double ionization cross section of He by anti-proton noted as this work 

(p ) is compared to the selected data. The references to the data are noted in the 
caption of Table 2. 

Fig. 3 Calculated ratio of double to single ionization cross sections for incident proton, 

noted as this work (p ) and incident anti-proton, noted as this work (p ) are 
compared to the corresponding incident proton data, as noted Shah et al (p ) and 
incident anti-proton data, noted as Aarhus group (p ) References to these selected 
data are noted in the caption of Table 3. 
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Table I : 


Double Ionization cross section by proton. 


Energy 

a (This work) 

a (Shah et al.) 

o (Puckett et al 

(keV) 

( 10 18 cm-) 

( 10 18 cm 2 ) 

(10' ,8 cm 2 ) 

50.00 

1.49580 



58.46 


0.920 ± 0.240 


60.00 

1.25340 



70.00 

.1*. 07890 



70.56 


0.9 10 ± 0. 190 


80.00 

0.94707 



80.67 


1.057 ± 0. 101 


84.67 


0.870 ± 0.180 


90.00 

0.84403 



100.00 

0.86345 



100.79 


0.890 ± 0 . 180 * 


150.00 



0.970 

200.00 

0.65856 


0.678 

201.59 


0.560 ± 0.017 


300.00 

0.35672 


0.356 

322.54 


0.315 ± 0.022 


400.00 

0.19341 


0.242 

403. 18 


0.236 ± 0.012 


500.00 

0.15492 


0.176 

503.97 


0 . 1 72 ± 0.008 


600.00 

0.12921 


0.146 

645.08 


0.127 ± 0.001 


700.00 

0.1 1082 


0.1 14 

800.00 

0.09701 


0.0977 

806.35 


0.096 ± 0.0051 


900.00 

0.086264 


0.0828 

1000.00 

0.075660 


0.0732 

1007.90 


0.0745 ± 0.0027 


1612.70 


0.0454 ± 0.0037 




0.0365 * 0.035 
0.0300 ± 0.02 1 


2000.00 

2015.90 

2398.90 

3000.00 

4000.00 

5000.00 

6000.00 

7000.00 

8000.00 

9000.00 

10000.00 

30000.00 

40000.00 

50000.00 


0.037883 


0.025934 

0.019455 

0.015566 

0.012973 

0.01 1 121 

0.0097309 

0.008650 1 

0.0077853 

0.0025956 

0.0019468 

0.0015574 



Table 2: 


Double Ionization cross section by anti-proton. 


Energy o (This work) a (Aarhus data) 


(keV) (10 l8 cm 


50.00 

3.5223 

53.90 

60.00 

2.9594 

67.10 


70.00 

2.4612 

80.00 

2.2260 

80.40 

90.00 

1.9555 

91.20 

100.00 

1.8366 

101.60 

200.00 

0.87739 

227.90 

270.10 

300.00 

0.68717 

400.00 

0.49127 

403.60 

500.00 

0.36346 

503.60 

600.00 

0.29580 

700.00 

0.25284 

800.00 

0.22134 

900.00 

0.19682 

1000.00 

0.17019 

2000.00 

0.088013 

3000.00 

0.059670 

4000.00 

0.045388 

5000.00 

0.036515 

6000.00 

0.029799 


) (10' lx cm-) 

2.01 1 ± 0.205 
1.923 ± 0.156 

1.543 ± 0.139 

1.657 + 0.145 

1.577 + 0.123 

0.720 + 0.039 
0.696 ± 0.035 

0.447 + 0.023 
0 . 4 10 ± 0.021 


0.025372 


7000.00 



8000.00 

9000.00 

10000.00 

30000.00 

40000.00 


0.023202 

0.019736 

0.017763 

0.0059220 

0.0045416 


50000.00 


0.0035534 



Table 3: Ratio of double to single ionization cross sections. 


Energy 

(keV) 

This work 

(P + ) 

(x 10 - 3 ) 

Aarhus group. 

(P + ) 

(X 10 ' 3 ) 

Shah et al. 

(P + ) 

(x10"3) 

This Work 
(P~) 

(xlO'3) 

Aarhus group 
(P-) 

(x 1 0*3) 

50.00 

20.2300 





67 10 





28.70 ± 2.30 

70.00 




28.5910 


80.67 



12.70 ± 1.20 



100.00 

10.1039 * 



21 .0229 


102.00 





23.10 ± 1.10 

200.00 




13.2250 


227.90 





14.46 ± 0.78 

300.00 

6.3632 



11.7990 


326.30 





13.00 ± 0.66 

400.00 

5.2785 



10.9873 


403.18 

5.2072 


5.36 ± 0.28 

10.5294 


403.60 





11.32 ± 0.58 

500.00 

4.7762 



9.9972 

10.00 ± 0.60 

503.97 

4.7153 


4.65 ± 0.27 

8.5462 


667.00 





8.50 ± 0.30 

700.00 

3.5231 



7.5163 


800.00 

3.4275 



7.2449 


806.35 

3.4144 


3.60 ± 0.20 

7.2074 


900.00 

3.2343 



6.9946 

7.15 ± 0.25 

1000.00 

3.1818 

3 . 22 - 0.10 


6.5322 


1007.90 

3.1627 


3.29 ± 0.16 

6.5243 


1200.00 

2.9642 

3.12 ± 0.20 


6.2312 


1370.00 





6 . 40 ± 0.19 

1612.70 

2.9263 


2.86 ± 0.21 

6.0021 


1900.00 

2.8859 

2.80 ± 0.20 


5.8296 

5.80 ± 0.25 

2000.00 

2.8419 

2.76 ± 0.06 


5.6404 


2015.90 

2.8326 


2.82 ± 0.26 

5.6212 


2100.00 

2.8201 



5.4923 



2260 00 


5.10 ± 0.41 



2398 90 

2.8189 


2.66 ±0.18 

5.3642 


2000 00 

2.7759 

2.64 ± 0.14 


5.1364 


3200 00 

2.6925 

2.50 ±0.20 


5.0044 

5.00 ±0.25 

* * 0J 00 

2.6473 

2.50 ±0.1 5 


4 7049 

4.50 ±0.25 

•uoo 00 

2.6001 

2.50 ±0.10 


4.5488 


7 000 00 

2.5230 

2.51 ±007 


3.8582 


' ( >000 00 

2.5148 

2.49 ± 0.10 


3 3233 

4 10 ±0 41 

20000 00 

2.4498 



2.3021 

2.90 ±0 20 



Double Ionization 



Energy (keV) 



Double Ionization 



Energy (keV) 



Ratio 



C IT, C l/“) O o 

m <n n — — 


(^oi) a 


Energy (keV) 



On the ionization of He by proton and anti-proton: I. Single ionization 


Tirthanath Das* and F. Bary Malik 

Physics Department, Southern Illinois University, Carbondale, IL 62901 

*Presently at division of mathematics and science, Shawnee Community College, Ullin, 
IL 62922, U.S.A. 

Abstract. We have presented here calculation of single ionization cross section of He by 
incident proton and anti-proton within the framework of a coupled channel theory 
restricted to two channels. The incident energy range considered is from 10 keV to 50 
MeV. Projectile wave functions in initial and final channel are taken to be plane waves. 
The initial ground state wave function of He is taken to be a correlated variational one 
and the final channel electronic wave function is represented as an antisymmetrized 
product of a hydrogen-like bound (Is) and continuum Coulomb functions. The calculated 
results are in reasonable agreement with the observed data and explain the observed 
difference in cross sections for single ionization by incident proton and anti-proton. 

These results have been used in calculating the ratio of double to single ionization cross 
section by proton and anti-proton in a subsequent paper. 



1. Introduction 

In a series of experiments Andersen et al (1986, 1987) have revealed that the ratio, R, 
of double to single ionization cross sections of He by incident proton, p* differs 
substantially from that by incident anti-proton p\ Subsequent measurements of anti- 
proton induced single and double ionization cross sections of He by Andersen et al 
(1990) and Hvelplund et al (1994) indicate these to differ from those measured for 
incident p" (Shah and Gilbody 1981, 1985 and Shah et al 1989). 

Before the series of investigations of single ionization by proton impact done by 
Shah and Gilbody using crossed-beam coincidence technique, there have been many 
previous measurements of proton induced single ionization of He using condenser plate 
technique which are referred to in the 1985 article of Shah and Gilbody. In particular, 
Rudd et al (1983) have measured single ionization cross section of He by incident protons 
in the energy range of a few keV to 4000 keV. Their measurements are in general 
agreement with those of Shah and Gilbody (1985) for energies above 200 keV. It is, 
therefore, well established that single ionization cross section of He is different for 
incident p~ and p\ These findings are significant because, within the framework of 
Lewis-Merzbacher’s theory (henceforth denoted as LM) using plane wave Bom 
approximation (Lewis and Merzbacher 1956) and considering only projectile-electron 
interaction as perturbation one expects single ionization cross section of He by p + and p‘ 
and ratio, R, to be the same. 

In Coulomb-Bom approximation, used by Trefftz (Trefftz 1963) to study ionization 
of OV and OV1 by electron impact, the wave function of projectile in the final channel 
for the single ionization case should be represented by Coulomb function. In that case 
one expects, in principle, a difference in ionization cross section by p* and p' because of 
the difference in signature of their charges in Coulomb function in the final channel. 
Similarly, in a distorted wave approximation used by Malik and Trefftz (1961) to study 
single ionization of O V by electron impact, the single ionization cross section of He by p* 
and p‘ should, in principle, differ. In fact, the study of Fainstein et al (1987) using 
distortion of electron wave function due to incident projectile in eikonal approximation 
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and Coulomb wave function for ejected electron as well as projectile in the final channel 
(termed henceforth as CDW-EISi Continuum-distorted-wave-eikonal-initial-state model) 
indicates different single ionization cross sections for incident proton and anti-proton, 
because of the difference of signature for charge of the projectile wave function in the 
final channel in the two cases. Their calculations for incident anti-proton reasonably 
account for the data from about 15 keV to 3000 keV (Hvelplund et al 1994). Their 
calculation for incident proton could reproduce the general trend of the observed cross 
section as a function of incident proton energy (Shah et al 1985, 1984) but differs by as 
much as 25% quantitatively (Andersen et al 1990). Unfortunately, no calculation of 
double ionization cross section of He in this approximation has been reported. 

The theoretical calculation based on classical-trajectory-Monte Carlo (termed as 
CTMC) method (Olson 1987, Shultz 1989) does not provide satisfactory agreement with 
the data (Hvelplund et al 1990) for proton induced single ionization of He. Calculations 
of single ionization cross section, based on an improved version of this model 
(Montemayor and Sciwietz 1989 and Meng et al 1993) that includes correlation between 
two electrons in terms of screening potentials (the model is termed as dCTMC), give 
0.404 and 0.293 A 2 for 300 and 500 keV incident proton, respectively, and 0.327 and 
0.25 1 A 2 for incident anti-proton, respectively. The observed data for incident proton 
energies of 260, 320 and 500 keV/amu are, respectively, 0.587±0.011, 0.516±0.09 and 
0.370±0.015 A 2 (Shah and Gilbody 1985) and for incident anti-proton energies of 270. 1, 
326.3 and 503.6 keV are, respectively, 0.471±0.08, 0.437±0.08 and 0.354±0.08 A 2 
(Hvelplund et al 1994). Calculated double-ionization cross sections in this method 
overestimate significantly the proton induced double ionization data at these energies and 
slightly under estimate the corresponding anti-proton data. 

Ford and Reading (1987a, b, 1989, 1994) have approached the problem from time- 
dependent view point. Their model, termed as forced-impulse method or FIM, breaks up 
the time development in small segments. Within a given segment, the electrons are 
treated in an uncorrelated fashion but between the segments the correlation is considered. 
Their earlier calculations done for incident energies higher than a few hundred ke V 
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including s and p electron orbitals are in qualitative agreement with the proton data but 
falls short on the ratio of double to single ionization by p'. However, their latest 
calculation done only at 300, 500 and 1000 keV incident energies (Ford and Reading 
1994) estimating the contributions of 9d orbitals to both single and double ionization are 
in reasonable quantitative agreement with the data around these energies. For example, 
their estimated cross sections for single ionization for incident proton are 0.498, 0.346 
and 0.203 A 2 at incident energies of 300, 500 and 1,000 keV, respectively, which are 
close to the observed values of 0.5 16±0. 009, 0.370±0.015 and 0.226±0.008 A 2 measured 
(Shah and Gilbody 1985) at incident energies of 320, 500 and 1000 keV/amu, 
respectively. Similarly, they estimate single ionization cross sections by incident anti- 
proton of energies 300, 500 and 1,000 keV to be 0.458, 0.331 and 0.188 A 2 , respectively. 
These are close to observed values (Hvelplund et al 1994, Andersen et al, 1990) of 
0.471±0.02, 0.437±0.02, 0.354±0.02 and 0.195±0.01 A 2 at incident energies of 270.1, 
326.3, 503.6 and 1 130 keV, respectively. So far, calculations at other energies, 
particularly lower ones have not been reported, although the measurements are available. 

The early calculations of Bell and Kingston (1969) using plane wave Bom 
approximation that approximates the electron-projectile interaction as dipole or the 
corresponding velocity operator (so called Bethe (1930) approximation) and uses a 
correlated ground state wave function do not agree with the data. In particular, the 
maximum of the cross section is not accounted for (Shah and Gilbody 1985). Of course, 
the anti-proton induced single ionization cross section is equal to that by proton in this 
approximation, and the observed difference in data in the two cases cannot be explained. 

Toburen et al (1978) have calculated single ionization cross section at 300, 500 and 
1000 keV incident proton energies in the LM approximation, except replacing initial 
hydrogenic He ground state wave function by Herman-Skillman’s (1963) parameterized 
Hartree-Fock wave function and the Sommerfeld Coulomb wave function for ejected 
electron by a distorted wave (Manson et al 1975). Their calculated result agrees with the 
data of Rudd et al (1983) at 300 keV but is slightly higher than the data at 500 and 1000 
keV, and differ significantly from the 100 keV data. Lewis and Merzbacher (1958) have 
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done calculation in the LM approximation up to about 200 keV incident proton energies, 
with limited agreement with the data. 

Das and Malik (1994, 1995) have analyzed single and double ionization of He by 
proton and anti-proton using coupled channel approach (Mott and Massey 1965 or Malik 
and Trefftz 1961) restricting it to two channels, initial and final, and representing 
projectiles in both channels by plane waves. Their calculations for a few very low and 
very high energies are in very good agreement with the data for single and double 
ionization. In this paper we extend the calculation from 10 keV to 50 MeV for proton 
and anti-proton induced single ionization of Helium and compare them with the available 
data. We also present calculation in the LM approximation in this energy range. In a 
subsequent paper calculations for double ionization cross section and the ratio R for 
incident p" and p' will be presented for incident energies ranging from a few tens of keV 
to 20 MeV. 

We present the theory in the following section and discuss the results in section 3. 
The conclusion and acknowledgment are presented in sections 4 and 5. 


2. Theory 

The Schroedinger equation for an incident projectile and a He-atom in the center of 
mass system of the incident particle and the nucleus of He may be written as 


— V 2 - — (V: ♦ V;) - ( + Ze 2 /R) - Ze 2 /r? -Ze 2 /r. >e 2 lr.-.— — 

2m R 2m 1 - 1 12 |rt-r,| 




- = o 


( 1 ) 


In (1) p and m are reduced masses of projectile and electrons with respect to He-nucleus, 
respectively, r,, r 2 and R are coordinates of electrons 1 and 2 and the projectile, 
respectively. Minus and plus sign refer, respectively, to incident proton and anti-proton. 
In (1) r, 2 = |r, - r 2 |, the interelectronic distance. 

One may expand i(f (R, r,, r 2 ) in the following orthonormal set 

Wr v rJ = XF n (R)Zj nm (r v r 2 ) (2) 
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In view of the fact that the two electrons are far apart in the final channel and do not 
interact, we define the orthonormal set to be the eigenfunctions of the following hermitian 
Hamiltonian: 


% 


2m 


r, 


Ze 


'P„J r v r 2) = E JnS r v r 2 ) 


(3) 


The above Hamiltonian represents accurately the asymptotic situation in the final 
channel where the ejected electron is far away from the He + ion and hence e‘/r 12 term is 
negligible. 

One may obtain the following equation for F n (R) by taking a scalar product of (1) 
with respect to iji^, and integrating over the coordinates r, and r 2; 


I ^ HE ~ E - )±? r FSR) 




= - 1 d ^ dr ii„$,J r v r 2 y 


Ze- 


Ze- 


I Rr, 


I R~r. 


\ln',f n i R )'LnA,,,'A r x r 2) 


(4) 


12 


The orthogonality condition on iji^ has been used in deriving (4). Using 



2m 


(£ - E J 


(5) 


one gets 


<Vj * *’ * ir )F ” (R) 

2m 


—f dr l dr zLM r v r 2> 


Ze 


Ze- 


\R-rJ | R-r 


21 r !2 


►L '.fn i R )L„ 4* ir v r 2 ) (6) 
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The equation (6) represents an infinite set of coupled channel equations and suitable 
approximations are needed before actual calculations could be earned out. 

The above equation is now to be solved with the following boundary condition in the 

final channel 


-r\tnR) 

lim F n (R) -/„(e,<P) £ — 

l*b- 


(7) 


with q = pZeVhk,,, Z being the projectile charge. f n (0,(p) is the scattering amplitude in the 
channel n at angles (0,(p) defined with respect to the direction of incident momentum. 
Following Mott and Massey (1965), one may obtain the expression of the differential 
cross section 


do 

JQ 


[/„M ! 


4rc 2 |i 2 k n 


^-fffdr t dr 2 dRFJR)l,^Jr v r 2 ) 


H,JR.r v r 1 )L^A R ».AVi)\ : 


( 8 ) 


with 


H J R d V r 2 ) 


Ze 2 

|r, - R \ 


Ze 2 + e^_ 
\r 2 ~ R \ r u 


(9) 


k Q in (8) is the incident channel wave number and F’ h (R) is the solution of the 
homogeneous part of (6) compatible with asymptotical boundary condition (7), i.e., 
outgoing Coulomb wave function for p + or p'. 

The total single ionization cross section, o + , is obtained by integrating over the solid 
angle dQ. 

It is important to note that the expression (8) is exact and the e 2 /r, 2 term is to be 
included in calculating matrix elements. In case electron-projectile interaction is 
switched off, the projectile wave function is zero and hence matrix elements of e 2 /r 12 and 
cross section are zero, i.e., there is no transition, as it should be. 


7 



Since the energy interval considered in this work is significantly larger than the 
ionization energy, we consider explicitly only two channels, the initial incident and the 
final. In this two-channel approximation the summation over n in (2) has two terms, the 
initial channel denoted by F U (R) and final channel noted F n (R). The summation over n' 
in (8) reduces to the following term: 




( 10 ) 


i|f g in ( 1 0) is the correlated ground state wave function of He. 

In this approximation the differential cross section is given by 


Ja_ = 4nY_ r d dr dRF 'jR^ (r r ) 

I j , 2 nhy P 27 

n o 

VJKr v r 2 )F o (Rn<r v r 2 )\ 2 


(ID 


In case F^ (R) and F U (R) are taken, respectively, to be Coulomb and plane waves this 
approximation is equivalent to Coulomb Bom approximation (CBA) used by Trefftz 
(1963) for single ionization of oxygen ions by electron. 

The two channel plane wave Bom approximation (2cPWBA) is obtained by 
replacing F 0 and F^ with plane wave function in respective channel. The differential 
cross section in the 2cPWBA is given by 

— = — — — -| f dr.dr 2 JRexdi(k n -k n)R] 

dd h * k o l J 1 2 n\<>o n ) \ (12) 

Total cross section is obtained by integrating over df). In (12) ft 0 and ft are unit vectors in 
the direction of momenta in incident and final channels, respectively. 

The LM approximation is obtained by (i) replacing tjt g with a product of two 
hydrogen-like (Is) wave functions, and (ii) neglecting e 2 /r 12 term in H int (R, r,, r 2 ). 

The Bethe approximation or sometime denoted as zero-order or lowest-order Bom 
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approximation is obtained by (i) omitting e : /r 12 from Hj n „ (ii) replacing t|/ g in (12) with a 
product of two (Is) hydrogen-like wave functions and (iii) keeping only terms up to dipole 
in the expansion of electron-projectile interaction. 

Because H int (R, r„ r 2 ) given by (9) is different for incident proton and anti-proton, 
the calculated single ionization cross section in the 2cPWBA for the two cases differ, in 
principle, from each other. On the other hand, in the LM and Bethe approximations, they 
should be the same. 

In the next section we present calculations in the LM approximation and the 
2cPWBA. 

Calculations in the LM approximation are done taking electronic wave function in 
the initial state to be a product of two hydrogenic (Is) functions, each having an effective 
charge 1.6875 and final state electronic wave function to be a product of a hydrogenic (Is) 
wave function with effective charge 2.0 and a continuum Coulomb wave function in 
Sommerfeld representation (Wentzel, 1929; Bethe, 1930; Sommerfeld, 1931) with an 
effective charge of 1.09. 

The initial state wave function used in the 2cPWBA is the following correlated wave 
function due to Hylleraas (1930): 

'\f g ( r l S 2 )=Nexp(-(y/a 0 )(r i +r 2 ))( 1 +cr? 2 ) (13) 

In (13) N and a 0 are, respectively, the normalization constant and the Bohr radius, y 
and c are determined from the Raleigh-Ritz variational principle to be 1.69 and 0. 142, 
respectively. The variational ground state energy is -78.28eV which is very close to the 
observed value of -78.62 eV. The finale state wave function in the 2cPWBA is taken to 
be the same as the one in the LM approximation. 
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3. Results and Discussion 

(a) The LM approximation 

In table 1 we present our calculated results in the LM approximation in the incident 
energy range of 10 keV to 50 MeV for both incident proton and anti-proton. In the 
energy range of 10 keV to 200 keV, our calculations are in agreement with those of Lewis 
and Merzbacher (1958) which also serve as a check to our numerical code. In Fig. 1 the 
data of Shah and Gilbody (1981, 1982, 1985 and 1989) and Rudd et al. (1983) for the 
single ionization cross section of He by proton have been presented as a function of 
incident energies and compared to our calculation in the LM approximation. The theory 
can account for the data satisfactorily for incident energies greater than a few hundred of 
keV. At lower energies there is room for improvement. 

(b) 2cPWBA for incident proton 

In table 1 we have also presented our calculation for single ionization cross section of 
He by incident protons in the energy range of 10 keV to 50 MeV in the two channel plane 
wave Bom approximation. In general the calculated cross sections in the 2cPWBA are 
lower than those in the LM approximation. The difference is about 10% below 50 keV 
and above 30 MeV and only a few percent at other energies. In Fig. 2, we have compared 
our calculation with the data and the agreement is excellent in the entire energy range. In 
particular, the agreement between the theory and experiment has improved significantly 
at lower energies compared to those calculated in the LM approximation. 

(c) 2cPWBA for incident anti-proton 

In table 2 we have presented our calculation in the 2cPWBA for incident anti-proton 
in the energy range of 10 keV and 50 MeV. They are compared to the observed single 
ionization cross section in Fig. 3. The calculation reproduces the general energy 
dependence of the cross section and the data quite well. However, the calculated values 
between 80 keV and 300 keV are somewhat higher than the data. Calculated cross 
section for incident anti-proton is always higher than those for proton but the difference is 
insignificant in the energy range of 70 to 2000 ke V. 
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(d) General discussion 

For the proton induced single ionization of He, our calculation can account for the 
observed cross section in the entire energy. The agreement between the data and our 
calculation is slightly better than the one between the data and the calculation of Fainstein 
et al. (1987) particularly below 150 keV and above 1000 keV. The Bethe approximation 
used by Bell and Kingston is inadequate at lower energies. The calculated values in the 
2cPWBA agree with the data somewhat better than the calculations in the dCTMC (Meng 
et al 1993) and in the FIM approximations (Ford and Reading 1994). The LM 
approximation used by us is very similar to the method of Toburen et al (1975, 1978), 
except for the use of wave functions for initial electronic states, and the numerical 
numbers in the two cases are very close. 

The contribution from the matrix elements of e 2 /r, 2 term to cross section for the 
proton is only significant at low energies, i.e., below 100 keV. In Table 3, we have 
compared calculations in the 2cPWBA done using plane wave for the ejected electron 
instead of Coulomb wave function to compute matrix elements of the e 2 /r 12 term and they 
do not differ significantly. It is, therefore, sufficient to evaluate its matrix element using 
plane wave for the ejected electron, instead of Sommerfeld representation of Coulomb 
wave function. The use of plane wave for the ejected electron simplifies the calculation 
considerably. 

For incident anti-proton the observed data agree very well with our calculation for 
incident energies greater than about 200 keV. At lower energies, the calculation of 
Fainstein et al (1987) are in somewhat better agreement with the data compared to ours. 
The data are in better agreement with our calculation compared to those done in the FIM 
(Ford and Reading 1994) and in dCTMC (Meng et al 1993). 

The 2cPWBA theory predicts that at very high incident energies i.e., energies above a 
few Me V to non-relativistic energy the single ionization cross section of He by anti- 
proton should be higher than that by incident proton. This is the consequence of the 
difference in the perturbed potential, (9), in two cases. The data are available up to about 
2.5 MeV for incident proton (Shah et al 1989) and 3.0 MeV for incident anti-proton 
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(Andersen et al 1990) and seem to borne out this theoretical calculations. It would 
certainly be most interesting to have data at higher incident energies. 

At incident energies lower than about 100 keV, our calculated cross section foi 
incident anti-proton lie a little higher than those for protons, whereas the data indicate the 
situation to be the opposite. At this incident energy region, the difference in Coulomb 
distortion in the two cases is likely to play an important role. In addition, the possibility 
of resonances for the anti-proton-Helium system may have to be considered. In the first 
approximation the location of such resonances is determined by the eigen-energies of the 
homogeneous part of (4) (Malik 1992). Clearly, there should be no resonance for the 
incident proton case since the (p~ - He) system is unbound. On the other hand, the (p - 
He) system has bound states, the energies of which are given by E = -(20.0 keV)Z 2 /n : in 
the hydrogen-like approximation. Hence, the resonances could contribute to the cross 
section in the region of a few tens of ke V for the anti-proton case and the spacing of these 
resonances should be from about a tenth to a few keV. It is interesting to note that the 
measured single-ionization cross section for the incident proton in the energy range of a 
few tens of keV is a smooth function of energy (Shah et al 1989) which is expected in the 
absence of resonances but the same for incident anti-proton seems to exhibit some 
structure (Hvelplund et al 1994) as expected in the presence of resonances. It would be 
most interesting to establish whether such structures are actually present for the anti- 
proton case. 

4. Conclusion 

The 2cPWBA can reasonably account for the energy dependence of the observed 
single ionization cross section of He by proton as well as anti-proton in the incident 
energy range of 10 keV to 3 MeV. The calculation has been presented up to 50 MeV 
incident energy with a view of using these results to calculate the ratio of double to single 
ionization in a subsequent paper. The data on the ratio are available in the MeV region. 
The inclusion of e : /r 12 term is not critical at energies higher than a few hundred keV and 
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lower than a few MeV but important at lower and higher energies. The success of the 
LM approximation, 2cPWBA and CDW-EIS in reproducing the data reasonably implies 
that the Coulomb distortion for the ejected electron must be included in calculation. 
Inclusion of ground state correlation in the initial electronic state of He is important at 
energies below a few hundred keV. 
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Table Caption 


Table 1 Comparison of calculated single ionization cross section of He^by proton in 
the LM approximation (column 2) and the 2cDWBA (column j) with the 
selected observed ones in experiments of Shah and Gilbody (1981, 1985), 
Shah et al (1989) and Puckett et al (1970). 

Table 2 Comparison of calculated single ionization cross section of He by anti- 
proton in the 2cDWBA with the selected observed data of the Aarhus group 
(Andersen et al, 1990 and Hvelplund et al 1994). 

Table 3 Comparison between calculations of single ionization cross section of He 
by p~ done in the 2cPWBA (column 2) and the approximation where the 
plane wave function is used for ejected electron in evaluating the matrix 
element of e 2 /r 12 (column 1). 
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Figure Caption 


Fig i Comparison between the calculated single ionization cross section of He by 

proton done in the LM approximation and selected data. Shah et al refer to 
data reported in (Shah and Gilbody 1981, 1985 and Shah et al 1989). Rudd 
et al and Puckett et al refer to data from (Rudd et al 1983) and (Puckett and 
Martin 1970), respectively. 

Fig. 2 Comparison between the calculated single ionization cross section of He by 

proton done in the 2cDWBA with selected data referred to in the caption of 
Fig. 1. 

Fig. 3 Comparison between the calculated single ionization cross section of He by 

anti-proton done in the 2cDWBA and selected data of the group primarily 
based at Aarhus. Aarhus group data refer to (Andersen et al 1990) and 
(Hvelplund et al 1994). 
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Table 1: 

Single ionization 

cross section by proton. 


Energy 

o (LM) 

a (This work) 

a (Shah et al.) 

a (Puckett e 

(keV) 

( 10' K> cm-) 

O 

O 

3 

to 

( 10 _ 16 cm 2 ) 

( 10' 16 cm 

10.00 

0.087826 

0.08098 



l 1 .09 



0.067 ± 0.003 


19.15 



0.192 ±0.006 


20.00 

0.31223 

0.29868 



28.22 

• 


0.340 ±0.014 


30.00 

0.51242 

0.49882 



50.00 

0.75246 

0.73939 



64.51 



0.750 ± 0.023 


70.00 

0.84746 

0.83285 



100.00 

0.85942 

0.84991 



100.79 



0.845 ± 0.046 


131.03 



0.809 ±0.012 


150.00 




0.880 

200.00 

0.70243 

0.69203 


0.718 

200.59 



0.693 ±0.01 7 


300.00 

0.57224 

0.55834 


0.510 

400.00 




0.429 

403.18 



0.441 ±0.009 


500.00 

0.41026 

0.40160 


0.350 

503.97 



0.370 ±0.015 


600.00 




0.314 

700.00 

0.32354 

0.31343 


0.273 

900.00 




0.227 

1000.00 

0.24879 

0.23698 


0.207 

1000.79 



0.226 ± 0.008 


2000.00 

0.14394 

0.13287 



2000.79 



0.1295 ±0.0011 


2398.90 



0.1 125 ±0.0013 


5000.00 

0.061243 

0.058680 



7000.00 

0.044078 

0.042975 





9000.00 

10000.00 

30000.00 

40000.00 

50000.00 


0 . 03974 1 

0.032416 

0.016347 

0.0099929 

0.0072450 


0.033956 
0.030742 
0.010692 
0.0080732 
0.006485 1 



Table 2: Single Ionization cross section by anti-proton. 


Energy 

o (This work) 

o( Aarhus data) 

(keV) 

(10' 16 cm 2 ) 

( 10 16 cm 2 ) 

10.00 

0.10570 


12.90 


0.199 ±0.037 

20.00 

0.32337 


24.20 


0.473 ± 0.038 

30.00 

0.52567 


34.60 


0.59 1 ± 0.027 

50.00 

0.76732 


53.90 


0.741 ±0.030 

70.00 

0.86083 


100.00 

0.87632 


101.60 


0.645 ± 0.026 

194.40 


0.528 ± 0.020 

200.00 

0.71770 


270.10 


0.471 ±0.019 

300.00 

0.58235 


326.30 


0.437 ±0.017 

500.00 

0.42529 


645.00 


0.3 10 ±0.022 

700.00 

0.33639 


1000.00 

0.26054 


1130.00 


0.195 ±0.014 

2000.00 

0.15604 


2470.00 


0.105 ±0.012 

2918.00 


0.099 ±0.012 

3000.00 

0.11617 


5000.00 

0.081557 


7000.00 

0.065762 


9000.00 

0.056724 


10000.00 

0.053450 


30000.00 

0.033468 


40000.00 

0.030741 


50000.00 

0.029239 




Table 3: 


Energy 

(keV^ 

10.00 

20.00 

30.00 

50.00 

70.00 

100.00 

200.00 

300.00 

500.00 

700.00 

1000.00 

2000.00 

3000.00 

5000.00 

7000.00 

9000.00 

10000.00 

30000.00 

40000.00 


o (plane wave) o (Sommerfeld) 

(10 16 cm-) (10 16 cm 2 ) 


0.08195 

0.30163 

0.50312 

0.74468 

0.83808 

0.85457 

0.69506 

0.56060 

0.40255 

0.31455 

0.23779 

0.13330 

0.093427 

0.058862 

0.043107 

0.034059 

0.030835 

0.010723 

0.0080959 

0.0065040 


0.08098 

0.29868 

0.49882 

0.73939 

0.83285 

0.84991 

0.69203 

0.55834 

0.40160 

0.31343 

0.23698 

0.13287 

0.093131 

0.058680 

0.042975 

0.033956 

0.030742 

0.010692 

0.0080732 

0.0064851 


50000.00 
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